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LASER  INTERFEROMETER  FOR  THE  STUDY  OF  PIEZOELECTRIC  AND 

ELECTROSTRICTIVE  STRAINS 

Q.M.  Zhang,  W.Y.  Pan,  and  L.E.  Cross 
Materials  Research  Laboratory 
The  Pennsylvania  State  University 
University  Park,  PA  16802 

ABSTRACT 

A  modified  Michelson  interferometer  is  used  to  study  the  strain  properties  of  piezoelectric 
and  electrostrictive  materials.  For  small  displacement,  a  feedback  loop  is  introduced  to  stabilize  the 
system  against  the  low  frequency  optical  ^path-length  drifting  and  the  system  is  capable  of  resolving 
displacements  of  the  order  otf&^CTor  the  strain  induced  by  domain  switching  a  dual-channel 
signal  detection  scheme  is  used  which  automatically  reads  out  the  displacement  of  the  sample.  The 
effect  on  the  measurement  of  the  sample  bonding  to  a  substrate  and  other  related  problems  are 
discussed. 
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I. 


INTRODUCTION 


The  study  of  the  piezoelectric  and  electrostrictive  properties  of  ferroelectric  and  other 
materials  involves  the  detection  of  the  strain  induced  by  the  external  driving  forces.  The 
displacement  related  to  this  land  strain  varies  from  the  order  of  10  A  to  more  than  1  pm.  Various 
methods  have  been  used  to  measure  these  strains  precisely.  Among  diem,  optical  methods  show 
great  promise  because  they  do  not  require  calibration  in  the  length  scale,  do  not  have  mechanical 
contact  and  also  can  be  used  to  probe  the  displacement  profile  of  the  sample  surface  at  various 
locations  and  under  different  conditions.1"5  The  basic  principle  of  the  laser  interferometer  is  that, 
for  a  monochromatic  light  of  wave  length  X  interfering  with  a  reference  beam  the  interference  light 
intensity  at  a  detection  point  is 

I-Ip  +  Ir  +  2IpIrcos(4iiAd/X)  (1) 

where  Ip  and  I,  are  the  light  intensity  for  the  probing  beam,  where  the  sample  is  placed,  and  the 
reference  beam  respectively,  Ad  is  the  optical  path-length  difference  between  the  two  beams.  A 
relative  change  Ad  between  the  two  beams  results  in  a  change  in  the  interference  light  intensity 
which  can  be  detected  by  different  methods  depending  on  the  magnitude  of  the  displacement 
change  Ad.  For  the  convenience,  eq.  (1)  can  be  rewritten  as: 

I  -  l/Zdaax  +  W  +  ltfdmw  *  W  cosf4xAd/X)  (2) 

where  Im»T  and  Im;n  are  the  maximum  and  minimum  interference  light  intensities,  which  can  be 
measured  from  the  interference  intensity  curve  conveniently. 

The  interferometer  which  is  discussed  below  makes  use  of  two  different  signal 

detection  schemes.  In  measuring  displacement  smaller  than  300  A,  a  dc  feedback  loop  is  used  to 


stabilized  the  system  against  thermal  and  mechanical  optical  path-length  fluctuation  and  it  is  able  to 
resolve  very  small  displacement  of  the  order  10*3  A  by  using  an  ac  signal  detection  method.  For 
the  large  displacement  measurement,  a  dual  channel  signal  detection  scheme  is  used  which  is  able 
to  automatically  and  precisely  determine  the  turning  point  of  the  motion  of  tire  detected  surface  and 
trace  out  the  displacement  The  error  associated  with  the  measurement  of  the  displacement  of 
several  pm  is  on  the  order  oflO  A,  which  is  a  negligibly  small  of  the  total  displacement  measured 
II.  DETECTION  OF  SMALL  DISPLACEMENT 

The  schematic  drawing  of  the  interferometer  as  well  as  the  electronic  system  for  detecting 
small  displacements  is  shown  in  Fig.  1.  In  this  simple  Michelson  interferometer,  the  incident  laser 
beam  (X  *  6328  A)  is  split  into  half-power  beams  by  a  cube  beamsplitter.  After  bouncing  back 
from  the  sample  surface  and  the  reference  mirror,  the  probing  beam  and  the  reference  beam  are 
recombined  by  the  same  beamsplitter  and  form  an  interference  pattern  at  the  detection  point  which 
is  described  by  eq.  (1).  In  the  real  experiment  situation,  due  to  the  finite  size  of  the  laser  beam,  if 
the  optical  path-length  difference  Ad  is  not  small  enough,  the  interference  pattern  will  show  a 
distribution  in  the  light  intensity  at  the  detection  plane  such  as  that  shown  in  Fig.  2.  If  the  photo¬ 
detector  aperture  is  larger  than  the  central  fringe  diameter,  the  sensitivity  of  the  system  can  be 
greatly  reduced.  Since  this  central  fringe  diameter  is  inversely  proportional  to  Ad,  effort  has  been 
made  in  die  system  set-up  process  to  reduced  this  Ad  so  that  in  the  detection  screen,  the  view  spot 
shows  uniform  tight  intensity  distribution,  indicating  that  the  central  fringe  diameter  is  much  larger 
than  the  aperture  of  the  photo-detector.  In  order  to  further  improve  the  sensitivity,  a  condenser 
lens  is  inserted  in  front  of  the  photo-detector  to  further  expand  the  diameter  of  the  central  fringe. 
Under  these  conditions,  the  tight  intensity  seen  by  the  detector  is  that  described  by  eq.  (1). 

For  small  displacement  measurement,  it  is  desirable  to  set  the  interference  between  the  two 
beams  at  a  point  about  (2n  +  l)X/8  (x/2  point)  where  the  light  intensity  change  will  be  maximized 
for  the  same  displacement  change  in  Ad.  Near  this  point,  we  can  write  Ad  as 
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Ad*dlc  +  (2n+ 1)V8  (3) 

and  eq.  (1)  can  be  reduced  to 

I-Ip  +  4±2IpIfSin4jid10/X*Ip  +  Ir±2IpIr4xd1,A  (4) 

In  eq.  (4),  the  sign  of  "+"  or  depends  on  whether  n  is  an  even  or  odd  number.  The 
approximation  of  Sinx  *  x  is  valid  for  small  x  and  in  our  case,  the  error  associated  is  less  than  1% 
if  /d^/  is  less  than  130  A.  Eq.  (4)  indicates  that  the  light  intensity  change  is  linearly  proportional  to 

the  change  d  for  small  displacement.  This  makes  it  possible  to  use  ac  method  to  determine  the 
strain  response  of  the  sample  to  the  external  driving  fence. 

The  stabilization  of  the  interferometex  at  the  x/2  point  is  achieved  by  introducing  a  feedback 
loop  in  the  reference  arm  as  shown  in  Fig.  1.  As  a  matter  of  fact,  the  configuration  used  here  can 
set  the  system  at  any  reference  point  by  adjusting  the  voltage  level  at  the  differential  amplifier 
reference  input  This  feedback  loop  respondes  only  to  low  frequency  signals  (<  10  HZ).  It  does 
not  affect  the  measurement  at  higher  frequencies.  In  the  same  manner,  it  will  not  reduce  the  high 
frequency  noise,  instead,  this  extra  loop  may  add  more  noise  at  high  frequency.  However,  due  to 
the  fact  that  the  system  can  stay  at  the  x/2  point  for  any  length  of  time,  dK  can  be  determined  by  a 

lock-in  amplifier  which  has  high  noise  rejection  ability  and  the  signal  can  be  averaged  over  any 
time  length  desired.  This  is  the  key  to  improve  the  sensitivity  of  the  system  for  the  detection  of 
small  displacements. 

For  a  sinusoidal  displacement,  d*.  a  c^cos  ox,  the  optical  signal  is  converted  by  the  photo¬ 
diode  to  a  voltage  change  and  then  pre-amplified.  This  signal  is  detected  by  an  Ithaco  lock-in 
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amplifier  (model  393)  as  which  is  a  rms  value  of  the  detected  signal  corresponding  to  the 
displacement  d</>/2.  From  eq.  (4),  it  can  be  shown  that 

d0-(VV2*)(V0ul/Vp.p)  (5) 

Where  VP-P  is  the  peak  to  peak  value  of  the  interference  signal,  which  corresponds  to  the  change  in 
the  interference  signal  (Im«  *  Imin)- 

III.  PIEZOELECTRIC  AND  ELECTROSTRICTIVE  MEASUREMENT 

Piezoelectric  and  electrostrictive  coefficients  were  measured  for  X-cut  quartz,  hard  PZT,  soft 
PZT,  and  PMN  samples.  In  these  measurements,  the  sample  was  polished  on  one  face  to  optical 
quality  in  order  to  perfectly  reflect  the  incident  laser  beam.  Since  it  was  difficult  to  polish  the 
surfaces  of  some  ceramic  samples  to  optical  flatness,  two  alternative  methods  were  used.  First, 
silver  paint  was  put  on  the  sample  surface  and  dried  at  600°.  The  silver  paint  film  was  polished  to 
a  minor  surface  which  also  serves  as  the  electrode.  In  the  second  method,  a  small  minor  was 
attached  to  the  sample  for  a  reflection  surface.  However,  this  latter  method  is  not  suitable  for  the 
measurements  at  high  frequency.  In  the  experimental  process,  it  was  found  that  the  quality  of  the 
reflection  surface  of  die  sample  is  very  important  for  accuracy.  A  poor  reflection  surface  can  cause 
distortion  of  the  interference  pattern  and  introduce  severe  errors  in  the  determination  of  small 
displacements  using  eq.  (3). 

Sample  bonding  is  another  important  factor  which  can  affect  the  strain  measurement  of  the 
sample  to  a  great  extent  In  our  measurement,  the  sample  size  may  vary  from  the  order  of  mm  in 
linear  dimension  to  1  cm.  The  sample  was  elec  traded  by  gold  sputtering.  The  back  face  of  the 
sample  was  attached  to  the  optical  sample  holder  using  silver  epoxy.  The  sample  holder  is  able  to 
move  in  X-Y-Z  directions  and  rotate  the  sample  about  the  vertical  and  horizontal  axes.  This  is 
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necessary  for  proper  adjustment  of  the  system.  Another  sample  bonding  method  was  proposed  in 
literature.2  This  method  is  not  suitable  for  the  high  frequency  measurements  and  may  introduce 
extra  complications  into  the  system.  In  using  the  bonding  material  to  attach  the  sample  to  the 
holder,  it  is  obvious  that  the  bonding  agent  will  have  a  clamping  effect  on  the  sample  (introducing 
stress  in  the  sample).  It  has  been  shown  that  for  die  same  driving  voltage,  the  strain  in  the  clamped 
sample  can  be  reduced  to  60%  of  that  measured  in  unclamped  samples.3  The  effects  of  the 
bonding  are  discussed  below. 

In  order  to  determine  the  piezoelectric  or  electrostricdve  coefficients  of  a  sample  by  this 
interferometer  detection  scheme,  we  have  assumed  that  the  strain  induced  by  the  driving  voltage 
causes  only  the  front  face  of  the  sample  to  move  (the  surface  facing  the  beam)  and  the  motion  of 
the  back  face  is  negligible.  This  assumption  may  not  be  valid  if  the  bonding  agent  is  not  rigid 
enough  to  hold  the  sample  or  if  the  sample  is  being  driven  at  a  very  high  frequency.  An  ideal 

|  bonding  material  should  have  an  infinite  elastic  constant  in  the  sample  motion  direction  (z  direction 

in  Fig.  3)  and  zero  elastic  constant  in  the  lateral  direction.  A  simple  model  can  be  used  to  analyze 
the  z  direction  motion  of  die  sample  as  shown  in  Fig.  3.  To  avoid  complications,  a  self-consistent 
method  is  employed  to  treat  the  problem.  In  Fig.  3,  the  bonding  agent  acts  as  a  spring  with  an 

i 

|  elastic  constant  k  (k  »  cA/1,  where  c  is  the  elastic  constant  of  the  bonding  agent  at  unit  area  and 

length,  A  the  area  and  1  the  thickness  ).  If  we  assume  that  the  back-face  of  the 

sample  is  perfectly  clamped,  the  motion  of  the  mass  center  of  the  sample  will  be  half  of  the  motion 
of  die  front  face.  For  a  displacement  with  a  frequency  f  and  an  amplitude  z0  of  the  front  face,  the 
force  needed  to  support  the  sample  would  be  equal  to  (2x2mz0f2).  For  a  typical  sample,  m  is  of 

- 

the  order  of  1  gram,  and  the  silver  epoxy  (ACME  3021)  has  c  (elastic  constant)  of  the  order  of  109 
newton/m2,  which  is  about  two  order  of  magnitude  smaller  than  that  of  the  sample. 6  The  contact 
area  of  the  silver  epoxy  in  our  experiment  is  of  the  order  of  1  mm2  and  the  thickness  is  of  the  order 
of  0.1  mm.  The  ratio  between  the  back  face  motion  (due  to  the  deformation  of  the  silver  epoxy  to 
produce  a  force  to  support  the  sample)  and  the  front  face  motion  z/Zq,  therefore,  will  be 


proportional  to  f2xlO’9.  Hence,  the  error  introduced  by  the  back  face  modon  of  the  sample  can  be 
smaller  than  1%  if  the  operating  frequency  of  the  system  is  lower  than  20  K  HZ 

The  clamping  effect  in  the  lateral  direction  is  negligible  due  to  the  fact  that  the  bonding 
material  used  is  about  two  orders  softer  than  the  sample.  Since  the  a  (Poisson  ratio)  is  in  the  region 
0  <  a  <  1/2,  the  influence  of  this  lateral  clamp  on  the  longitudinal  motion  of  the  sample  is  very 
small  as  long  as  the  system  is  far  away  from  a  resonance  frequency  and  the  sample  strain  in  the 
lateral  direction  is  small  (<0.2%). 

The  measurements  on  X-cut  quartz,  hard  PZT  and  soft  PZT  support  the  points  discussed 
above.  At  frequencies  below  2  K  HZ,  the  dn  value  for  quartz  is  2.26  x  10*12  within  1%  error. 

The  frequency  scan  for  the  quartz  sample  is  shown  in  Fig.  4.  The  same  result  can  be  seen  for  hard 
PZT  which  is  shown  in  the  same  figure.  Mechanical  resonance  occurs  at  frequencies  above  5  K 
HZ  This  kind  of  resonance  was  also  observed  in  the  measurement  by  another  group. 7  From  the 
model  in  Fig.  3,  the  mechanical  resonance  due  to  bonding  for  the  system  would  occur  at  frequency 
near  or  above  SO  K  HZ  We  do  not  fully  understand  why  this  resonance  occurs  at  such  a  low 
frequency.  It  may  be  due  to  the  thin  electrode  layer  between  the  silver  epoxy  and  die  sample  or  the 
sample  holding  system  which  has  a  large  mechanical  mass.  For  very  high  frequency  strain 
measurements,  any  mechanical  method  of  bonding  the  sample  is  obviously  inadequate.  Further 
effort  will  be  made  to  separate  the  motion  for  the  sample  as  a  whole  body  and  the  relative  motion 
(the  strain)  between  the  front  and  back  faces  of  the  sample. 

In  Ftg.  4,  we  also  show  the  voltage  scan  of  the  hard  PZT  sample  at  a  driving  frequency  of 
900  HZ  It  can  be  seen  that  the  strain  response  of  the  sample  to  the  external  driving  voltage  is  still 
linear  for  a  field  as  weak  as  5mV/cm  which  will  be  discussed  in  a  separate  paper.  The  resolution 
of  our  system  at  this  frequency  is  on  the  order  of  10*3  A.  Such  kind  of  resoltion  can  only  be 
achieved  when  the  experiment  is  done  under  a  quiet  environment  and  the  signal  is  averaged  in  the 
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lock-in  over  a  relatively  long  time  period.  It  is  clear  that  the  resolution  of  the  displacement  should 
be  frequency  dependent  since  the  noise  spectrum  is  a  function  of  frequency.  At  very  high 
frequencies,  the  only  limit  to  the  resolution  of  the  system  will  be  that  from  electronic  noise.  At 
frequencies  below  100  HZ,  any  outside  mechanical  disturbance  can  make  the  measurement  difficult 
and  the  resolution  of  the  system  is  not  better  than  5  x  10-3  A. 

Electrostrictive  coefficients  were  also  measured  for  a  PMN  ceramic  sample.  The 
electrostricdve  coefficient  is  defined  as 

S  -ME2  (6) 

where  E  is  the  field  strength.  In  this  measurement,  a  voltage  V  is  applied  across  a  sample  with 
thickness  L,  and  E  -  V/L.  Strain  amplitude  SQ  -  dJL,  where  dQ  is  the  displacement  measured  by 

the  lock-in  amplifier  at  2f  mode  (the  driving  voltage  is  used  as  the  reference  signal  which  is  at 
frequency  f).  For  a  driving  voltage  V  -  Vgcosox,  V2»(V02/2Xl-cos2<nc).  Eq.  (6)  can  be  written 

as 


S  -  (MV02/2L2)(l-cos2ax)  (7) 

Using  eq.  (5)  to  find  d^, 

M-2doiyV02  (8) 

For  a  PMN  ceramic  sample,  Mu  is  measured  at  a  frequency  of  250  HZ  and  is  found  to  be  1.03  x 
KH^mVV2.  The  same  PMN  sample  was  measured  by  the  strain  gauge  method.8  Qj  j  obtained  is 
1.9  x  10~2,  therefore  M  »  (EEq)2  x  Q  *  1.02  x  10-16m2/V2  (e  =  8300  is  used).  This  is  in  good 
agreement  with  our  value. 


IV.  DUAL  CHANNEL  DETECTION  AND  THE  DETERMINATION  OF  THE 
STRAIN  DUE  TO  DOMAIN  SWITCHING  IN  FERROELECTRICS 


In  the  measurement  of  strain  induced  by  domain  switching  the  maximum  strain  at  saturation 
field  may  vary  from  less  than  10*3  to  more  than  5  x  10-3.  For  the  sample  used  in  this  measurement 
the  thickness  ranges  from  0.15  mm  to  1  mm.  The  mechanical  displacement,  therefore,  is  on  the 
order  of  iim  or  larger  which  correspond  to  several  wavelength  of  the  laser  light  (X  =  6328  A).  It  is 
clear  that  the  ac  measurement  scheme  described  above  is  not  suitable  here.  The  determination  of 
the  displacement  must  involve  recording  the  interference  pattern  with  time.  However,  in  these 
measurements  the  displacement  will  have  two  opposite  directions.  The  reversal  of  the 
displacement  direction  will  cause  a  corresponding  reversal  of  the  interference  pattern  (shown  in 
Fig.  5)  which  may  coincide  with  the  maximum  or  minimum  point  of  the  interference  pattern.  This 
makes  the  data  treatment  difficult  and  confusing.  A  single  channel  signal  detection  is  not  adequate 
and  it  is  impossible  to  read  out  the  strain  automatically.  A  dual  channel  detection  scheme  which 
will  be  described  below  can  overcome  this  difficulty  and  our  experimental  results  show  that  this 
technique  is  much  more  sensitive  than  conventional  techniques  such  as  the  linear  variable 
differential  transformer  (LVDT)  and  the  strain  gauge. 

A  schematical  drawing  for  this  part  of  the  set-up  is  shown  in  Fig.  6.  In  order  to  have  dual 
channel  detection,  the  polarized  nature  of  the  laser  beam  is  used.  A  polarized  beam  splitter  reflects 
the  light  with  the  polarization  in  the  horizontal  plane  and  transmits  the  light  vertically  polarized. 
The  polarizations  of  the  probing  beam  and  the  reference  beam  are  perpendicular  to  each  other.  In 
both  the  probing  arm  and  the  reference  arm,  a  quarter-wave  plate  is  inserted  which  rotates  the 
polarization  direction  in  both  arms  by  90°  before  the  lights  come  back  to  the  polarized  beam 
splitter.  A  neutral  beam  splitter  divides  the  combined  beam  again  into  two  detection  channels.  In 
one  of  the  detection  path,  a  phase  plate  is  used  to  adjust  the  phase  of  the  light  in  one  polarization 


direction  by  X/4  with  respect  to  the  other  polarization  direction.  The  final  interference  signals 
detected  by  the  two  detectors,  therefore,  have  a  90°  phase  difference  as  shown  in  Fig.  5(d).  It  is 
clear  that  in  an  X-Y  plot,  the  combined  signal  from  these  two  detectors  will  be  &  circle  and  the 
angular  position  of  the  signal  represents  the  displacement  due  to  the  strain.  The  two  directions  of 
the  rotation  of  the  circle  correspond  to  the  forward  and  backward  motions  of  die  detected  surface  of 
the  sample.  The  statements  above  can  be  summarized  as 

tan0(t)  =  Vy(tWx(t)  (9) 

and  the  displacement  versus  time  d(t)  is 

d(t)  =  (A/47i)0(t)  (10) 

where  Vy  and  Vx  are  the  signals  from  the  two  photo-detectors  and  t  is  time.  A  computer  program 
can  be  madf.  to  simultaneously  read  out  the  displacement  versus  time  from  eqs.  (9)  and  (10). 
There  is  no  need  to  determine  the  inflection  point  of  the  displacement.  The  rotation  direction  of 
0(t)  automatically  indicates  the  displacement  directions  once  it  is  been  calibrated. 


The  experiment  was  done  on  PLZT  8/65/35  ceramic  samples  supplied  by  Shanghai  Institute 
of  Ceramics  (China).  These  samples  were  hot-pressed  to  transparency  and  the  average  grain  size 
is  23  iim.  The  sample  preparation  is  the  same  as  that  for  the  small  displacement  measurement  In 
the  measurement  the  fringe  signals  from  the  main  detection  channel  and  the  quadrature  channel  as  a 
function  of  time  were  simultaneously  recorded  by  a  digital  oscilloscope  (Nicolet  model  204A). 
The  triangle  driving  voltage  was  also  recorded.  These  data  are  shown  in  Fig.  7.  In  the 
experimental  process  we  used  a  condenser  lens  to  focus  the  light  on  the  sample  surface  as  shown 
in  Fig.  6.  This  lens  is  essential  for  a  proper  measurement  Since  the  movement  of  the  sample 
surface  is  not  always  strictly  perpendicular  to  the  incident  beam,  this  may  cause  a  change  in  I  p  and 
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also  the  relative  phase  of  the  return  beam.  This  will  result  in  an  error  in  the  data  reading.  In  the 
case  that  the  sample  reflection  surface  is  in  the  focal  plane  of  the  condenser  lens,  the  reflected  beam 
is  approximately  from  the  focal  point  of  that  lens.  The  deviation  of  the  return  beam  from  the 
incident  beam,  after  passing  the  condenser  lens,  is  greatly  reduced  due  to  the  fact  that  all  the  light 
originated  from  one  point  of  the  focal  plane  will  propagate  parallel  to  the  corresponding  optical  axis 
of  that  lens.  In  our  case,  it  is  clear  that  this  optical  axis  is  on  the  line  of  the  incident  laser  beam. 
This  insures  that  Ip  is  approximately  constant  during  the  measurement 

The  final  strain  hysteresis  loop  for  the  PLZT  sample  is  plotted  in  Fig.  8.  Initial 

switching  is  not  observed  because  the  signals  were  recorded  after  a  few  switching  cycles.  We 
repeated  the  experiment  at  frequencies  below  10  HZ  for  the  sample  with  different  bonding 
conditions,  that  is,  with  different  contact  area  and  thickness  of  the  silver  epoxy.  For  the  case  that 
the  silver  epoxy  totally  covers  the  back  face  of  the  sample  (maximum  clamp  in  the  measurement) 
and  the  case  that  the  bonded  area  is  on  die  order  of  1/10  of  the  total  area  of  die  sample,  no  obvious 
difference  in  the  strain  hysteresis  loop  and  the  magnitude  of  the  strain  was  observed  which  agrees 
with  the  discussion  above  regarding  the  clamping  effect  in  our  measurements. 


Several  variables,  such  as  thermal  and  mechanical  drift  of  the  optical  path  length,  may 
affect  the  absolute  accuracy  of  the  measurement3  We  believe,  in  our  experiment  that  the  main 
error  is  introduced  by  the  limited  resolution  of  the  data  recording  system,  the  digital  oscilloscope. 
For  time  resolution  1  ms  under  1  HZ  driving  field  (1/1000  a  cycle  resolution),  the  uncertainty  of 
the  displacement  for  PLZT  sample  is  about  30  A.  Compared  with  the  total  displacement  (--1  pm) 
for  the  sample,  this  error  is  actually  negligible. 

In  order  to  further  improve  the  bonding  condition  to  reduce  possible  clamping  effect  for 
extra  large  strain  measurement  and  to  eliminate  the  sample  reflection  face  tilting  due  to  mechanical 
imbalance  of  the  bonding  surface,  a  vacuum  sample  holder  was  designed  and  used  in  some  of  the 
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experiments.  It  is  shown  schematically  in  Fig.  9,  the  front  face  of  the  holder  is  polished  and  an 
array  of  small  holes  with  the  diameter  of  1  mm  or  less  are  drilled  through  the  front  plate.  Several 
such  plates  were  made  with  different  number  of  holes  for  different  sample  sizes.  The  back  face  of 
the  sample  is  weakly  bonded  to  the  holder  with  a  very  thin  layer  of  vacuum  grease,  which  serves  to 
keep  vacuum  in  the  holder.  A  mechanical  rotary  pump  is  used  to  produce  the  vacuum  to  hold 
die  sample.  Therefore,  the  sample  holding  force  is  atmospheric  pressure.  This  arrangement  makes 
the  sample  under  test  almost  free  from  clamping  and  due  to  the  fact  that  the  whole  back  face  of  the 
sample  is  held  tighdy  against  the  holder,  the  possibility  of  the  sample  reflecting  face  tilting  is 
almost  eliminated.  The  experimental  results  on  several  PZT  family  samples  show  this  is  indeed  the 
case. 
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FIGURE  CAPTIONS 


Figure  1.  Schematic  drawing  of  the  ac  interferometer.  The  lock-in  can  be  operated  at  either  If  or 
2f  mode  to  measure  the  piezoelectric  or  electrostrictive  coefficients. 

Figure  2.  Interference  pattern  in  the  detection  plane  for  a  finite  beam  size  of  the  laser  and  a  large 
Ad.  I  is  the  light  intensity  and  L  is  the  distance  from  die  center  of  the  spot 

Figure  3.  A  model  system  to  analyze  the  motion  of  the  sample  in  the  z  direction. 

Figure  4.  (a)  Frequency  scan  of  dlt  for  X-cut  quartz.  F  is  the  operation  frequency. 

(b)  Frequency  scan  of  d33  for  hard  PZT. 

(c)  Voltage  scan  for  hard  PZT  at  900  HZ  driving  frequency.  dQ  is  the  amplitude  of  the 
displacement  and  V  is  die  amplitude  of  the  driving  voltage. 

Figure  5.  (a)  The  interference  light  intensity  versus  die  displacement  (from  eq.  (1».  This  is  the 
signal  detected  by  the  photo-detector. 

(b)  The  change  of  the  displacement  direction  causes  an  inflection  in  the  interference 
signal.  Here,  the  inflection  coincides  with  the  maximum  point  of  the  interference 
pattern.  The  detector  can  not  determine  the  turning  point. 

(c)  The  displacement  inflection  at  position  other  than  and  1^.  The  inflection  can 
be  seen  clearly. 

(d)  The  signals  from  the  main  detection  channel  and  the  quadrature  channel.  They  are 


Figure  6.  Schematic  drawing  of  the  dual-channel  detection  interferometer.  F  is  the  filter,  PS  is 
the  phase  plate,  W.G.  is  the  waveform  generator,  PBS  is  the  polarized  beam  splitter, 
ty4  is  the  quarter-wave  plate,  P  the  polarizer. 

Figure  7.  Recorded  signals  for  PLZT  8/65/35  samples,  (a)  is  the  signal  from  the  main  channel 
and  (b)  is  from  the  quadrature  channel  The  driving  triangle  voltage  is  also  recorded 
and  shown  as  (c).  All  the  signals  are  recorded  versus  time. 

Figure  8.  Strain  hysteresis  loop  for  PLZT  8/65/35. 


Figure  9.  Vacuum  sample  holder  for  the  strain  measurement.  The  holder  base  is  attached  to  an 
optical  stage  which  can  move  along  x,  y,  z  axis  and  rotate  above  x,  z  axis. 
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ISOTHERMAL  AGING  OF  THE  COMPLEX  DIELECTRIC  PERMITTIVITY  IN  MnO 
DOPED  LEAD  MAGNESIUM  NIOBATE:  LEAD  TITANATE  RELAXOR 
FERROELECTRIC  CERAMICS 


W.Y.  Pan,  Q.Y.  Jiang,  and  L.E.  Cross 
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ABSTRACT 

Isothermal  aging  of  the  complex  permittivity  was  investigated  in  MnO  doped  lead  magnesium 
niobate:  lead  titanate  (0.9  PMN:  0.1  PT)  relaxor  ferroelectric  ceramics.  Unlike  most  normal 
fenoelectrics,  both  real  and  imaginary  parts  of  the  permittivity  deviate  strongly  from  the  logarithmic 
aging  law,  both  curves  tending  to  saturate  at  longer  times.  The  effect  is  most  pronounced  in  the 
imaginary  part  of  the  permittivity  which  approaches  a  nearly  constant  value  quite  quickly  whilst  the 
real  part  continues  to  age.  The  results  can  be  explained  using  a  model  for  the  relaxor  in  which  there  is 
a  distribution  of  micro-polar  regions  of  varying  size  which  contribute  to  the  relaxation  component  of 
the  polarizability.  During  aging  it  is  postulated  that  defect  dipoles  induced  by  the  aliovalent  MnO 
doping  reorient  so  as  to  stabilize  the  polar  vector  in  the  micro-polar  regions.  Larger  regions  stabilize 
first,  removing  the  longer  relaxation  time  processes  and  leading  to  a  distribution  richer  in  shorter 
relaxation  times.  That  the  vs.  plot  taken  over  very  long  times  shows  no  evidence  of 
intersecting  the  K*  axis  at  a  large  value  suggests  that  the  micro-polar  regions  still  contribute  a  major 
part  of  the  polarizability  even  after  very  long  aging  times.  This  datum  is  consistent  with  evidence  from 
microwave  and  Raman  spectroscopy  that  conventional  soft  modes  are  not  major  contributors  to  the 
dielectric  response  of  relaxor  ferroelectric. 

Introduction 

In  normal  ferroelectric  ceramics,  at  temperatures  below  the  ferroelectric  Curie  point, 
homogeneity  of  composition  gives  rise  to  a  structure  of  polar  macrodomains.  For  many  such 
dielectrics  the  weak  field  dielectric  permittivity  is  not  strongly  dispersive  and  under  isothermal 
conditions  the  real  part  of  the  pe»  nittivity  decays  linearly  with  the  logarithm  of  time,  the  so-called 
logarithmic  aging  law. 

For  lead  zirconate:  titanate  ceramics,  Arlt  has  shown  that  there  is  a  strong  domain  wall 
contribution  to  the  permittivity,  and  that  it  is  this  dispersive  component  which  carries  the  aging 
phenomenon.  Plotting  real  (K^)  versus  imaginary  (K^)  components  as  a  function  of  time  exhibits  a 


series  of  stnight  lines  for  different  frequencies  which  extrapolate  for  •  time  to  a  constant  value  on  the 
id  axis.  It  is  suggested  that  this  value  is  the  average  of  the  intrinsic  single  domain  permittivities  (D, 
hi  re  lax  or  ferroelectrics,  the  situation  is  rather  different.  Compositional  heterogeneity  on  a  very 
fine  scale  gives  rise,  on  lowering  temperature,  to  a  sub-structure  of  polar  micro-regions  (2^  Earlier 
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studies  of  aging  in  PbMg^Nb^Oj:  PbTiOj  with  a  manganese  oxide  dopant  (3)  showed  that  for 

frequencies  in  the  relaxation  range  the  aging  is  not  a  linear  function  of  die  logarithm  of  time. 

This  paper  attempts  a  more  thorough  study  of  the  aging  process  in  this  relaxor  and  postulates  a 
model  to  account  for  the  behavior. 

Experimental  Procedure 

The  ceramic  samples  of  the  0.9  PMN:  0.1  PT  composition  used  in  this  study  were  prepared 
using  mixed  oxides  in  the  manner  described  earlier  by  Swartz  and  Shrout  (4,5).  In  this  processing 

MgO  and  Nt^O^  are  first  reacted  in  a  pre -calcine  to  form  the  columbite  structure  MgNt^O^  which  is 
then  reacted  with  PbO  to  for  the  PbjMgNl^Og  composition.  For  doping,  the  MnO  was  introduced 

using  a  dilute  solution  of  Mn  (NOj^  added  to  the  calcine.  The  final  sintered  samples  were  all  above 

97%  theoretical  density  and  were  shown  by  X-ray  diffraction  to  be  free  from  contamination  with  the 
pyrochlore  phase. 

Dielectric  constant  values  measured  on  freshly  de-aged  samples  showed  little  variation  with 
different  MnO  concentrations  in  the  range  up  to  0. 1%  (by  weight)  used  in  this  study. 

Aging  was  accomplished  in  an  isothermal  chamber  maintained  at  constant  temperature  by 
circulating  water  from  a  temperature  controlled  bath.  For  temperatures  above  22°  C  fluctuations  in  the 
aging  chamber  were  less  than  0. 1°  C  over  the  1 ,000  minute  aging  cycle.  Below  22°  C  control  was  by 
mixing  heated  and  cooled  water  and  was  better  than  0.2°  C.  Dielectric  properties  were  measured 
using  a  Hewlett  Packard  LCR  meter  (Model  427 A)  under  full  computer  control.  Before  measurement 
all  samples  were  de-aged  above  150°  C  and  field  levels  used  during  the  study  were  maintained  at  10 
v/cm  on  all  samples. 

Results  and  Discussion 

Figure  1  and  Figure  2  show  data  for  the  isothermal  aging  of  a  0.9  PMN:  0.1  PT  ceramic  with  a 
.1  (wt)  %  doping  of  MnO  measured  at  1, 10,  and  100  KHz  and  a  constant  37°  C  temperature.  Clearly 


series  of  straight  lines  for  different  frequencies  which  extrapolate  for  •  time  to  a  constant  value  on  the 
K1  axis.  It  is  suggested  that  this  value  is  the  average  of  the  intrinsic  single  domain  permittivities  (D, 

In  relaxor  ferroelectrics,  the  situation  is  rather  different.  Compositional  heterogeneity  on  a  very 
fine  scale  gives  rise,  on  lowering  temperature,  to  a  sub-structure  of  polar  micro-regions  (2^  Earlier 
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studies  of  aging  in  PbMgi/jNb^Cfy  PbTiOj  with  a  manganese  oxide  dopant  (3)  showed  that  for 

frequencies  in  the  relaxation  range  the  aging  is  not  a  linear  function  of  the  logarithm  of  time. 

This  paper  attempts  a  more  thorough  study  of  the  aging  process  in  this  relaxor  and  postulates  a 
model  to  account  for  the  behavior. 

Experimental  Procedure 

The  ceramic  samples  of  the  0.9  PMN:  0. 1  PT  composition  used  in  this  study  were  prepared 
using  mixed  oxides  in  the  manner  described  earlier  by  Swartz  and  Shrout  (4,5).  In  this  processing 

MgO  and  Nl^Oj  are  first  reacted  in  a  pre -calcine  to  form  the  columbite  structure  MgNt^O^  which  is 
then  reacted  with  PbO  to  for  the  Pt^Mghft^Og  composition.  For  doping,  the  MnO  was  introduced 

using  a  dilute  solution  of  Mn  (NC>3)2  added  to  the  calcine.  The  final  sintered  samples  were  all  above 

97%  theoretical  density  and  were  shown  by  X-ray  diffraction  to  be  free  from  contamination  with  the 
pyrochlore  phase. 

Dielectric  constant  values  measured  on  freshly  de-aged  samples  showed  little  variation  with 
different  MnO  concentrations  in  the  range  up  to  0.1  %  (by  weight)  used  in  this  study. 

Aging  was  accomplished  in  an  isothermal  chamber  maintained  at  constant  temperature  by 
circulating  water  from  a  temperature  controlled  bath.  For  temperatures  above  22°  C  fluctuations  in  the 
aging  chamber  were  less  than  0. 1°  C  over  the  1,000  minute  aging  cycle.  Below  22°  C  control  was  by 
mixing  heated  and  cooled  water  and  was  better  than  0.2°  C.  Dielectric  properties  were  measured 
using  a  Hewlett  Packard  LCR  meter  (Model  427A)  under  full  computer  control.  Before  measurement 
all  samples  were  de-aged  above  150°  C  and  field  levels  used  during  the  study  were  maintained  at  10 
v/cm  on  all  samples. 
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Results  and  Discussion 

Figure  1  and  Figure  2  show  data  for  the  isothermal  aging  of  a  0.9  PMN:  0.1  PT  ceramic  with  a 
.1  (wt)  %  doping  of  MnO  measured  at  1 , 10,  and  100  KHz  and  a  constant  37°  C  temperature.  Clearly 


neither  K*  nor  are  linear  functions  of  the  logarithm  of  time  and  this  is  typical  for  all  MnO  doped 
samples.  An  alternative  presentation  of  the  same  data  is  given  in  Figure  3,  where  K*  vs.  for  the 
three  frequencies  again  taken  at  37°  C  are  plotted.  At  the  lower  1  KHz  frequency  it  is  evident  that  the 
sample  loses  rapidly  with  only  a  small  change  in  and  that  for  all  frequencies  if  there  is  an 
intersection  of  the  axis  on  extrapolation  beyond  the  present  time  scale,  it  must  occur  at  much  lower 

K*  values.  Plots  of  the  relative  changes  AKVk^  vs  AK^/K^  (Figure  4)  again  confirm  the  non-linearity 
for  all  three  frequencies. 

That  the  aging  process  is  defect  controlled  is  confirmed  by  the  data  in  Figure  5  which  examines 
the  time  dependence  of  id  as  a  function  of  MnO  concentration.  Clearly,  the  undoped  samples  show 
no  aging  and  qualitatively  the  aging  rate  increases  with  increasing  MnO  concentration.  From 
isothermal  aging  studies  at  different  temperature  (Figure  6)  the  behavior  is  similar  to  that  for  room 
temperature  but  with  slightly  faster  aging  rates  at  the  lower  temperatures.  In  Figure  7  the  data  are 

replotted  so  as  to  emphasize  die  relative  changes  AK^/K^  vs  AK^/K^  bringing  the  curves  to  a 

common  origin  and  underscoring  the  decrease  of  die  change  in  loss  for  a  given  change  in  permitdvity 
as  the  temperature  increases. 

To  explain  the  observed  behavior,  we  speculate  that  the  aging  mechanism  is  the  same  as  that  in 
conventional  ferroelectric  dielectrics  and  is  due  to  the  slow  reorientation  of  defect  dipoles  under  the 
effective  field  of  the  polar  domains  stabilizing  the  existing  domain  configuration.  We  believe  that  the 
obvious  differences  in  the  behavior  of  the  telaxor  ferroelec tries  are  due  to  their  unusual  sub-structure 
of  polar  micro-regions  in  which  the  polar  vectors  are  being  constantly  reoriented  by  thermal 
fluctuations. 

The  proposed  explanation  is  then  as  follows:  in  the  temperature  range  above  the  dielectric 
maximum  Tc  the  PMN:  PT  is  made  up  of  small  polar  micro-regions  in  a  non-polar  but  very  high  K 

matrix  phase.  Because  the  regions  are  very  small  (~  200  A0  diameter)  their  Ps  vectors  are  unstable 

against  thermal  motion  and  are  constantly  switching  between  alternative  orientation  states.  It  is 
modulation  of  these  orientations  by  the  weak  external  measuring  field  which  contributes  a  major  part 
of  the  dielectric  permittivity. 
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The  polar  sub-domain  structure  is  engendered  by  the  local  compositional  heterogeneity  (local 
Mg:  Nb  ratio),  but  because  of  this  local  heterogeneity  translational  symmetry  is  lost  and  the  local 
(micro-region)  symmetry  is  lower  than  the  global  prototype  crystal  symmetry.  Under  this 

circumstance  the  various  orientations  for  Ps  in  the  micro-region  will  now  be  of  slightly  differing 

energies.  Now,  although  die  polar  vector  in  the  micro-region  will  still  suffer  step  wise  reorientation 
by  thermal  vibrations,  because  of  the  inequivalence  in  energy  states  it  will  tend  to  spend  longer  times 
in  the  lowest  energy  minimum,  thus  imparting  a  weak  time  averaged  polarization  to  the  micro-region. 

In  the  doped  ceramic,  this  time  average  polarization  will  act  to  reorient  the  Mn  induced  defect 
dipoles  so  as  to  slowly  stabilize  this  lowest  energy  state.  As  the  energy  of  this  state  is  lowered,  the 
ability  of  the  weak  measuring  field  to  modulate  the  orientation  is  lost  and  thus  the  dispersive 
component  of  the  polarizability  ages  out 

Exploring  now  the  effects  on  the  frequency  spectrum  of  permittivity,  on  the  model  of  flipping 
micro-polar  regions,  the  energy  barrier  over  which  Ps  must  be  carried  for  reorientation  (H)  must  be 

proportional  to  the  volume  of  the  micro-region.  Since  the  flipping  frequency  \>  is  given  by 

\)  *  i)0  e'H/^T 

where  \)0  is  a  lattice  vibration  frequency, 

K  is  the  Boltzmann  constant, 

H  is  the  local  energy  barrier. 

The  larger  micro-regions  will  be  slower  switching  but  will  also  have  the  larger  local  dipole 

moment  since  this  local  moment  n  is  simply  given  by  v  Ps  where  v  is  the  polar  volume. 

In  normal  dipolar  dielectrics,  the  energetics  are  such  that 

p£  «  kT  «  H. 

Under  these  circumstances,  the  fractional  ordering  imposed  by  the  E  field  is  trivial  and  the 
polarizability  a  linear  function  of  E  up  to  very  high  fields. 


For  the  ferroelectric  relaxor  the  local  moment  ji  is  much  larger  and  for  the  larger  micro-regions 


pE  may  approach  kT  at  modest  E  fields.  Evidence  for  this  energy  balance  may  be  adduced  from  the 

strong  non-linearity  of  the  polarization,  particularly  for  temperatures  close  to  Tc. 

In  the  doped  ceramic,  the  larger  slower  switching  regions  will  spend  longer  times  in  the 
preferred  orientation  position  and  thus  will  begin  to  reorient  the  dipolar  defect  structure.  The  orienting 

defects  will  generate  an  effective  dipolar  field  E<j  further  stabilizing  Ps,  as  the  preferred  well  deepens 

the  pE<j  term  will  slowly  override  kT  rendering  the  longer  relaxation  times  inoperative.  Thus  we 

expect  on  aging  that  the  distribution  of  relaxation  times  will  change  so  as  to  favor  the  shorter  times. 

This  scenario  is  depicted  schematically  in  Figure  8.  As  the  distribution  function  changes  due  to 
the  aging  process,  we  expect  the  trace  for  a  single  frequency  to  follow  a  path  such  as  that  delineated  in 
Figure  8,  Le.,  a  continuous  change  of  which  at  longer  times  shows  little  or  no  change  in  K^. 

At  present  the  explanation  proposed  above  is  largely  speculation,  however  the  pattern  is 

consistent  with  more  recent  studies  on  aging  under  DC  bias,  which  permit  an  assessment  of  Ej  the 

effective  defect  field  (6)  and  there  is  now  clear  evidence  of  polar  micro-regions  in  pure  PMN  which  are 
stabilized  at  nitrogen  temperature  and  observed  in  transmission  electron  microscopy  (7). 
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Figure  8.  The  proposed  trace  of  aging  in  the  Cole-Cole  plots  for  the  system  containing  a  distribution 
of  domain  size 
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THERMODYNAMIC  THEORY  OF  THE  LEAD  ZIRCONATE-TTTANATE 
SOLID  SOLUTION  SYSTEM,  PART  I:  PHENOMENOLOGY 

M.  J.  Haun^,  E.  Furman^,  S.  J.  Jang,  and  L.  E.  Cross 

Materials  Research  Laboratory 
The  Pennsylvania  State  University 
University  Park,  PA  16802 

Compositions  within  the  lead  zirconate-titanate  (PZT)  solid  solution  system  have 
been  extensively  used  in  polycrystalline  ceramic  form  in  a  wide  range  of  piezoelectric 
transducer  applications.  However,  the  growth  of  good  quality  PZT  single  crystals  for 
compositions  across  the  entire  phase  diagram  has  not  been  accomplished.  Due  to  the  lack  of 
single-crystal  data,  an  understanding  of  the  properties  of  polycrystalline  PZT  has  been 
limited.  If  the  single  domain  (intrinsic)  properties  of  PZT  could  be  determined,  then  the 
extrinsic  contributions  (e.g.  domain  wall  and  defect  motions)  to  the  polycrystalline  properties 
could  be  separated  from  the  intrinsic  contributions.  The  purpose  of  this  research  has  been  to 
develop  a  thermodynamic  phenomenological  theory  to  model  the  phase  transitions  and 
single-domain  properties  of  the  PZT  system. 

This  paper  is  the  first  of  a  series  of  five  papers  describing  the  thermodynamic  theory 
of  PZT  that  has  been  developed  for  the  entire  solid  solution  system.  In  this  paper  the 
previous  work  that  led  to  the  present  theory  will  first  be  reviewed,  followed  by  the 
presentation  of  the  energy  function  for  PZT,  and  the  solutions  and  property  relations  that  can 
be  derived  from  this  energy  function. 
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I.  INTRODUCTION 


The  lead  zirconate-titanate  (PZT)  phase  diagram 1  is  shown  in  Figure  1 .  A  cubic 
paraelectric  phase  (P^)  occurs  at  high  temperatures  and  has  the  perovskite  crystal  structure 
ABO3.  On  the  lead  titanate  (PbTi03)  side  of  the  phase  diagram,  a  ferroelectric  tetragonal 
phase  (Fy)  exists  with  a  spontaneous  polarization  along  the  pseudocubic  [001  ]  direction.  A 
morphotropic  boundary  separates  the  tetragonal  phase  from  a  ferroelectric  high-temperature 
rhombohedral  phase  (Fr(HT))- 

Another  ferroelectric  to  ferroelectric  phase  transition  occurs  between  the 
high-temperature  rhombohedral  phase  and  a  low-temperature  rhombohedral  phase  (Fr^lT))- 
Both  of  these  rhombohedral  phases  have  a  spontaneous  polarization  that  occurs  along  the 
[111]  direction.  The  low -temperature  rhombohedral  phase  has  a  tilting  or  rotation  of  the 
oxygen  octahedra  about  the  [1 1 1]  axis,  which  does  not  occur  in  the  high-temperature  phase. 

On  the  lead  zirconate  side  of  the  phase  diagram  antiferroelectric  tetragonal  (Ay)  and 
orthorhombic  (Aq)  phases  are  present.  These  antiferroelectric  phases  are  composed  of  two 
sublattices  with  equal  and  opposite  polarization,  resulting  in  zero  net  polarization. 

More  recent  data^  have  shown  that  the  Ay  phase  does  not  occur  in  pure  lead 
zirconate,  but  will  occur  with  small  dopants  or  impurities  such  as  strontium.  These  data  have 
also  shown  that  the  high-temperature  rhombohedral  phase  extends  over  to  lead  zirconate. 

The  theory  described  in  this  paper  will  be  developed  to  model  the  phase  diagram  according  to 
these  recent  data. 

PZT  compositions  have  important  technological  applications  as  piezoelectric 
transducers,  pyroelectric  detectors,  electro-optic  devices,  and  explosively  induced  charge 
storage  devices.  The  main  application  of  PZT  has  been  as  piezoelectric  transducers  with 
compositions  near  the  tetragonal-rhombohedral  morphotropic  boundary.  These  applications 
include  phonograph  pickups,  buzzers,  relays,  accelerometers,  igniters,  micropositioners, 
hydrophones,  sonar,  wave  Filters,  earphones,  delay  line  welders,  cutters,  and  high  voltage 
sources.  I  Recently  piezoelectric  motors  have  even  been  made.^ 


Lead  titanate  has  become  important  for  possible  application  in  hydrophones  and 
medical  ultrasonic  imaging,  because  of  the  large  electromechanical  anisotropy  that  can  be 
obtained  in  lead  titanate  polycrystalline  ceramics^.  The  theory  developed  in  this  paper  has 
been  used  to  explain  this  large  anisotropy,  which  occurs  in  lead  titanate  in  polycrystalline 
form,  but  not  in  single-crystal  fomA 

PZT  compositions  with  the  high  to  low  temperature  rhombohedral  phase  transition 
have  been  used  as  pyroelectric  detectors^.  When  cooling  through  this  transition,  the 
polarization  increases  suddenly  due  to  a  contribution  from  the  tilting  of  the  oxygen  octahedra 
in  the  low-temperature  phase.  This  results  in  a  large  pyroelectric  coefficient.  The  dielectric 
constant  is  relatively  low,  and  only  has  a  very  slight  change  when  going  through  this 
transition.  A  large  pyroelectric  coefficient  p  (=  dPg/dT)  and  small  dielectric  constant  K 
result  in  a  large  pyroelectric  figure  of  merit  (=  p  /K^). 

PZT  is  also  an  important  material  for  electrooptic  applications,  when  doped  with 
lanthanum  to  form  PLZT.  According  to  Haertiing  and  Land**,  PLZT  has  advantages  in 
electrooptic  applications  because  of  ”1)  high  transparency  essentially  independent  of  grain 
size,  2)  ease  of  fabrication,  and  3)  electrically  variable  optical  retardation,  including 
electrooptic  memory  and  linear  or  quadratic  modulation  capabilities". 

Explosively  induced  charge  storage  devices  have  also  been  fabricated  from  PZT 
compositions  utilizing  the  ferroelectric  to  antiferroelectric  transition.  In  these  applications  a 
stress  is  applied  to  a  poled  ceramic  sample  with  the  ferroelectric  structure,  causing  it  to 
transform  to  the  antiferroelectric  state  and  to  suddenly  release  the  stored  charge. 

The  applications  described  above  show  that  PZT  is  an  important  material.  These 
applications  use  PZT  in  polycrystalline  ceramic  form.  The  growth  of  good  quality  PZT 
single  crystals  for  compositions  across  the  entire  phase  diagram  has  not  been  accomplished.^ 
Due  to  the  lack  of  single-crystal  data,  an  understanding  of  the  properties  of  ceramic  PZT  has 
been  limited.  If  the  single-domain  properties  of  PZT  could  be  determined,  then  these 
intrinsic  contributions  to  the  ceramic  properties  could  be  separated  from  the  extrinsic 


contributions  (e.g.  domain  wail  and  defect  motions).  The  purpose  of  this  research  has  been 
to  develop  a  thermodynamic  phenomenological  theory  to  predict  the  phase  transitions  and 
single-domain  properties  of  the  PZT  system.  This  theory  could  then  be  used  to  further  the 
understanding  of  the  properties  of  polycrystalline  materials. 

This  paper  is  the  first  of  a  series  of  five  papers  covering  the  details  of  the 
thermodynamic  theory  of  PZT  that  has  been  developed  for  the  entire  solid  solution  system. 

In  this  paper  the  previous  work  that  led  to  the  present  theory  will  first  be  reviewed,  followed 
by  the  presentation  of  the  energy  function  for  PZT,  and  the  solutions  and  property  relations 
that  can  be  derived  from  this  energy  function.  In  the  next  three  papers  the  coefficients 
of  the  energy  function  will  be  determined  from  experimental  data.  Finally  in  the  last  paper  of 
the  series^  theoretical  calculations  of  the  PZT  single-domain  properties  will  be  presented, 
along  with  a  discussion  of  the  applications  of  the  theory. 

D.  DEVELOPMENT  OF  A  THERMODYNAMIC  THEORY  OF  PZT 

The  Landau-G  insburg  -  Devonshire  type  phenomenological  theory  was  used  to 
explain  the  transitions  and  properties  at  the  tetragonal-rhombohedral  morphotropic  phase 
boundary  in  the  PZT  system.  ^21  However,  due  to  the  lack  of  single-crystal  data  the 
coefficients  of  the  energy  function  could  not  be  adequately  determined. 

Over  the  past  several  years  at  the  Materials  Research  Laboratory  of  the  Pennsylvania 
State  University,  the  coefficients  of  an  energy  function  for  the  PZT  system  have  been 
determined  using  indirect  methods.  The  development  of  this  theory  occurred  in  a  series  of 
steps,  where  different  parts  of  the  phase  diagram  were  modeled  separately.  This 
development  will  be  described  below. 

Amin ^  developed  a  modified  Devonshire  form  of  the  elastic  Gibbs  free  energy 
function  for  the  simple  proper  ferroelectric  phases  of  the  PZT  system.  The  spontaneous 
strain  in  the  ferroelectric  phases  was  determined  using  high-temperature  x-ray  diffraction. 
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and  used  to  calculate  the  spontaneous  polarization  through  the  electrostrictive  constants. 
These  data  were  used  to  determine  the  higher  order  dielectric  stiffness  coefficients  (aj  j  and 
am).  The  morphotropic  boundary  between  the  tetragonal  and  rhombohedral  phases  was 
then  used  to  determine  the  remaining  higher-order  dielectric  stiffness  coefficients  (a^. 
dj  i2»  ®»d  0^23).  and  allowed  for  the  first  time  the  calculation  of  the  dielectric,  piezoelectric, 
elastic,  and  thermal  properties  of  the  ferroelectric  single-domain  states.  This  theory  was 
developed  for  the  region  of  the  PZT  system  from  lead  titanate  over  to  the  morphotropic 
boundary  between  the  tetragonal  and  rhombohedral  phases,  and  could  be  extrapolated  into 
the  high-temperature  rhombohedral  phase  field. 

In  this  initial  theory  the  Curie  constant  was  assumed  to  be  independent  of 
composition.  Amin  et  al.^  later  found  from  a  combination  of  calorimetric  and 
phenomenological  data  that  the  Curie  constant  was  dependent  on  composition  with  a  peak 
forming  near  the  morphotropic  boundary.  The  theory  was  then  modified  to  account  for  the 
compositional  dependence  of  the  Curie  constant.24  This  resulted  in  better  agreement 
between  the  theoretical  and  experimental  dielectric  data  near  the  morphotropic  boundary. 

To  account  for  the  tilting  of  the  oxygen  octahedra  in  the  low-temperature 
rhombohedral  phase  Halemanc  et  al.2®  expanded  the  energy  function  in  a  one-dimensional 
power  series  of  the  polarization  (P)  and  tilt  angle  (8)  assuming  isothermal  and  zero  stress 
conditions.  The  possible  phase  transition  sequences  and  solutions  to  the  energy  function 
were  then  investigated.  This  theory  was  applied  to  the  Pb(ZrQ  gTiQ  ^03  composition,  first 
by  assuming  2nd  order  phase  transitions,2®  and  later  by  assuming  1st  order 
transitions2® ,2^. 

The  energy  function  was  then  expanded  in  a  three-dimensional  power  series  of 
P  and  0  and  used  to  describe  all  of  the  ferroelectric  phases  of  the  PZT  system  assuming  1st 
order  phase  transitions.2^  The  coefficients  of  the  energy  function  were  determined  from  the 
phase  boundary  and  equilibrium  conditions  combined  with  experimental  data,  resulting  in 
smooth  continuous  functions  across  the  phase  diagram.  The  experimental  and  theoretical 
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phase  diagrams  were  shown  to  agree  very  well.  Using  this  theory  the  spontaneous 
polarization  and  tilt  angle  were  calculated  as  a  function  of  composition  and  temperature. 

A  tricritical  point,  where  a  phase  transition  changes  from  first  to  second  order,  was 
found  to  occur  between  the  cubic  and  rhombohedral  phases  at  the  PZT  94/6  composition 
(94%  PZ  and  6%  PT).29-3 1  From  lead  zircon  ate  to  the  tricritical  point  the 
cubic-rhombohedral  transition  was  shown  to  be  first  order,  and  then  to  change  to  second 
order  from  the  tricritical  point  over  to  at  least  the  PZT  88/12  composition.^  Lead  titanate 
has  a  first  order  transition  from  cubic  to  tetragonal,  and  thus  a  second  tricritical  point  should 
occur  between  the  PZT  88/12  and  lead  titanate  compositions,  where  the  transition  would 
change  back  to  first  order. 

To  provide  additional  data  to  determine  the  coefficients  of  the  energy  function  and  to 
locate  the  second  tricritical  point,  pure  homogeneous  sol-gel  derived  PZT  powders  were 
prepared  for  several  compositions  in  the  rhombohedral  phase  field. The  lattice  parameters 
of  these  compositions  were  determined  from  high-temperature  x-ray  diffraction,  and  used  to 
calculate  the  spontaneous  strain.^  By  using  these  data  to  determine  the  higher-order 
dielectric  stiffness  coefficients  of  the  energy  function,  the  second  order  transition  region  was 
found  to  extend  over  to  near  or  possibly  at  the  morphotropic  boundary.^4  Additional  details 
of  the  tricritical  behavior  in  PZT  will  be  described  in  the  second  paper  in  this  series.  ^ 

Using  the  x-ray  data  described  above  a  phenomenological  theory  was  developed  to 
account  for  the  2nd  order  transition  region  of  the  PZT  system. ^  In  this  theory  additional 
terms  were  added  to  the  energy  function  to  account  for  the  rotostrictive  coupling  between  the 
stress  and  the  square  of  the  tilt  angle  (analogous  to  the  electrostrictive  coupling).  A 
rotostrictive  constant  was  determined  from  experimental  spontaneous  strain,  polarization, 
and  tilt  angle  data  using  a  relation  derived  from  the  energy  function  for  the  low-temperature 
rhombohedral  phase.  The  rotostrictive  contribution  to  the  spontaneous  strain  was  found  to 
be  opposite  in  sign  to  the  electrostrictive  contribution. 
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A  more  complete  phenomenological  theory  of  lead  titanate  was  recently  developed 
independently  of  the  PZT  theory  using  the  available  single  crystal  data.  35  The  theoretically 
predicted  spontaneous  polarization  and  strains,  and  dielectric  and  piezoelectric  properties  are 
in  good  agreement  with  the  experimental  data.  The  results  of  this  theory  were  used  to  show 
that  the  large  electromechanical  anisotropy  found  in  lead  titanate  ceramics,  but  not  present  in 
the  single  crystal  properties,  is  simply  due  to  the  intrinsic  averaging  of  the  electrostrictive 
constants. 

A  phenomenological  theory  for  lead  zirconate  was  also  developed  independently  of 
the  PZT  theory  using  the  available  single-crystal  data.  37  This  theory  was  used  to  calculate 
the  phase  stability,  antiferroelectric  spontaneous  strain,  and  ferroelectric  polarization  and 
strain.  The  polycrystalline  ceramic  dielectric  constant  was  calculated  from  the  single-crystal 
constants  in  this  theory,  and  found  to  be  in  good  agreement  with  microwave  dielectric  data. 
This  lead  zirconate  theory  was  also  extended  into  the  the  PZT  system  to  account  for  the 
antiferroelectric  region  of  the  phase  diagram. 

As  described  above  phenomenological  theory  has  been  developed  for  different 
sections  of  the  PZT  phase  diagram.  The  purpose  of  this  senes  of  papers  is  to  present  the 
recent  work  that  has  been  completed  to  combine  these  separate  theories  into  a  complete 
theory  of  PZT,  where  a  single  energy  function  and  set  of  coefficients  can  be  used  to  calculate 
the  phase  stability  and  properties  of  the  entire  PZT  system.  In  the  next  section  this  energy 
function  will  be  presented. 


ED.  FREE  ENERGY  FUNCTION  FOR  PZT 
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In  developing  a  phenomenological  theory  for  the  PZT  system  the  order  parameters 
that  cause  the  phase  transitions  must  be  accounted  for.  The  tetragonal  and  high-temperature 
rhombohedral  phases  undergo  proper  ferroelectric  transitions  from  the  paraelectric  cubic 
state,  where  the  spontaneous  polarization  is  the  order  parameter  causing  the  phase  transition. 
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An  improper  ferroelectric  phase  transition  occurs  between  the  high  and  low  temperature 
rhombohedral  phases ,  where  the  spontaneous  tilting  of  the  oxygen  octahedra  causes  the 
phase  transition  and  contributes  to  the  spontaneous  polarization.  An  antiferroelectric  type 
polarization  develops  when  transforming  from  the  ferroelectric  to  antiferroelectric  phases  and 
is  therefore  the  order  parameter  for  the  antiferroelectric  state. 

Thus  to  account  for  all  of  the  phase  transitions  in  the  PZT  system  the  energy 
function  should  include  the  ferroelectric  and  antiferroelectric  measures  of  the  polarization, 
and  the  tilting  of  the  oxygen  octahedra  as  order  parameters.  A  two-sublattice  modeP^  was 
used  to  derive  the  ferroelectric  and  antiferroelectric  measures  of  the  polarization.  To  account 
for  the  tilting  of  the  oxygen  octahedra  an  additional  term  was  added  to  the  elastic  Gibbs  free 
energy  as  described  in  References  25-28. 

The  following  energy  function  for  the  PZT  system  was  then  derived  from  the 
symmetry  of  the  paraelectric  phase  (m3m)  assuming  isothermal  conditions,  using  reduced 
notation,  and  expanding  the  energy  function  in  powers  of  the  ferroelectric  (P,)  and 
antiferroelectric  (pj)  polarizations,  oxygen  octahedral  tilt  angle  (0j),  and  including  couplings 
between  these  order  parameters  and  between  the  stress  (Xm)  and  the  order  parameters: 

A  G  =  04  [P^  +  P22  +  P32]  +  04  l  [Pi4  +  P24  +  P34] 

+  <*12  Pi2  p22  +  p22  P32  +  P32  Pi2)  +  <*l 1  l[Pl6  +  p26  +  p36l 
+  an2  [Pi4  (P22  +  P32)  +  P24  (Pi2  +  P32)  +  P34  (Pi2  +  P22)J 
+  al23  pl2  p22  p32  +  al  CPl2  +  P22  +  P32l  +  al  l  tPl4  +  P24  +  P34l 
+  ai2[Pi2P22  +  P22P32  +  P32Pl2]  +01U  [Pl6  +  P26  +  P36l 
+  12  [Pl4  (P22  +  P32)  +  P24  (Pl2  +  P32)  +  P34  (Pl2  +  P22>J 

+  °123  Pl2  P22  P32  +  P-l  1  [Pl2  Pl2  +  P22  P22  +  p32  P32l 
+  M-12  [pl2  (P22  +  P32)  +  p22  (Pl2  +  P32)  +  P32  (Pl2  +  P22)l 
+  M-44  [P 1  P2P  1 P2  +  P2P3P2P3  +  P3P1P3P1I  +  Pl  (®12  +  ®22  +  0321 

+  Pj  1  [9i4  +  e24  +  e34]  -*-Yn  [p12012  +  p22  e22  +  p32  e32i 
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+  Yl2  [Pi2  (0(2  +  022)  +  P22  (®I2  +  032)  +  P32  Oi2  +  022^  j 

+  Y44  [Plp2el92  +  p2p3®2e3  +  P3Ple30lJ 

-  1/2  Su[X!2+  X22+  X32J  -  S12  [XjX 2+  X2X3+  XjXJ 

-  1/2  S44  [X42+  X62+  X^2]  -  Qn  [X ! P ! 2  +  X2P22  +  X3P32] 

-  Ql2  [Xj  (P22  +  P32)  +  X2  (Pt2  +  P32)  +  X3  (Pt2  +  P22)] 

-  Q44  [X4  P2p3  +  X5  pip3  +  X*  PlP2]  -  Zj !  [XlPl2  +  X2p22  +  X3p32] 

-  Z[2  [Xi  (p22  +  p32)  +  X2  (pi2  +  p32)  +  X3  (Pl2  +  p22)] 

-  Z44  [X4P2P3  +  X5  Plp3  +  X^  p1p2]  -R[|  [Xj0j2  +  X2022  +  X3032] 

-  Rj2  [Xj  (022  +  032)  +  X2  (0j2  +  032)  +  X3  (0j2  +  022)] 

-  [X4  0203  +  X3  0j03  +  X^  0j02]  (1) 

The  coefficients  of  this  energy  function  are  defined  in  Table  I.  The  energy  function  includes 
all  possible  ferroelectric  and  antiferroelecthc  polarization  terms  up  to  the  sixth  order,  tilt  angle 
terms  up  to  the  first  fourth  order  term,  and  only  the  first  order  coupling  terms. 

IV.  SOLUTIONS  TO  THE  ENERGY  FUNCTION 

The  following  solutions  to  the  energy  function  (Equation  1 )  are  of  interest  in  the 
PZT  system: 

Paraelectric  Cubic  (Pq) 

pl  =  p2  =  p3  =  0.  Pl  =  P2  =  P3=0’  0j=02  =  03  =  O  (2) 

Ferroelectric  Tetragonal  <Fj) 

pi  =  p2  =  0,  p32*o,  pi  =  p2  =  p3  =  0,  el=e2  =  e3  =  o  (3) 

Ferroelectric  Orthorhombic  (Fq) 

pl  =0,  P22  =  P32*0,  p1=p2  =  p3=0,  0j=02  =  03  =  O  (4) 


WWW 


Ferroelectric  High-temperature  Rhombohedral  (Fpffff)) 

P12  =  P22  =  P32*0.  Pi  =P2  =  P3  =  °>  0j  =  02  =  63  =  0  (5) 

Ferroelectric  Low-temperature  Rhombohedral  (F^j_j)) 

Pj2  =  P22  =  P32^0.  Pi  =  P2  =  P3  =  ©l2  *  ®22  *  932  *  0  (6) 

Antiferroelectric  Orthorhombic  (Aq) 

Pl  =  P2  =  P3  =  0,  pt  =0,  p22  =  p32*0,  01  =02  =  03  =  O  (7) 

All  of  these  solutions,  except  for  the  ferroelectric  orthorhombic  solution,  are  stable  in  the 
PZT  system.  The  ferroelectric  orthorhombic  solution  was  also  included  here,  because  the 
coefficients  necessary  to  calculate  the  energy  of  this  phase  can  be  determined.  An 
independent  check  of  the  coefficients  can  then  be  made  by  checking  if  this  phase  is 
metastable  across  the  PZT  system,  as  it  should  be  experimentally. 

Applying  these  solutions  to  Equation  ( 1 )  under  zero  stress  conditions  results  in  the 
following  relations  for  the  energies  of  each  solution: 

Pc  AG  =  0  (8) 

F'p  AG  =  <*1^3^  +  1^3^ +  1 1^*3^  (9) 

Fq  AG  =  2 04P32  +  (2  otj2)  P34  +  2  (aju+  aj  j2)  P3**  (10) 

Fr(HT)  AG  =  3  a^P32  +  3  (a^+ a^2)  P3^  +  (3  aj  j  j+ 6  aj  i2  +  a^23>  P3^  (11) 

fR(LT)  ag  =  3  alp32  +  3  (al  1+  al 2)  p34  +  (3  all  l+  6  04  12  +  a123)  P36 

+  3(3l032  +  3  p11034  +  3  (Yii+2Y12  +  Y44)P32032  (12) 

Aq  AG  =  2  0[P32  +  (2  <Jj2)  P34  +  2  (a^  j  Oj  12)  P3^  (13) 

The  spontaneous  ferroelectric  and  antiferroelectric  polarizations  (P3  and  P3)  and  tilt 
angle  (©3)  in  the  above  equations  can  be  found  from  the  first  partial  derivative  stability 


conditions  (dAG/dPj,  3AG/dp3<  and  9AG/d03)  as  shown  below: 


Ff  9AG/3P3  =  0  =  3  j?34  +  2  aj  ^3^  -t-  <X[  (14) 

Fq  9AG/3P3  =  0  =  3  (a^ u+ cxj^) p34  +  (2 ct[ j+ a^) P32  +  ct[  (15) 

fR(HT)  aAG/aP3  =  0  =  (3  a! ! !+  6  aj  12  +  a123>  P34 

+  2(a11+a12)P32  +  al  (16) 

FR(LT)  3AG/3P3  =  0  =  (3  aj  6  aj  j2  +  ai23^  P34 

+  2(all+a12)P32  +  a1+Yn032  (17) 

dAG/003  =  0  =  (3 1  +  2(3j[032-*-  Y11P32  (18) 

Aq  aAG/dp3  =  0  =  3  (ajj  1+ 0|  12)  P34  +  (2  <Jn+ 0[2)  P3^  +  ^1  (19) 


The  polarizations  and  tilt  angle  can  be  calculated  by  solving  these  quadratic  equations. 
Equations  (9)  -  (19)  relate  the  energies  of  each  solution  to  the  coefficients  of  the  energy 
function.  Thus  by  determining  these  coefficients,  the  energies  of  each  phase  can  be 
calculated. 

V.  SPONTANEOUS  ELASTIC  STRAINS 


The  spontaneous  elastic  strains  x,  (=  9AG/5X|)  under  zero  stress  conditions  can  be 
derived  from  Equation  ( 1 )  as  follows: 


pc 

*1  ~  *2  =  x3  ~  x4  —  *5  =  *6  —  ^ 

(20) 

ft 

xl  =x2  =  Ql2P32-  *3  =  Qi1P32’  x4  =  x5  =  x6  =  ° 

(21) 

Fo 

xl  =  2Ql2P32’  x2  =  x3  =  ^Qll+Ql2>p32 

x4  =  Q44  P32’  x5  =  x6  =  ® 

(22) 

FR(HT)  X1  =*2  =  *3=<Qll+2Ql2>p32>  *4  =  *5  =  x6  =  Q44  p32 


(23) 


(24) 


FR(LT)  X1  =  x2  =  x3  =  (Ql  1+  2  Ql2>  p32  +  <R1 1+  2  R12}  032’ 
x4  =  x5  =  x6  =  Q44P32  +  R44  032 

A0  x1  =  2Zl2P32.  x2  =  x3=(Z11+Z12)P32 

x4  =  Z44  P32,  X5  =  xg  =  0  (25) 

In  the  next  three  papers  of  this  series  2  these  spontaneous  strain  relations  will  be 
shown  to  be  very  important  in  determining  the  coefficients  of  the  energy  function. 
Spontaneous  strain  data  will  be  determined  from  x-ray  diffraction  of  PZT  powders,  and  used 
with  the  electrostrictive  constants  to  calculate  the  spontaneous  polarization,  which  is  needed 
to  determine  coefficients  of  the  energy  function. 


VI.  DIELECTRIC  PROPERTIES 

Relations  for  the  relative  dielectric  stiffnesses  Xij  (=  32AG/3P-dPj)  were  derived 
from  Equation  ( 1 )  for  the  six  solutions: 

PC  Xn=X22  =  X33  =  2eOab  Xl2  =  X23=X31=0  U6> 

FT  Xu  =  X22  =  2  e0  lal  +  al  lp32  +  al  1 1P34J- 

X33  =  2e0lal +6allp32+ l5alllp34^  Xl2  =  X23  =  X31  =°  (27> 

F0  X 1 1  =  2  £q  [«i  +  2  a12  P32  +  (2  a112  +  a123)  P34], 

X22  =  X33  =  2€0taI  +  (6«ll  +aI2)P32  +  (15ani  +7au2>P34J. 

Xl2  =  X31=°-  X23=4e0ta12p32  +  4all2p34^  (28) 

FR(HT)  Xl  1  =  X22  =  X33  =  2  e0  £al  +  al  1+  2  a12>  P32 
+  (15am  +  14aI12  +  al23  )  P34]. 

Xl2  =  X23  =  X31  =  4eo£a12P32  +  £4all2  +  a123  )P341  (29) 

FR(LT)  Xu  =X22  =  X33  =  2eo£al +£6all+2a12)p32 

+  (15am  +  14a112  +  al23  )P34  +  (YU+2Y12)e32J, 

Xl2  =  X23  =X31  =  4e0la12P32  +  £4a112  +  a123  >  p34  ■*"  Y44  0321  (30) 


xnyji  kt.  Xy^  vj~  vw  k i  *o  n*  \j* 
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A0  Xu=2e0[a1  +  2^12P32],  X22  =  X33  =  2e0tal +^11 +^12)P32!- 

Xl2  =  X31 =0-  X23=£0^44P32  <3l> 

The  multiplication  by  the  permittivity  of  free  space  £q  in  these  equations  was  required  to 
convert  from  absolute  to  relative  dielectric  stiffnesses.  Equations  (26)  -(31)  can  be  used  to 
calculate  the  relative  dielectric  stiffnesses  for  each  phase  based  on  the  original  cubic  axes. 

In  the  orthorhombic  state  the  polarization  can  be  along  any  of  the  <1 10>  directions 
of  the  original  cubic  axes.  The  polarization  of  the  rhombohedral  state  can  be  along  any  of  the 
<  1 1 1>  directions.  By  rotating  these  axes  so  that  for  both  states  the  new  xj  axis  is  along  the 
polar  directions,  diagonalized  matrices  will  result.  The  new  dielectric  stiffness  coefficients 
(indicated  by  a  prime)  can  be  related  to  the  old  coefficients  [defined  by  equations  (28)  -(31)1 
with  the  following  relations: 

F0andA0  Xu'=Xll.  X22  =  X33  '  X23 

X33  =  X33  +  X23-  X12  =  X23  =  X3l'  =  0  <32> 

FR(HT)andFR(LT)  Xir  =  X22  =  Xll -Xl2’  X33'  =  Xll  +  2  Xl2 

Xl2  =  X23  =  X31  =  0  (33> 

These  equations  can  be  used  to  calculate  the  dielectric  stiffnesses  of  the  orthorhombic  and 
rhombohedral  phases  parallel  and  perpendicular  to  the  polar  axes. 

The  dielectric  susceptibility  coefficients  (rjy)  can  be  determined  from  the  reciprocal 
of  the  dielectric  stiffness  matrices  (X[j)  using  the  following  relation3^: 

r|ij  =  Ajj/A,  (34) 

where  Ajj  and  A  are  the  cofactor  and  determinant  of  the  Xy  matrix.  Using  this  relation  results 
in  the  following  relations  for  the  dielectric  susceptibility  coefficients  (r|ij); 

pc  Thi=sn22  =  n33  =  l&ii-  ni2  =  n23  =  Ti3i  =  °  (35> 

fT  ^11  =  n22  =  */Xl  1-  H33=l/X33  (36) 
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I 

•0 

I 


Bunrow. 


F0andA0  Tlii  =  l/Xll* 

ni2*n3i =0* 


(37) 


n22  =  ,l33  =  X3y(X332  -X232)* 
n23  =  *X23/(X332-X232) 


nil’*  1/Xll*.  H22  =  l& 22-  ^33  =  l/X33' 

ni2,  =  n23'  =  n3i'  =  °  (38> 

FR(HT)andFR(LT) 

nil  =  H22  =  H33  =  ^Xll2  -  X122) / (Xl l3  -  3  Xl  1X122  +  2  Xl23) 
ni2  =  n23  =  n3l  =<Xl22-XllXl2)/(Xll3-3XllXl22  +  2Xl23)  09) 

nn,=n22,=  1/xii'.  n33  =  I/X33' 

ni2,i=n23,  =  n3i'  =  0 

These  equations  can  be  used  to  calculate  the  dielectric  susceptibilities  of  each  phase  from  the 
coefficients  of  the  energy  function. 


VH.  PIEZOELECTRIC  PROPERTIES 

Relations  for  the  piezoelectric  by  coefficients  (=  d2AG/9Pj5Xj)  were  derived  from 
Equation  ( 1)  for  the  tetragonal  and  rhombohedral  states  as  shown  below: 

Fj  b33  =  2QllP3’  b31  =  b32  =  2Ql2P3’ 

b15  =  b24  =  044  P3’  bll  =bi2  =  bi3  =  b14  =  b16  =  0' 
b21  =b22  =  b23  =  b25  =  b26  *  b34  =  b35  *  b36  =  0 

FR(HT) and  FR(LT) 

bll  =  b22  =  b33  =  2  Ql  lp3>  b14  =  b25  =  b36  =  0 
b12  =  bi3  =  b2i  =  t>23  =  b3 1  =  b32  =  2  Q 1 2  P3- 
b15  *  b16  =  b24  =  b26  =  b34  =  b35  =  $44  p3’ 

Since  a  coupling  term  of  the  type  XjPjQj  was  not  included  in  Equation  ( 1 ).  the  by 
relations  [Equation  (42)]  for  the  high  and  low  temperature  rhombohedral  phases  are  of  the 
same  form.  However,  the  spontaneous  polarizations  P3  are  defined  by  different  relations  for 


(41) 


(42) 
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the  high  and  low  temperature  rhombohedrai  phases,  and  thus  different  values  would  result 
for  these  coefficients. 

The  piezoelectric  dy  coefficients  are  defined  by4®: 

dij'bkjHik  <43> 

Using  this  relation  for  the  tetragonal  and  rhombohedrai  states  results  in  the  following 
relations: 


Ft  d33  =  2eO®33Qllp3-  d31  =  d32  =  2e0Tl33Ql2p3’ 

d15  =  d24=  Co  ^1 1  Q44P3*  dn  =dl2  =  d13  =d14  =  d16  =  0) 

d21  =  d22  =  d23  ~  d25  =  d26  =  d34  =  d35  *  d36  =  0  ^44) 

FR(HT)andFR(LT) 

dll  =d22  =  d33  =  2e0^llQll  +2Tll2Ql2>  p3> 

d12  =  dl3  =  d21  =d23  =  d31  =d32  =  2e0^11  Ql2  +  Tll2^Qll  +  Ql2>I  p3« 

d14  “  d25  =  d36s2e0  ^12  044  p3  • 

d15  =  dl6  =  d24  =  d26  =  d34  =  d35=:e0(T>ll  +  r'l2)044p3-  (45) 

The  multiplication  by  the  permittivity  of  free  space  £q  in  these  equations  was  required  to 
convert  the  dielectric  susceptibilities  from  relative  to  absolute.  Equations  (41),  (42),  (44), 
and  (45)  can  be  used  to  calculate  the  piezoelectric  by  and  dy  coefficients  of  the  tetragonal  and 
rhombohedrai  phases  from  the  coefficients  of  the  energy  function. 

Vin.  SUMMARY 

The  applications  of  compositions  of  the  PZT  solid  solution  system  as  piezoelectric 
transducers,  pyroelectric  detectors,  electro-optic  devices,  and  explosively  induced  charge 
storage  devices  were  described  in  the  introduction  to  demonstrate  the  technological 
importance  of  PZT.  These  applications  use  PZT  in  polycrystalline  ceramic  form,  and  thus 
the  properties  of  these  ceramics  are  well  established  in  the  literature.  However,  the 


mechanisms  contributing  to  these  outstanding  polycrystaliine  properties  of  PZT  are  not  well 
understood,  because  of  the  complexity  of  the  interactions  within  the  polycrystaliine  material. 

A  first  step  in  the  analysis  of  a  ferroelectric  polycrystaliine  material  is  to  separate  the 
intrinsic  and  extrinsic  contributions  to  the  properties.  The  intrinsic  contributions  result  from 
the  averaging  of  the  single-domain  single-crystal  properties,  while  the  extrinsic  contributions 
arise  from  the  interactions  at  grain  or  phase  boundaries  and  from  the  domain  wall  or  thermal 
defect  motions.  Unfortunately,  due  to  the  difficulty  of  growing  good  quality  single  crystals 
of  PZT,  very  little  single-crystal  data  is  available.  Thus  the  goal  of  this  project  was  to 
develop  a  thermodynamic  phenomenological  theory  to  calculate  the  single-domain  properties 
of  PZT.  This  theory  could  then  be  used  to  separate  the  intrinsic  and  extrinsic  contributions  to 
the  polycrystaliine  properties.  In  addition  there  are  several  other  applications  of  this  theory, 
which  will  be  in  the  last  paper  of  this  series.  ^  ^ 

A  two-sublattice  theory,  where  each  sublattice  has  a  separate  polarization,  was  used 
to  account  for  the  ferroelectric  and  antiferroelectric  phases  of  the  PZT  system.  An  additional 
order  parameter  was  also  included  to  account  for  the  tilting  of  the  oxygen  octahedra  in  the 
low-temperature  rhombohedral  phase.  The  resulting  energy  function  can  be  used  to  model 
the  phase  transitions  and  single-domain  properties  of  the  PZT  system.  Solutions  to  this 
energy  function  were  used  to  derive  relations  for  the  energies;  spontaneous  polarizations, 
strains,  and  tilt  angles;  and  dielectric  and  piezoelectric  properties  corresponding  to  the 
different  phases  in  the  PZT  system. 

The  coefficients  needed  to  calculate  the  energies  and  properties  of  the  solutions  will 
be  determined  in  the  next  three  papers  in  this  series.  1  ^  Theoretical  calculations  and 

comparisons  with  experimental  data  will  then  be  presented  in  the  final  paper  of  this  series.^ 
Additional  details  of  this  theory,  including  tables  of  all  the  experimental  data  collected,  can  be 

found  in  Reference  4 1 . 
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Table  I.  Coefficients  of  the  PZT  Energy  Function. 
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ferroelectric  dielectric  stiffnesses  at  constant  stress 
anriferroelectric  dielectric  stiffnesses  at  constant  stress 
coupling  between  the  ferroelectric  and  anriferroelectric  polarizations 
octahedral  torsion  coefficients 

coupling  between  the  ferroelectric  polarization  and  tilt  angle 
elastic  compliances  at  constant  polarization 

electrostrictive  coupling  between  the  ferroelectric  polarization  and  stress 
electrostrictive  coupling  between  the  anriferroelectric  polarization  and  stress 
rotostricrive  coupling  between  the  tilt  angle  and  stress 
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Two  tricritical  points,  where  a  phase  transition  changes  from  first  to  second  order, 
were  found  to  occur  in  the  lead  zirconare-titanate  (PZT)  solid  solution  system. 
High-temperature  x-ray  diffraction  data  on  sol-gel  derived  PZT  powders  were  used  to 
calculate  the  cell  parameters  and  spontaneous  strain  in  the  ferroelectric  phases.  These  data 
were  used  with  addition  data  from  the  literature  to  determine  values  of  the  higher-order 
dielectric  stiffness  coefficients,  which  were  then  used  to  locate  the  tricritical  points.  The 
values  of  the  coefficients  are  also  needed  in  the  development  of  a  thermodynamic  theory  of 
the  PZT  system. 
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I.  INTRODUCTION 

Tricriticai  behavior,  where  a  phase  transition  changes  from  first  to  second  order,  has 
been  shown  to  occur  tn  the  lead  zirconate-titanate  (PZT)  solid  solution  system.1'4  The  end 
members  lead  titanate  (PbTiC^)  and  lead  zirconate  (PbZrC^)  both  have  large  first-order 
phase  transitions  from  a  paraelectric  cubic  phase  at  high  temperatures  to  ferroeiectnc 
tetragonal  and  rhombohedral  phases,  respectively,  at  lower  temperatures.  By  forming  a  solid 
solution  between  these  fust -order  end  members  a  second-order  transition  region  develops  in 
the  middle  of  the  phase  diagram.4 

The  degree  of  first  order  behavior  has  been  shown  to  decrease  from  lead  zirconate  to 
the  PZT  94/6  (94%  PZ  and  6%  PT)  composition,  where  a  tricriticai  point  occurs  and  the 
transition  changes  to  second  order.  ^  With  increasing  titanium  content  the  second  order 
transition  region  was  found  to  extend  from  the  tricriticai  point  over  to  the  PZT  88/12 
composition.^  Since  lead  titanate  has  a  first  order  transition,  a  second  tricriticai  point  should 
occur  between  the  PZT  88/12  and  lead  titanate  compositions. 

To  provide  additional  data  to  locate  the  second  tricriticai  point,  pure  homogeneous 
sol-gel  derived  powders  were  prepared  for  several  PZT  compositions.^  The  lattice 
parameters  of  these  compositions  were  determined  from  high-temperature  x-ray  diffraction, 
and  used  to  calculate  the  spontaneous  strain.^  By  using  these  data  to  determine  the  higher 
order  dielectric  stiffness  coefficients,  the  second-order  transition  region  was  found  to  extend 
over  to  the  morphotropic  boundary  between  the  tetragonal  and  rhombohedral  phases.4 

In  this  paper  additional  data  will  be  presented  to  show  that  the  second  tricriticai  point 
occurs  in  the  tetragonal  phase  field.  The  compositional  dependence  of  the  higher-order 
dielectric  stiffness  coefficients  will  be  determined  from  high-temperature  x-ray  diffraction 
data,  and  used  to  locate  the  tricriticai  points.  Values  of  the  dielectric  stiffness  coefficients  are 
also  needed  in  the  development  a  thermodynamic  theory  to  model  the  phase  transitions  and 
single-domain  properties  of  the  entire  PZT  system.^'  ^ 
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In  the  next  section  experimental  high- temperature  x-ray  diffraction  data  will  be 
presented.  This  data  will  be  used  to  calculate  the  spontaneous  polarization  for  tetragonal  and 
rhombohedral  compositions  in  Sections  HI  and  IV,  respectively.  This  polarization  data  will 
then  be  used  to  determine  the  higher-order  dielectric  stiffness  coefficients,  and  to  establish 
the  location  of  the  tricntical  points.  In  Section  V  the  morphotropic  boundary  will  be  used  to 
complete  the  evaluation  of  these  coefficients,  followed  by  a  summary  of  this  paper  in  Section 
VI. 


H.  HIGH-TEMPERATURE  X-RAY  DIFFRACTION 

High-temperature  x-ray  diffraction  was  used  to  determine  the  cell  parameters  of 
sol-gel  derived  lead  titanate  and  PZT  powders.  The  lead  titanate  data  was  previously 
presented  in  Reference  1 1.  and  used  to  develop  a  thermodynamic  theory  of  lead  titanate. 
X-ray  data  on  PZT  90/10,  80/20,  70/30  and  60/40  were  previously  presented  in  Reference  6 
In  this  section  high-temperature  x-ray  diffraction  data  on  sol-gel  derived  powder  of  PZT 
32/68  will  be  presented.  The  sol-gel  procedure  used  to  prepare  this  powder  was  described  in 
Reference  5.  The  x-ray  diffractometer  setup  and  data  analysis  procedure  that  were  used  were 
described  in  References  6  and  1 1 . 

The  splitting  of  the  002/200  and  123/312/321  peaks  for  the  PZT  32/68  composition 
are  shown  in  Figure  1.  These  x-ray  data  were  used  to  calculate  the  cell  constants  using  the 
Cohen  least-squares  refinement  method^  as  shown  in  Figure  2.  The  data  in  these  figures 
indicate  that  the  transition  from  cubic  to  tetragonal  in  the  PZT  32/68  composition  is  either 
second  order,  or  only  slightly  first  order.  This  is  more  obvious  when  plotting  the  cell 
volume  (=  aj^cj)  versus  temperature  as  shown  in  Figure  3,^  A  continuous  change  in  the 
volume  would  indicate  that  the  transition  is  second  order. 
iikMijii'iiwniri  m  hit  fii'V  oHbi 


4a*  Additional  data  will  be  presented  later  tn  this  paper  to  show  that  the  cubic  to  tetragonal 


transition  in  PZT  32/68  is  probably  second  order. 

The  spontaneous  strains  x  j  and  X3  can  be  calculated  from  the  tetragonal  cell 
constants  ay  and  cy  using  the  following  relations: 

aT  -  a^'  CT  "  *C 

Xj  =  -  X3  =  -  (l) 

ac’  ac’ 

is  the  cubic  cell  length  ex  -apolated  into  the  tetragonal  region. 

Different  procedures  have  been  used  to  calculate  the  spontaneous  strains,  depending 
on  the  procedure  used  to  determine  a^’.  Haun  et  al.  ^  extrapolated  the  cubic  cell  constant 
into  the  tetragonal  region  by  assuming  that  the  electrostrictive  ratio  Qj  ^/Q 1 2  ^  independent 
of  temperature,  which  is  a  fairly  good  assumption  in  PZT.  ^  This  ratio  can  be  determined 
from  the  ratio  of  the  spontaneous  strains  xj  and  X3  from  Equation  (21)  in  Reference  7: 


xl  Ql2 


If  the  Qj  1/Q12  ratio  is  known,  then  a^’  can  be  determined  by  combining  Equations  (l )  and 

(2): 

CT  “  iQll/Ql2)aT 

ar'  = -  (3) 

l  -  Q11/Q12 

By  determining  the  Qj  j/Q  1 2  rati°'  aC>  can  ^  calculated  and  used  with  the  ceil 
constants  to  calculate  the  spontaneous  strains.  Haun  et  al. 1 1  used  a  procedure  of 
extrapolating  the  cubic  cell  constant  data  to  below  the  transition  to  calculate  the  Q  j  j/Q  \  2 
ratio  for  PbTiC>3.  Due  to  the  second  order  transition  behavior  of  the  PZT  32/68  composition, 
this  procedure  was  not  posssible.  However,  the  electrostrictive  constants  of  PZT  have  been 
recently  approximated  from  a  combination  of  single-crystal  and  polycrystalline  data.  *  ^  A 

t  "1 

value  of  -  3. 166  for  the  Q 1  [/Q  j  2  ratio  was  determined  for  the  PZT  32/68  composition. 

This  value  was  used  in  Equation  (3)  to  calculate  a^’  as  shown  by  the  data  plotted  in 
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Figure  2. 

The  spontaneous  strains  x  ^  and  xj  were  then  calculated  for  the  PZT  32/68 
composition  from  the  data  in  Figure  2  using  Equation  1.  The  results  are  shown  in  Figure  4, 
along  with  the  lead  titanate  strain  data  from  Reference  11.  The  numerical  values  of  this  data 
are  listed  in  Reference  14.  These  strain  data  will  be  used  in  the  next  section  with  the 
electrostrictive  constants  to  calculate  the  spontaneous  polarization,  which  will  then  be  used  to 
determine  the  higher-order  dielectric  stiffness  coefficients. 


HI.  TETRAGONAL  COMPOSITIONS 


All  of  the  dielectric  stiffness  coefficients  in  the  PZT  energy  function  [Equation  [1 )  in 
Reference  7]  were  assumed  to  be  independent  of  temperature,  except  the  dielectric  stiffness 
constant  which  was  given  a  linear  temperature  dependence  based  on  the  Curie-Weiss  law. 


2e0C 

C  is  the  Curie  constant,  Eq  is  the  permittivity  of  free  space,  and  Tq  is  the  Curie-Weiss 
temperature.  By  finding  values  of  Tq  and  C,  aj  can  be  calculated  as  a  function  of 
temperature. 

The  simplest  and  usual  method  of  finding  values  of  Tq  and  C  is  to  fit  the  inverse  of 
the  dielectric  constant  in  the  paraelectric  state  using  the  Curie-Weiss  law.  Unfortunately, 
very  little  dielectric  data  on  PZT  in  the  high  temperature  cubic  state  exists,  due  to  the  lack  of 
good  quality  single  crystals.  Measurements  on  polycrystalline  ceramic  samples  typically  give 
similar  values  of  Tq  and  C  as  from  single  crystals,  but  these  measurements  are  difficult  in 
pure  PZT  due  to  the  large  electrical  conductivity  at  high  temperatures.  In  addition, 
considerable  discrepancies  exist  in  the  measurements  that  have  been  made.  For  these  reasons 
indirect  methods  were  used  to  determine  Tq  and  C.  Tq  will  be  determined  by  fitting 
spontaneous  strain  data  as  described  later  in  this  section.  C  will  be  found  from 
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low-temperature  dielectnc  data  the  third  paper  of  this  series.8 

The  Curie  constant  acts  as  a  scaling  constant  in  the  dielectnc  properties  and  energies 
of  the  phases.  The  dielectric  stiffness  coefficients  are  inversely  proportional  to  the  Curie 
constant,  as  shown  by  Equation  ( 1)  for  aj.  However,  the  spontaneous  polarization,  strain, 
and  tilt  angle  are  ail  independent  of  the  Curie  constant.  Since  the  Cune  constant  was  not 
known  at  this  point,  new  constants  were  formed  by  multiplying  the  dielectric  stiffness 
coefficients  by  the  Cune  constant.  These  new  constants  can  be  determined  from  the 
spontaneous  polarization  data  independent  of  the  Curie  constant.  Equation  (4)  then  becomes: 


A 


T-T, 


a^C  = 


0 


2e0 


(5) 


Now  if  Tq  can  be  determined,  then  the  new  constant  aiC  can  be  calculated. 

The  spontaneous  polarization  of  the  tetragonal  state  can  be  related  to  the  dielectric 
stiffness  coefficients  by  solving  the  quadratic  relation  formed  from  the  first  partial  derivative 
stability  condition  [Equation  (14)  in  Reference  7]: 


P32  = 


-aH±[ali2-3a1alu] 


1/2 


(6) 


3a 


ill 


Only  the  +  solution  to  the  quadratic  equation  is  considered  here,  because  the  +  solution 
corresponds  to  a  free  energy  minimum,  while  the  -  solution  corresponds  to  a  maximum  of 
the  free  energy. 

Multiplying  the  numerator  and  denominator  of  Equation  (6)  by  the  Curie  constant 
results  in  the  following  relation  which  relates  P3  to  the  new  ajC,  a  j  [C.  and  j  jC 
constants: 


•  a.  ,C  +  [(a.  [Cl2  -  3  a,C  0|||C11/2 
*  - 
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(7) 


-i 


a 


I 


h 

<•> 

T*. 

f, 


1 


3 


rasasMa^^ 


,  p  _  » 

A.'t'J 


■'-‘A--- 


I  liFHWW  krtwovtKn  ww  K.uuuuumivwnAA  wwww’wwvjJWy^A WlArtWWJTWWJVWJVUUW 


H 


If  the  cubic -tetragonal  transition  is  second  order  then  Tq  =  Tq,  and  the  following 
relation  results  from  Equations  (5)  and  (7): 


P32  =  a  { L  -  [1  -  b  (T -TC)J1/2}, 

where  a=  -  cij  |C  / (3  aj  j  jC),  and  b=  3  04  1  jC  /  [2  £q  (aj  | C >2)  (8) 

The  a  and  b  relations  can  also  be  rearranged  in  terms  of  the  aj  and  j  jC  constants: 

allC  =  -  1/(2  eq  C  ab)  aj !  [C  =  1/(6  £q  C  a2  b)  (9) 

Since  jC  and  aj  j  jC  are  assumed  to  be  independent  of  temperature,  a  and  b  will  also  be 
independent  of  temperature.  Equation  (8)  shows  that  for  a  second  order  transition  the 
spontaneous  polarization  P3  goes  to  zero  at  the  transition  temperature  Tq. 

The  spontaneous  strains  x  j  and  *3  of  the  tetragonal  state  are  related  to  P3  through 
the  electrostrictive  coefficients  (see  Equation  (21 )  in  Reference  7].  By  substituting  Equation 
(8)  into  Equation  (21 )  from  Reference  7  the  following  relations  result: 

*1=  aQ12{l-[l-b(T-Tc)]1/2}  (10) 

*3=  aQn  (I -[1  -b(T-Tc)]1/2}  (11) 

These  equations  can  be  used  to  fit  the  tetragonal  strain  data,  if  the  cubic -tetragonal  transition 
is  second  order. 

If  the  cubic -tetragonal  transition  is  first  order,  then  Tq  *  Tq,  and  the  spontaneous 
polarization  changes  discontinuousiy  at  the  transition.  At  Tq  two  relations  must  be  satisfied: 

0  =  alC  +  al  l  +  al  1 1  (12) 

0  =  a1c  +  2a1[P3C2  +  3a111P3C4  ,  (13) 

where  ct[Q  and  P3Q  are  aj  and  P3  at  Tq.  Equation  (12)  was  derived  from  the  requirement 
that  the  AG’s  of  the  cubic  and  tetragonal  phases  [Equations  (8)  and  (9)  in  Reference  7]  must 
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be  equal  ax  T(;.  Equation  (13)  is  the  first  partial  derivative  stability  condition  [Equation  ( 14) 
in  Reference  7]  at  T^,  which  must  be  satisfied  so  that  the  stable  state  corresponds  to  the 
minima  of  the  energy  function. 


From  Equation  (4)  04  at  is: 


ctir  - 


TC-T0 

2  £q  C 


Substituting  this  equation  into  Equations  (12)  and  (13),  and  solving  for  the  temperature 


independent  coefficients  and  otju  results  in: 
-(TC-T0) 

all=  - a 

e0  C  P3C2 


TC‘T0 

:111  =  - - ——j 

2  ^  C  P3C4 


By  multiplying  these  relations  by  the  Curie  constant  C,  equations  result  for  the  new  constants 


a^Canda^jC: 


-  (TC  -  T0) 

«llc=  - — 2*  • 

*0  P3C2 


alllc  = 


TC-T0 

2  ^0  P3C 


Substituting  Equations  (16)  and  (5)  into  Equation  (7),  or  substituting  Equations  (13)  and  (4) 
into  Equation  (6),  results  in  the  following  relation: 

2  2  2  1  r  3  (T-T0) 

P32  =  yP3C2,  where  y  =  -(l+[l - J  1/2  }  (17) 


4(TC-T0) 


This  equation  can  be  used  to  calculate  the  spontaneous  polarization  P3  from  P3C,  Tq,  and 
TC  -  Tq,  when  the  paraelectric  -  ferroelectric  transition  is  first  order.  Comparing  this 
equation  with  the  second  order  P3  relation  [Equation  (S)]  shows  that  a  change  in  sign 
occurred  after  the  first  term.  This  is  because  for  a  first -order  transition  is  negative. 


Similar  relations  can  be  derived  for  the  spontaneous  strains  x^  and  X3  by 
substituting  Equation  (17)  into  Equation  (20)  from  Reference  7: 

*1  =  V*1C’  where  ^ic  =  Ql2^3C  (18) 

X3  =  y*3c,  where  X3^  =  Qj^  j  ^3C  (19) 

The  above  equations  were  used  to  fit  the  experimental  tetragonal  spontaneous  strain 
data  that  was  determined  from  high-temperature  x-ray  diffraction  data.  In  each  equation  there 
are  three  independent  unknown  constants  [x^,  Tq,  and  Tq  -  Tq  in  Equation  (18);  and  X3 q, 
Tq,  and  T(3  -  Tq  in  Equation  ( 19)1.  With  three  unknown  constants  many  combinations  of 
the  values  of  these  constants  will  give  similar  fits  of  the  experimental  data.  For  this  reason  a 
value  ofTc  was  first  determined  from  the  experimental  phase  diagram,  which  reduced  the 
number  of  unknown  constants  to  two.  With  only  two  unknown  constants,  the  combination 
of  values  that  gave  the  best  least  squares  fit  of  the  data  could  be  easily  found. 

*  Amin  et  al.  ^  fit  T^  of  the  experimental  phase  diagram  with  the  following 
polynomial  equation: 

Tc  =  (21 1 .8  +  486.0  x  -  280.0  x2  +  74.42  x3)  °C  ,  (20) 

where  x  is  the  mole  fraction  PbTiC^  in  PZT.  This  equation  will  be  used  to  calculate  T^ 
versus  composition  for  the  evaluation  of  the  dielectric  stiffness  constants  in  this  section,  but  a 
new  equation  [Equation  (42) j  will  be  given  in  Section  IV,  which  will  be  used  to  calculate  T^ 
in  References  8-10. 

Haun  et  al.  ^  used  a  computer  program  to  determine  values  of  x  x-^q,  and  Tq 
that  gave  the  best  least-squares  fit  of  the  lead  titanate  strain  data  (shown  in  Figure  4)  using 
Equations  ( 18)  and  ( 19)  with  T^  equal  to  492. 2°C  [calculated  from  Equation  (20)).  The 
electrostrictive  constants  of  lead  titanate  were  determined  using  Gavrilyachenko  et  als.  ^ 
room  temperature  spontaneous  polarization  value  of  0.75  C/m2,  and  used  with  the  strain  data 
to  calculate  the  spontaneous  polarization,  as  shown  in  Figure  5.  The  values  of  all  of  the 


constants  used  in  these  calculations  are  Listed  in  Table  I. 

The  same  procedure  was  used  to  fit  the  PZT  32/68  spontaneous  strain  data  from 
Figure  4.  However,  for  this  composition  the  least-squares  error  continually  became  smaller 
as  the  Tc  -  Tq  difference  was  reduced,  indicating  that  the  transition  was  actually  second 
order.  Therefore,  Equations  (10)  and  (11)  were  used  to  fit  this  data.  Values  of  the  constants 
aQi  i*  aQi2»  and  b  were  then  found  that  gave  the  best  fit  of  the  strain  data.  The 
spontaneous  polarization  was  calculated  by  combining  the  strain  data  with  the  electrostricuve 
constants  from  Reference  13.  The  experimental  polarization  data  and  theoretical  fit  of  the 
data  are  shown  in  Figure  5.  The  values  of  the  constants  used  in  these  calculations  are  listed 
in  Table  I. 

IV.  RHOMBOHEDRAL  COMPOSITIONS 

The  spontaneous  polarization  of  the  high-temperature  rhombohedral  state  can  be 
solved  from  the  quadratic  first  partial  derivative  stability  condition  [Equation  (16)  in 
Reference  TJ: 

-^K2- 9  a,  51 1/2 

p,2  =  - - -  ,  (21) 

3  % 

where  £  =  3  (a^  +  a12),  and  %  -  3  j  j  +  6  ctj  l2  +  al23  (22) 

As  in  the  tetragonal  case  the  numerator  and  denominator  of  this  equation  were  multiplied  by 
the  Curie  constant  C  to  relate  P3  to  the  new  constants: 

,  -;C  +  («;C)2-9alC$C]1''2 

P32  =  -  .  (23) 

3$C 

where  =  3  (at  jC  +  a12C  ).  and  =  3  0^  t  jC  +  6  04  j2C  +  a^C  (24) 


If  the  cubic -rhombohedral  transition  is  second  order  then  Tq  =  Tq,  and  the 
following  relation  results  from  Equations  (5)  and  (23): 


P32  =  a  |1  -[1  -b(T-Tc)J1/2). 


where  a  =  -  C,C  /  (3  qC).  and  b  =  9  qC  /  [2  E0  (CC)21 


The  a  and  b  relations  can  also  be  rearranged  in  terms  of  qC  and  ^C: 


=  -  3/(2  £0  C  a  b) 


qC  =  1/(2  £<)  C  a2  b) 


A  similar  relanon  can  be  found  for  the  spontaneous  strain  X4  by  substituting 


Equation  (23)  into  Equation  (24): 


x4  =  a  Q44  { 1  -  [1  -  b  (T  -  Tq)]*'2) 


As  in  the  tetragonal  case  at  Tq,  if  the  cubic-rhombohedral  transition  is  first  order, 
then  a  different  procedure  must  be  used  to  evaluate  the  coefficients.  At  Tq  the  energies  of 


the  cubic  and  rhombohedral  phases  must  be  equal  [Equations  (8)  and  (10)  in  Reference  7], 


and  the  first  partial  derivative  stability  condition  [Equation  ( 16)  in  Reference  7]  must  be 


satisfied: 


0  =  3  aic  +  ;  P3C2  +  %  P3C4  .  and 

o  =  alc  +  w;p3c2  +  5P3C4  , 


where  (Xjq  and  P3q  are  ctj  and  P3  at  Tq. 


Substituting  Equation  ( 14)  into  Equations  (28)  and  (29)  and  solving  for  the 


temperature  independent  coefficients  £  and  %  results  in: 


-  3  (Tq  -  T0) 


3  (Tq-Tq) 


e<)CP3C‘ 


2*0CP3C4 


By  multiplying  these  relations  by  the  Curie  constant  C,  equations  result  for  the  new  constants 
£C  and  %C: 


-3(Tc-T0)  ^  3(TC-T0) 

p3C2  2  ®0  P3C4 


(31) 


Substituting  Equations  (31)  and  (5)  into  Equation  (7)  or  substituting  Equations  (30)  and  (4) 
into  Equation  (6)  results  in  the  following  relation: 


2  2  2  I  r  3  (T'T())  1  1/2  ( 

Pi2  =  yPiC2.  where  y  =  -  (  1+  [  l - J  1/2  ) 

3  4(TC-T0) 


(32) 


This  equation  is  the  same  as  Equation  ( 17)  derived  for  the  spontaneous  polarization  of  the 
tetragonal  state,  except  that  P3  and  P^q  in  this  case  refer  to  the  rhombohedral  phase. 

A  similar  relation  can  be  derived  for  the  spontaneous  strain  X4  by  substituting 
Equation  (32)  into  Equation  (26): 


*4  *  V  X4C  •  where  *4C  =  044  P3C2 


(33) 


Either  Equation  (27)  or  (33)  was  used  to  fit  the  experimental  high-temperature 
rhombohedral  spontaneous  strain  data  from  Reference  6.  The  transition  temperature  T^  was 
first  determined  from  the  fit  of  the  experimental  phase  diagram  [Equation  (20)],  and  then 
values  of  the  remaining  two  unknown  constants  were  found  that  gave  the  best  least-squares 
fit  of  the  data.  The  best  fit  of  the  PZT  90/10  data  was  found  to  be  slightly  first  order,  while 
the  best  fits  of  the  PZT  80/20,  70/30,  and  60/40  compositions  were  second  order.  The  stain 
data  was  then  used  to  calculate  the  spontaneous  polarization  through  the  electrostrictive 
constants  from  Reference  13,  as  shown  in  Figure  6.  The  values  of  the  constants  used  in  the 
calculations  are  listed  in  Table  II. 
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V.  FITTING  THE  MORPHOTROPIC  PHASE  BOUNDARY 

In  the  last  two  sections  constants  involving  the  product  of  the  Cune  constant  and  the 
fourth  and  sixth  order  tetragonal  (ot[  jC  and  ctj  j  ^C)  and  rhombohedral  (£C  and  ^C) 
dielectric  stiffness  constants  were  determined  from  spontaneous  strain  and  electrostrictive 
data,  ai  jC  and  ^  jC  were  determined  for  PbTi03  and  PZT  32/68,  and  and  for 
PZT  90/10,  80/20.  70/30,  and  60/40.  In  addition  to  these  data,  values  of  and  qC  were 
calculated  for  PbZrOj  ^ .  A  value  of  was  also  determined  for  PbTiC^  * 

From  the  above  data  the  compositional  dependences  could  be  estimated  for  from 
PbZrOj  to  PbTiC>3,  for  qC  only  across  the  rhombohedral  phase  field,  and  for  jC  and 
aj  j  jC  only  across  the  tetragonal  phase  field.  Additional  data  were  therefore  needed  to 
determine  the  compositional  dependence  of  all  of  these  constants  across  the  entire  PZT 
system. 

The  £C  [=  3  (otj  a^)  C]  data  was  plotted  in  Figure  7  (c)  versus  composition. 
This  data  was  fit  with  the  following  equation: 

£C  =  [(a  +  bx)  exp^*  +  dx  +  el  10^  ,  (34) 

where  a  =  *  9.6,  b  =  -0.012501,  c  =  12.6,  d  =  0.42743,  e  =  2.6213,  exp  is  the  exponential 
function,  and  x  is  the  mole  fraction  of  PbTiC>3  in  PZT.  This  equation  was  then  used  to 
calculate  £C  versus  composition. 

The  rhombohedral  sixth-order  dielectric  stiffness  (%)  cannot  become  negative  when 
the  rhombohedral  phase  is  metastable,  and  the  stable  tetragonal  phase  undergoes  a 
second-order  transition  to  the  cubic  state.  The  [=  3  aj  j  6  aj  ^2  +  a123^  data  ^rom 
the  last  section  showed  that  this  constant  does  decrease  across  the  rhombohedral  phase  field 
towards  the  lead  titanate  composition,  and  may  level  off  without  becoming  negative.  To 
keep  this  constant  from  becoming  negative,  was  assumed  to  be  equal  to  j  jC  at  the 
lead  titanate  composition.  This  is  the  same  as  assuming  that  cx^C  is  equal  to  six  times 
al  12^-  ^or  PbTiC^.  This  data  point  was  plotted  in  Figure  8  (d)  along  with  the  rest  of  the  qC 


data,  which  was  then  fit  with  the  following  equation: 

*C  =  [(a  +  bx)  exp'cx  +  dx  +e|  10 ^  ,  (35) 

where  a  =  16.225,  b  =  -  0.088651,  c  =  21.255,  d  =  -  0.76973,  e  =  0.887,  and  x  is  the  mole 
fraction  of  PbTi03  in  PZT.  This  equation  was  then  used  to  calculate  versus 
composition. 

The  morphotropic  phase  boundary,  where  the  energies  of  the  tetragonal  and 
high-temperature  rhombohedral  phases  must  be  equal,  was  used  to  extrapolate  the  a  j  [  C  and 
<Xi  i  {C  constants  into  the  rhombohedral  phase  field.  The  energies  of  the  tetragonal  and 
high-temperature  rhombohedral  phases  were  defined  by  Equations  (9)  and  (1 1)  in  Reference 
7.  Multiplying  these  equations  by  the  Curie  constant  results  in  the  following  relations: 

Ft  AGC  =  alCP32  +  a11CP34  +  alllCP36  (36) 

fR(HT)  AC  C  =  3  a,C  P32  +  C  C  P34  ♦  \  C  P36  (37) 

Since  the  Curie  constant  is  the  same  for  the  tetragonal  and  rhombohedral  phases,  the  product 
of  AG  and  C  for  these  two  phases  should  be  equal  at  the  morphotropic  boundary. 

From  the  data  presented  in  the  last  two  sections,  the  paraelectnc-ferroelectnc 
transition  appears  to  be  second  order  on  both  sides  of  the  morphotropic  boundary,  and  thus 
T0  =  T(^.  The  relation  for  T^  [Equation  (20)]  was  therefore  used  to  calculate  Tq,  which  is 
needed  to  determine  04C  versus  composition  [see  Equation  (5)].  Using  this  procedure  to 
calculate  ajC  with  £C  and  ^C  determined  from  Equations  (34)  and  (35),  the  product  of  the 
energy  AG  and  the  Cune  constant  for  the  rhombohedral  state  could  be  calculated  using 
Equation  (37). 

A  computer  program  was  written  to  extrapolate  the  aj  jC  and  (X|  j  jC  constants  into 
the  rhombohedral  phase  field  by  fitting  the  morphotropic  boundary.  The  first  step  was  to 
make  initial  guesses  for  the  values  of  aj  jC  and  otj  j  jC  for  PbZi03.  A  quadratic  fit  of  this 
data  and  the  values  of  jC  and  ctj  |  ^C  determined  in  Section  in  for  the  PZT  32/68  and 


PbTi03  compositions  were  then  made  using  the  following  equation: 

a||C  or  otij[C  =  (a  +  bx  +  c^)  lO1^  (38) 

The  product  of  the  energy  AG  and  the  Curie  constant  of  the  tetragonal  state  could  now  be 
calculated  using  Equation  (36).  The  morphotropic  boundary  was  then  calculated  from  the 
cross  over  of  the  AG  C  products  of  the  two  phases,  and  a  least-squares  error  with  the 
experimental  data  was  calculated.  New  values  of  dj  jC  and  dj  y  jC  for  PbZrC^  were 
chosen,  and  the  procedure  repeated  until  the  best  least-squares  fit  of  the  morphotropic 
boundary  was  obtained.  The  final  fit  is  shown  in  Figure  2  of  Reference  10  by  the  solid 
curve.  The  values  and  final  fittings  of  the  d^C  and  d^C  constants  are  plotted  in  Figures 
7  and  8.  In  Equation  (38)  the  final  constants  for  dj  jC  were  a  =  10.612,  b  =  -22.655,  and  c 
=  10.955;  and  for  d^  [C  were  a  =  12.026,  b  =  -17.296,  and  c  =  9.1790. 

This  program  visually  demonstrated  how  the  shape  of  the  morphotropic  boundary 
could  be  changed  by  varying  the  dielectric  stiffness  coefficients.  Depending  on  the  values  of 
the  coefficients  chosen,  the  phase  boundary  could  be  shifted  either  towards  PbZr03  or 
PbTiC^,  and  the  shape  could  be  changed  from  vertical  to  curved.  The  bending  of  the 
boundary  could  also  be  made  to  go  towards  PbZrC>3  or  PbTiC^. 

The  d^C  constant  can  be  calculated  from  the  dj  yC  and  [=  3  (dj  j+  d^)  C] 
constants: 

di2C  =  ;c/3  -  dUC  (39) 

Using  this  relation  and  Equations  (34)  and  (38)  for  £C  and  dj  ^C,  d^C  was  calculated 
versus  composition  as  shown  in  Figure  7. 

When  dj  iC  or  C,C  change  sign  cricritical  points  occur,  and  the  cubic-tetragonal  or 
cubic-rhombohedral  transition  changes  from  first  to  second  order.  As  shown  in  Figure  7 
tricritical  points  (labeled  TCP)  occur  at  the  Pb(Zrj.xTix)03  compositions  with  values  of  x 
equal  to  0. 102  and  0.7 1 7.  Between  these  tricritical  points  the  paraelectric-ferroelectric 
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transition  is  second  order,  and  first  order  regions  exist  near  the  end  members  PbZr03  and 
PbTi03-  In  the  first  order  regions  the  Tc  -  Tq  difference  and  spontaneous  polarization  at  Tc 
(PjC)  can  be  calculated  from  the  following  relations  derived  from  Equations  (16)  and  (31): 

Ft  P3C  =  [-anC/(2a12C)l1/2.  Tc-T0  =  E0(allC)2/(2a12C)  (40) 

FR(HT)  P3C  =  [-sC/(2qC)i1/2,  Tc-T0  =  £0(CC)2/(6qC)  (41) 

Using  Equations  (34).  (35).  and  (38)  to  calculate  the  constants  in  these  equations.  P3C  and 
TC  -  Tq  were  calculated  as  shown  in  Figure  9. 

T(3  data  from  the  experimental  phase  diagram  were  used  with  the  above  -  Tq 
calculations  to  determine  Tq  data  in  the  first-order  regions.  These  data  were  then  combined 
with  the  T(3  Hata  from  the  second-order  region,  which  is  equal  to  Tq.  and  fit  with  the 
following  polynomial  equation: 

Tq  =  a  +  b  x2  +  c  x3  +  d  x4  +  e  x3  +  f  ,  (42) 

where  a  =  189.48.  b  =  843.40,  c  =  -  2105.5,  d  =  4041.8.  e  =  -  3828.3,  f  =  1337.8,  and  x  is 
the  mole  fraction  of  PbTiC^  in  PZT. 

T(3  was  then  calculated  from  the  following  relation: 

TC  =  (TC'T0)  +  T0>  (43) 

where  Tq  and  T^  -  Tq  were  calculated  from  Equations  (40)  -  (42).  T^  was  calculated  from 
this  equation  and  plotted  in  in  Figure  2  of  Reference  10  along  with  the  experimental  data. 
Equation  (42)  and  (43 )  were  used  to  calculate  Tq  and  T^  throughout  the  rest  of  this  series  of 
papers. 


VI.  SUMMARY 


High-temperature  x-ray  diffraction  data  on  sol-gel  derived  PZT  32/68  powder  was 
presented,  and  used  with  electrostrictive  data  to  calculate  the  spontaneous  polarization. 

These  data  were  used  with  previously  published  data  to  determine  the  values  of  the 
higher-order  dielectric  stiffness  coefficients  at  several  different  compositions.  The 
compositional  dependence  of  each  coefficient  was  determined  by  fitting  these  values  with  an 
equation.  A  set  of  equations  was  then  established  that  can  be  used  to  calculate  values  of  the 
coefficients  at  any  composition. 

All  of  the  coefficients  were  assumed  to  be  independent  of  temperature,  except  the 
dielectric  stiffness  coefficient  a^,  which  was  given  a  linear  temperature  dependence  based  on 
the  Curie- Weiss  law.  The  experimental  phase  diagram  was  extensively  used  in  the 
evaluation  of  the  coefficients,  by  requiring  that  the  energies  of  the  adjacent  phases  of  a  phase 
transition  to  be  equal  at  the  transition.  The  first  partial  derivative  stability  conditions  were 
also  used  as  additional  relations  in  the  evaluation  procedure. 

The  fourth-order  tetragonal  (<Xj  j)  and  rhombohedral  (aj  j  -t-  a^)  dielectric  stiffness 
coefficients  were  found  to  change  sign  as  a  function  of  composition,  indicating  that  two 
tricritical  points  occur  in  the  PZT  system,  where  the  paraelectric-ferroelectric  phase  transition 
changes  from  first  to  second  order.  By  extrapolating  between  the  data  points,  the  tricritical 
points  were  found  to  occur  at  Pb(Zrj.xTix)03  compositions  with  x  equal  to  0.102  and 
0.717.  The  transition  was  first  order  from  the  end  members  to  the  tricritical  points,  and  then 
a  large  second  order  region  occurs  across  the  phase  diagram  between  the  tricritical  points. 

The  tricritical  point  on  the  lead  zirconate  side  of  the  phase  diagram  was  previously 
found  to  occur  at  the  PZT  94/6  composition,^  four  percent  closer  to  lead  zirconate  than 
where  the  data  in  this  study  indicates  that  it  should  occur.  This  difference  may  simply  be  due 
to  the  experimental  error  involved  in  the  measurements  and  theoretical  fitting,  or  may  be 
related  to  differences  in  homogeneity  of  the  powders  used  (mixed-oxide  versus  sol-gel). 
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Haun  et  al."*  indicated  that  the  second  tricritical  point  occurred  near  or  possibly  at  the 
morphotropic  boundary.  With  the  more  recent  data  presented  in  this  paper  and  a  more 
careful  analysis  of  the  previous  data,  the  second  tricritical  point  appears  to  be  located  at  the 
PZT  28/72  composition  as  described  above.  Differential  Scanning  Calorimetric  (DSC)  data 
on  the  sol-gel  powders  prepared  in  this  project  also  indicated  that  the  second  tricritical  point 
occurs  near  the  PZT  28/72  composition.  The  shape  of  the  DSC  peaks  at  T^>  changed  from 
sharp  narrow  peaks  (first  order)  to  very  broad  peaks  (second  order)  at  this  composition. 
Additional  work  on  preparing  homogeneous  PZT  powders,  and  then  determining  the  lattice 
constants  as  a  function  of  temperature  from  high-temperature  x-ray  diffraction  is  needed  to 
locate  the  tricritical  points  more  precisely. 
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Figure  I.  The  angles  of  the  002/200  and  123/312/321  diffraction  peaks  of  PZT  32/68 
plotted  versus  temperature. 

Figure  2.  The  lattice  constants  of  PZT  32/68  plotted  versus  temperature,  a-p  and  cp  are  the 
lattice  constants  of  the  tetragonal  structure,  a^  and  a^'  are  the  lattice  constants 
of  the  cubic  structure  above  T^  and  extrapolated  into  the  tetragonal  region, 
respectively. 

Figure  3.  The  unit  cell  volume  of  (a)  PbTi03  and  (b)  PZT  32/68  plotted  versus  temperature. 

Figure  4.  The  spontaneous  strains  x  j  and  X3  plotted  versus  temperature  for  (a)  PbTiC^ 
and  (b)  PZT  32/68.  The  data  points  were  calculated  from  the  lattice  constant 
data  shown  in  Figure  2. 

Figure  5.  The  spontaneous  polarization  P3  of  PbTiC^  and  PZT  32/68  plotted  versus 
temperature.  The  data  points  were  calculated  from  the  experimental 
spontaneous  strain  data  plotted  in  Figure  4.  and  the  electrostrictive  constants 
from  Reference  16.  The  solid  curves  are  theoretical  fits  of  the  data. 

Figure  6.  The  spontaneous  polarization  P3  (the  component  of  P5)  of  the  high~temperarure 
rhombohedral  phase  for  the  PZT  90/10,  80/20,  70/30,  and  60/40  compositions 
plotted  versus  temperature.  The  data  points  were  calculated  from  the  experi¬ 
mental  spontaneous  strain  data  Reference  6,  and  the  electrostrictive 
constants  from  Reference  16.  The  solid  curves  are  theoretical  fits  of  the  data. 

Figure  7.  The  product  of  the  fourth-order  dielectric  stiffness  coefficients  and  the  Curie 
constant.  The  data  points  were  determined  from  the  experimental 
measurements.  The  solid  curves  are  theoretical  fits  of  the  data. 

Figure  8.  The  product  of  the  sixth-order  dielectric  stiffness  coefficients  and  the  Curie 
constant.  The  data  points  were  determined  from  the  experimental 
measurements.  The  solid  curves  are  theoretical  fits  of  the  data. 
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Figure  9.  The  spontaneous  polarization  at  (P3C)  [(a)  and  (b)j  and  the  Tq  -  Tq 

difference  [(c)  and  (d)j  plotted  versus  composition  in  the  first  order  regions  of 
the  phase  diagram. 
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Table  I.  Constants  used  to  Fit  the  Spontaneous  Strain  Data  and  to  Calculate  the 
Spontaneous  Polarization  Data  for  the  Tetragonal  Compositions. 


PbTiQ3 


TC(°C)  492.2 

T0(°C)  478.8 

TC-To(°C)  13.4 

XLC  ( 10-2)  -0.362 

x3c(10-2)  1.24 

Qu  (10-2m4/C2)  8.9 

Q12  (10*2  m4/C2)  -2.6 

P3C  (C/m2)  0.373 

a{  {C  (1012  m5  °C/C2F)  -  10.878 

al  i  lC  (101^  m9  °C/C4F)  3.909 
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Table  D.  Constants  used  to  Fit  the  Spontaneous  Strain  Data  and  to  Calculate  the 
Spontaneous  Polarization  Data  for  the  Rhombohedrai  Compositions 


PZT  90/10 


TC(°C) 

257.674 

T0(°C) 

257.655 

Tc  -  T0  CO 

0.019 

X4C  do*5) 

3.5328 

Q44(10-2  m4/C2) 

4.9 

P^c  (C/m2) 

0.0329 

;C(1014m5  °C/C2F) 

-  0.059390 

§C(1015m9aC/C4F) 

2.7458 

PZT  80/20 

PZT  70/30 

PZT  6040 

Tc  =  T0(oC) 

298.40 

334.41 

366.16 

a  Q44  (10"^) 

-2.2600 

-1.7487 

-5.4780 

b(10-2  0C-L) 

2.2808 

2.8542 

0.55338 

Q44  ( lO'2  m4/C2) 

5.17 

5.52 

6.29 

a  (lO'2  C2/m4) 

-4.3713 

-3.1679 

-8.7091 

;C(1014m5oC/C2F) 

1.6992 

1.8737 

1.9715 

§C(1015m9oC/C4F) 

3.5153 

3.5153 

1.3454 
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Values  of  the  Curie  constant  (C)  and  sixth-order  polarization  interaction  dielectric 
stiffness  coefficients  ((Xj  j2  and  0^23)  are  needed  for  the  development  of  a  thermodynamic 
theory  for  the  entire  lead  zirconare-titanate  (PZT)  solid  solution  system.  Low-temperature 
dielectric  data  measured  on  pure  homogeneous  polycrystalline  P2T  samples  were  used  to 
determine  values  of  these  coefficients  at  several  compositions  across  the  phase  diagram. 
Equations  were  then  fit  to  this  data  to  determine  the  compositional  dependence  of  the 
coefficients.  The  Curie  constant  was  found  to  form  a  peak  in  the  middle  of  the  phase 
diagram  at  the  PZT  50/50  composition. 
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I.  INTRODUCTION 


This  paper  is  the  third  paper  in  a  series  of  five  papers  *  describing  the  development 

of  a  thermodynamic  theory  for  the  entire  lead  zirconate-titanate  (PZT)  solid  solution  system. 
Values  of  the  Curie  constant  (C)  and  sixth-order  polarization  interaction  dielectric  stiffriess 
coefficients  (ctj  [2  and  <*123) 310  needed  for  the  development  of  this  theory. 

Due  to  the  lack  of  experimental  data,  the  Curie  constant  was  originally  assumed  to  be 
independent  of  composition  in  the  theory  developed  for  the  single-cell  region  of  the  PZT 
system.^  Amin  et  al.^  later  found  from  a  combination  of  calorimetric  and  phenomenological 
data  that  the  Curie  constant  was  dependent  on  composition  with  a  peak  forming  near  the 
morphotropic  boundary  between  the  tetragonal  and  rhombohedral  phases.  The  theory  was 
then  modified  to  account  for  the  compositional  dependence  of  the  Curie  constant.^  The 
Curie  constant  data  determined  by  Amin  et  al.^  extends  from  lead  titanate  to  the  morphotropic 
boundary  between  the  tetragonal  and  rhombohedral  phases.  Additional  data  were  needed  to 
complete  the  compositional  dependence  of  the  Curie  constant  from  the  morphotropic 
boundary  to  lead  zirconate. 

In  this  paper  values  of  C,  j2«  and  ®123  W^1  be  determined  as  a  function  of 
composition  from  low-temperature  polycry staiiine  dielectric  data.  In  the  next  section  the 
results  of  the  low-temperature  dielectric  measurements  will  be  presented.  These  data  will 
then  be  used  in  Section  m  to  determine  values  of  C,  (X[  j2-  and  0^23.  A  summary  of  this 
paper  will  be  presented  in  Section  IV. 


0.  LOW-TEMPERATURE  DIELECTRIC  MEASUREMENTS 

Additional  experimental  data  was  needed  to  determine  values  of  the  Curie  constant 
and  aj  1 2  and  23  coefficients.  Low -temperature  polycrystalline  dielectric  data  was  chosen 
for  this  purpose.  At  low  temperatures  the  thermally  activated  contributions  to  the 
polycrystalline  dielectric  properties  should  "freeze  out”.  The  remaining  dielectric  properties 


were  assumed  to  be  due  to  an  averaging  of  the  single-domain  properties. 

Zhuang  et  ai.^  fabricated  pure  homogeneous  polycrystalline  ceramic  PZT  samples 
from  sol-gel  derived  powders.  The  same  procedure  was  used  to  prepare  disc  shaped 
samples  for  use  in  this  study.  The  samples  were  cut  and  polished,  and  sputtered  with  gold 
electrodes.  The  samples  with  thermal-resistance  wire  attached  as  leads,  were  then  shielded  in 
a  copper  enclosure  in  an  Air  Products  and  Chemicals  model  LT-3-1 10  cryogenics  system  to 
cool  the  temperature  to  4.2  K.  The  dielectric  constant  and  loss  were  measured  at  1  KHz  on  a 
Hewlett  Packard  4270A  automatic  digital  capacitance  bridge. 

After  measuring  the  dielectric  properties  on  the  unpoled  samples  at  4.2  K,  the  same 
samples  were  poled  with  electric  fields  of  20  to  40  KV/cm  for  4  to  30  minutes.  The 
piezoelectric  strain  coefficient  djj  was  measured  using  a  Berlincourt  Piezo-d33  meter  to 
determine  the  completeness  of  poling.  The  poled  samples  were  then  cooled  back  down  to 
4.2  K,  and  the  dielectric  properties  were  remeasured.  The  unpoled  (£33),  and  poled  (£33  , 
parallel  to  the  poling  direction)  dielectric  constant  data  is  listed  in  Table  I. 

The  Bruggeman  formula^  was  used  to  relate  the  poled  (£33**)  and  unpoled  ( £33 ) 
dielectric  constants  of  the  tetragonal  and  rhombohedral  polycrystalline  samples  to  the 
single-domain  constants: 

£3/=  l/4(eu  -2fr(EU  +e33)+  {[en  -2  fr(en  +E33)]2  +  8eu£33l1/2)  0> 

£33=  l/4(en  +  |£h2  +  8  £i  1^33 1 l/2)  (2) 

Ell  and £33  are  the  dielectric  constants  perpendicular  and  parallel  to  the  poling  direction,  and 
d  is  the  fraction  of  90°  or  71°  (109°)  domain  alignment.  £j  j  and  £33  were  assumed  to  be 
equal  to  the  dielectric  susceptibilities  T)j  j  and  1)33  (actually  Ey  =  rjjj  +  1)  for  the  tetragonal 
compositions  and  equal  to  q  n'  and  1)33’  for  the  rhombohedral  compositions  [see  Equations 
(27)  and  (29)  in  Reference  l]. 

In  addition  to  the  poled  and  unpoled  ceramic  dielectric  constants,  the  fraction  of  90° 
or  71°(  109°)  domain  alignment  0  is  required  to  use  Equations  ( 1 )  and  (2)  to  calculate  the 


single-domain  constants.  Turik  et  al.9  used  an  x-ray  method  to  determine  0  as  a  function  of 
composition.  These  data  were  used  to  approximate  values  of  0  as  Listed  in  Table  l. 

Hie  differences  between  the  poled  and  unpoled  ceramic  dielectric  constants  at  low 
temperatures  (listed  in  Table  I)  showed  a  similar  behavior  as  was  seen  at  room 
temperature.9^®  By  poling  the  samples  the  dielectric  constant  decreased  for  the 
rhombohedral  compositions,  and  slightly  increased  for  the  tetragonal  compositions.  Turik  et 
al 9  attributed  these  changes  from  poling  as  being  due  to  the  reduction  of  the  clamping  effect 
from  the  180°  domains.  When  180°  domains  are  present,  the  dielectric  constant  is  lowerei 
from  a  clamping  effect  between  these  domains.  By  poling  a  ceramic  sample  virtually  all  of 
the  180°  domains  reorient  closest  to  the  poling  direction,  and  the  clamping  effect  is  reduo 
causing  the  dielectric  constant  to  increase. 

However,  90°  or  71°(109°)  domain  reorientation  also  occurs  during  poling,  which 
decreases  the  dielectric  constant  if  £33  is  less  than  Ej  j.  Thus  the  dielectric  constant  will 
increase  or  decrease  after  poling  depending  on  which  of  these  competing  mechanisms 
dominates.  In  the  tetragonal  state  the  small  fraction  of  the  90°  domains  that  realign  from 
poling  does  not  quite  cancel  out  the  increase  in  dielectric  constant  from  the  reduction  of  the 
clamping  effect  of  the  180°  domains,  and  therefore  the  dielectric  constant  increases  slightl; 
In  the  rhombohedral  state  a  large  fraction  of  the  71°  (109°)  domains  align  from  poling  and 
dominate  the  increase  from  the  reduction  of  the  clamping  effect,  causing  the  dielectric 
constant  to  decrease.  In  addition,  a  smaller  fraction  of  180°  domains  exist  in  a  rhombohedra 
ceramic  compared  to  a  tetragonal  one. 1 1 

Because  of  the  clamping  effect,  as  described  above,  Equations  (1)  and  (2)  were  only 
used  to  calculate  the  single-domain  dielectric  susceptibilities  n  1 1*  and  ^33’  for  the 
rhombohedral  PZT  94/6  through  54/46  compositions  using  the  data  from  Table  I  at  4.2  K  . 
This  data  is  listed  in  Table  Q.  The  single-domain  dielectric  susceptibilities  of  the  PZT  52/48, 
50/50,  and  40/60  compositions  will  be  calculated  using  a  different  procedure  in  the  next 
sectioa  The  single-domain  data  will  then  be  used  to  calculate  values  of  the  Curie  constant. 
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and  the  ct£  £2  and  a123  coefficients. 

m.  EVALUATION  OF  THE  CURIE  CONSTANT  AND  SIXTH-ORDER 
POLARIZATION  INTERACTION  DIELECTRIC  STIFFNESS  COEFFICIENTS 

In  the  second  paper  in  this  series^  the  compositional  dependences  of  the  jC, 
a i  £  £C,  CC.  and  %C  constants  were  determined.  These  constants  with  the  ot£C  constant  can 
be  used  to  calculate  the  spontaneous  polarization  and  AG  C  of  the  tetragonal  and  high 
temperature  rhombohedral  phases  [Equations  (7),  (23),  (36),  and  (37)  in  Reference  2].  In 
addition  to  these  constants,  the  Curie  constant  C,  and  the  polarization  interactions  coefficients 
a12*  al  12*  aiK*  <*123  are  required  to  calculate  the  second  derivative  properties,  such  as  the 
dielectric  properties  [see  Equations  (27)  and  (29)  in  Reference  1]. 

The  a^C  constant  was  calculated  from  the  ot£  £C  and  £C  [=  3  (a£  £+  ct£2>  C] 
constants  in  Reference  2.  A  combination  of  the  ct£  £2C  and  ot£23C  constants  can  be 
calculated  from  the  aj  £  £C  and  %C  [=  (3  ct£  £  £+  6  ct£  £  2  +  <X£  £  £ )  C]  constants  with  the 
following  relation: 

1 

(6 an2  +  ai23>C  =  | C  -  3ajj£C  (3) 

Using  this  equation  allowed  (6  aj  £2  +  ot£  23)  C  to  be  calculated  from  the  values  of  the  %C 
and  «£  £  £C  that  were  determined  in  Reference  2,  but  additional  data  were  still  needed  to 
separate  the  ot£  £2  C  and  <*123  C  constants. 

The  low-temperature  dielectric  data  presented  in  the  last  section  were  used  to 
calculate  the  Curie  constant  and  to  separate  the  dielectric  stiffness  coefficients  «£  12  a™* 
a123-  T°  accomplish  this  the  high-temperature  rhombohedral  equations  were  used  at  low 
temperatures,  where  actually  the  low-temperature  rhombohedral  phase  is  stable,  because  of 
the  following  reasons. 
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To  calculate  the  dielectric  susceptibility  coefficients  of  the  low-temperature 
rhombohedral  phase,  the  polarization-tilt  angle  coupling  coefficients  must  be  determined,  in 
addition  to  the  dielectric  stiffness  coefficients  [see  Equation  (30)  in  Reference  1  j.  The 
dielectric  properties  have  been  experimentally  shown  to  only  change  very  slighdy  at  the 
transition  between  the  high  and  low  temperature  rhombohedral  phases.  ^  A  fairly  good 
assumption  can  then  be  that  the  dielectric  susceptibility  coefficients  of  the  high  and  low 
temperature  rhombohedral  phases  are  equal  at  the  transition  between  these  phases.  Using 
this  assumption  with  the  methods  described  in  this  paper  and  the  second  and  fourth  papers  of 
this  series  to  determine  the  other  coefficients,  the  constants  needed  to  calculate  the  dielectric 
susceptibilities  of  die  low-temperature  rhombohedral  phase  can  be  solved  for. 

However,  when  using  the  low-temperature  rhombohedral  dielectric  susceptibility 
equations  and  die  experimental  ceramic  data  at  low  temperatures  in  the  present  theory  (tilt 
angle  coefficients  are  independent  of  temperature  and  only  go  up  to  the  fourth  order),  the 
resulting  values  of  the  Curie  constant  will  not  agree  with  the  available  experimental  data.  A 
temperature  dependence  was  added  to  the  second  order  tilt  angle  related  coefficient  [(3^  in  j 

Equation  (1)  of  Reference  1],  but  this  still  did  not  resolve  the  problem.  Probably  what  is 
needed,  in  addition  to  this  temperature  dependence,  is  to  add  a  sixth-order  tilt  angle  term  to  < 

the  energy  function.  But  if  this  is  done  the  equations  for  the  spontaneous  polarization  and  tilt 
angle  will  change  from  quadratic  to  quartic.  At  this  point  the  additional  complexity  resulting 
from  a  sixth-order  tilt  angle  term  is  probably  not  warranted,  and  there  is  not  enough 
experimental  data  available  to  properly  determine  these  additional  constants.  ] 

The  question  now  is  how  can  the  low -temperature  ceramic  dielectric  data  be  used  to 
determine  the  Curie  constant  and  to  separate  the  sixth  order  polarization  interaction 
coefficients,  if  the  low-temperature  rhombohedral  equations  will  not  give  reasonable  results 
with  the  present  theory.  The  dielectric  properties  of  the  high  and  low  temperature 
rhombohedral  phases  were  experimentally  found  to  be  very  similar  at  the  transition  between 
these  phases.  If  the  dielectric  properties  of  these  two  phases  remain  similar  down  to  low  ' 
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temperatures,  then  the  high-temperature  rhombohedral  dielectric  susceptibilities  relations 
could  be  used  at  low  temperatures. 

Using  the  high-temperature  rhombohedral  relations  at  low  temperatures  produced 
very  reasonable  results.  A  Curie  constant  of  2  x  10^  °C  was  obtained  for  the  PZT  90/10 
composition,  which  is  in  excellent  agreement  with  experimental  single-crystal 
measurements.1  ^  Thus  this  indicates  that  the  high  and  low  temperature  rhombohedral 
dielectric  properties  are  probably  similar  down  to  even  very  low  temperatures  (4.2  K).  The 
high-temperature  rhombohedral  dielectric  susceptibility  relations  were  therefore  used  for  the 
dielectric  properties  of  both  the  high  and  low  temperature  rhombohedral  phases,  which 
appears  to  be  a  fairly  good  assumption. 

The  Tin'  and  1133'  data  for  the  rhombohedral  PZT  94/6  through  54/46  compositions 
were  used  to  calculate  the  Curie  constant,  and  to  separate  the  (X|  12  and  a  ^3  constants  from 
the  following  high-temperature  rhombohedral  dielectric  stiffness  relations: 

nii'  =  C/{E0[(T-T0)/e0+l2auCP32  +  (30aniC+12a112C-2al23C)  P34]| 

(4) 

T133*  =  C/I  eo  [CT  -  T0)  /  eo  +  4  ;c  P32  +  10  $C  P34] }  (5) 


These  equations  from  derived  from  Equations  (29),  (33),  and  (40)  from  Reference  1. 

Equation  (5)  was  used  to  calculate  the  Curie  constant  from  the  1133'  data  in  Table  II. 
and  from  values  of  the  £C  and  %C  constants  from  Equations  (34)  and  (35)  from  Reference  2. 
Equations  (4)  and  (5)  were  then  combined  to  calculate  the  aj  j2C  and  C423C  constants 
using  Equations  (38)  and  (42)  from  Reference  2  to  determine  values  of  Tq,  a^C,  and 
aj  j  jC.  The  C,  aji2C,  and  aj^C  data  are  listed  in  Table  II. 

The  04  |2C  data  listed  in  Table  II  with  the  lead  titanate  value  from  Reference  13 
were  fit  with  the  following  relation: 

ot  1 1 2C  =  a  exp'15*  +  cx  +  d  ,  (6) 
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where  a  =  58.804,  b  =  29.397,  c  =  -3.3754,  d  =  4.2904,  exp  is  the  exponential  function, 
and  x  is  the  mole  fraction  PbTi03  in  PZT. 

The  compositional  dependence  of  the  a^C  constant  was  determined  from  the 
following  relation  derived  from  equation  (3): 

a123c  =  '  3alllc  *  6a112’ 

where  and  j  were  determined  from  Equations  (35)  and  (38)  from  Reference  2,  and 
al  12^  was  determined  from  Equation  (6)  above.  The  data  from  Table  II  and  the 
calculations  from  Equations  (6)  and  (7)  are  shown  in  Figure  1. 

The  Curie  constants  of  the  PZT  52/48  and  50/50  compositions  were  determined  by 
substituting  Equations  (4)  and  (5)  into  Equation  (1),  and  using  the  data  from  Table  I, 
along  with  the  dielectric  stiffness  constants  calculated  from  the  previously  given  equations. 
Values  of  T|n'  and  1333'  were  then  calculated  from  Equations  (4)  and  (5)  and  are  listed  in 
Table  n.  This  procedure  was  also  used  for  the  tetragonal  PZT  40/60  composition  using 
Equations  (27)  and  (36)  from  Reference  1  for  the  tetragonal  dielectric  susceptibilities.  The 
values  of  C,  and  ^33  are  listed  in  Table  II. 

The  best  least-squares  fit  of  the  following  relations  were  used  to  fit  the 
compositional  dependence  of  the  Curie  constant  using  the  data  listed  in  Table  n  along  a  value 
of  1.5  x  105  °C  for  PbTi03: 

For  x  £  0.5  :  C  =  [2.1716/(1  +  500.05  x)2  +  0.131  x  +  2.01]  105 

For  x£  0.5:  C  =  [2.8339/(1  +  126.56  x)2  +  1.4132]  105  (8) 

The  Curie  constant  data  and  fitting  from  this  equation  are  plotted  in  Figure  2.  Amin  et  al7 
also  found  a  peak  to  occur  in  the  Curie  constant  at  the  PZT  50/50  composition,  although  their 
peak  value  was  somewhat  larger  than  the  results  presented  here  (7.7  versus  4.3  x  10^  °C). 
Tsuzuki  et  al.  ^  experimentally  found  Curie  constant  values  of  5.8  and  6.0  x  10^  °C  from 
single-crystal  dielectric  measurements  on  PZT  50/50  and  51/49  compositions. 
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The  dielectric  <jata  listed  in  Table  II  are  plotted  in  Figure  3,  along  with  the  theoretical 
calculations  using  the  constants  determined  in  this  paper.  The  cross  over  of  the  q  1 1 '  and 
1133  coefficients  on  the  rhombohedral  side  of  the  diagram  resulted  from  the  extrapolation  of 
the  constants  between  the  compositions  where  the  data  points  occur,  and  may  or  may  not  be 
a  teal  effect  in  the  material.  In  any  case  the  dielectric  anisotropy  becomes  very  small  in  this 
region.  A  more  detailed  discussion  of  the  dielectric  anisotropy  will  be  provided  Reference  4. 

IV.  SUMMARY 

The  dielectric  properties  were  measured  on  unpoled  and  poled  polycrystalline 
samples  at  4.2  K.  These  data  were  used  to  calculate  the  single-domain  dielectric 
susceptibility  coefficients  using  the  Bruggeman  relation.  Values  of  the  Curie  constant  and 
<Xn2  ai23  coefficients  were  calculated  from  the  single-domain  data.  Equations  were 
then  used  to  fit  these  values  to  determine  the  compositional  dependence  of  the  coefficients  for 
use  in  the  development  of  a  thermodynamic  theory  for  the  PZT  system.  As  a  summary  of 
this  paper  and  the  second  paper  of  the  series,*"  Table  HI  gives  values  of  the  coefficients  of  the 
.energy  function. 

The  Curie  constant  was  found  to  form  a  peak  at  the  PZT  50/50  composition,  similar 
to  the  compositional  dependence  of  the  electrostrictive  constants.  ^  This  type  of  behavior  is 
in  agreement  with  earlier  data  in  the  literature,' ^>7.14  ^  point  it  is  not  clear  whether 

the  peak  is  occuring  due  to  the  motphotropic  boundary  between  the  tetragonal  and 
rhombohedral  phases  or  due  to  some  other  effect.  The  Curie  constant  is  a  constant  from  the 
high-temperature  cubic  state,  and  thus  should  probably  not  be  related  to  the  morphotropic 
boundary.  By  studying  the  compositional  dependence  of  these  constants  in  other  solid 
solution  systems  where  a  morphotropic  boundary  occurs  away  from  the  50/50  composition 
[e.g.  Pb(Mg  1/3^2/3)03  -  PbTi03  systemj  may  help  to  resolve  some  of  these  questions. 


„  _ _ _ 


REFERENCES 


*M.  J.  Haun,  E.  Furman,  S.  J.  Jang,  and  L.  E.  Cross,  "Thermodynamic  Theory  of 

the  Lead  Zirconate-Titanate  Solid  Solution  System,  Part  I:  Phenomenology," 

J.  Appl.  Phys.  (submitted). 

J.  Haun,  E.  Furman,  H.  A.  McKinstry,  and  L.  E.  Cross,  "Thermodynamic 

Theory  of  the  Lead  Zirconate-Titanate  Solid  Solution  System,  Part  II:  Tricriticai 
Behavior,"  J.  Appl.  Phys.  (submined). 

-*M.  J.  Haun,  E.  Furman,  T.  R.  Halemane,  and  L.  E.  Cross,  "Thermodynamic  Theory  of 
the  Lead  Zirconate-Titanate  Solid  Solution  System,  Part  IV:  Tilting  of  the  Oxygen 
Octahedra,"  J.  Appl.  Phys.  (submitted). 

^M.  J.  Haun,  E.  Furman,  S.  J.  Jang,  and  L.  E.  Cross,  "Thermodynamic  Theory  of 
the  Lead  Zirconate-Titanate  Solid  Solution  System,  Part  V:  Theoretical 
Calculations,"  J.  Appl.  Phys  (submitted). 

^A.  Amin,  Phenomenological  and  Structural  Studies  of  Lead  Zirconate-  Titanate 

Piezoceramics,  Ph.D.  Thesis  in  Solid  State  Science,  The  Pennsylvania  State 
University  (1979). 

^A.  Amin,  L.  E.  Cross,  and  R.  E.  Newnham,  Ferroelectrics  65,  107  (1985). 

^A.  Amin,  M.  J.  Haun,  B.  Badger,  H.  A.  McKinstry,  and  L.  E.  Cross,  Ferroelectrics  65, 
107(1985). 

^Z.  Q.  Zhuang,  M.  J.  Haun,  S.  J.  Jang,  and  L.  E.  Cross,  Am.  Ceram.  Soc.,  (submitted). 

^A.  V.  Turik,  M.  F.  Kupriyanov,  E.  N.  Sidorenko,  and  S.  M.  Zaitsev,  Zh.  Tekh.  Fiz.  50. 
2146  (1980)  [Sov.  Phys.  Tech.  Phys.  25,  1251(1980)]. 

*®Z.  Q.  Zhuang,  M.  J.  Haun,  S.  J.  Jang,  and  L.  E.  Cross,  Proceedings  of  the  6th  IEEE 

International  Symposium  on  the  Applications  of  Ferroelectrics,  Lehigh,  PA.  edited 
by  Van  Wood  (IEEE,  New  York,  1986)  pp.  394-397. 

*  *D.  Berlincourt  and  H.  A.  Krueger,  J.  Appl.  Phys.  30.  1804  (1959). 


1 1 

l2R.  W.  Whatmore,  R.  W.,  R.  Garke,  and  A.  M.  Glazer,  J.  Phys.  C:  Solid  State  Phys.. 
11,3089(1978). 

J.  Haun,  Thermodynamic  Theory  of  PbTiO^,  J.  Appl.  Phys.  62,  3331  (1987). 
l4K.  Tsuzuki,  K.  Sakata,  and  M.  Wada,  Ferroelectrics,  8,  501  ( 1974). 

15m.  J.  Haun,  Z.  Q.  Zhuang,  E.  Furman,  S.  J.  Jang,  and  L.  E.  Cross,  Am.  Ceram.  Soc. 
(submitted). 


UST  OF  FIGURES  AND  TABLES 


1  2 


Figure  1. 


Figure  2. 


Figure  3. 
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Product  of  the  sixth-order  polarization  interaction  dielectric  stiffness  coefficients 
and  the  Curie  constant  plotted  versus  composition.  The  data  points  were 
determined  from  experimental  data.  The  solid  curves  are  fits  of  the  data  using 
Equations  (6)  and  (7). 

The  Curie  constant  plotted  versus  composition.  The  data  points  are  from 
Table  Q.  The  solid  curve  is  a  fitting  of  the  data  using  Equation  (8 ). 

The  dielectric  susceptibility  coefficients  r|  y  y  and  TI33  in  the  tetragonal  state  and 
rju'  and  ^33'  in  the  rhombohedral  state  plotted  versus  composition  at  4.2  K. 
The  data  points  are  from  Table  II,  except  for  the  points  at  PbTiC>3  which  are 
from  Reference  13.  The  solid  curves  are  the  theoretical  calculations  after 
determining  the  compositional  dependence  of  the  coefficients. 
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Table  I.  The  Dielectric  Constant  at  4.2  K  on  Poled  and  Unpoled  Ceramic  Samples  and  the 
Fraction  of  90°  (or  71,  109°)  Domain  Alignment  d. 

Table  II.  Calculations  from  the  Low  Temperature  Ceramic  Dielectric  Data. 

Table  III.  Values  of  the  Coefficients  of  the  Energy  Function  at  25  °C. 


Table  I.  The  Dielectric  Constant  at  4.2  K  on  Poled  and  Unpoled  Ceramic  Samples 
and  the  Fraction  of  90°  (or  7 1 ,  109°)  Domain  Alignment  0. 


Zr/Ti 

e33P 

e33 

O 

94/6 

93.2 

87.2 

0.65 

90/10 

106.2 

101. 1 

0.65 

70/30 

139.4 

135.3 

0.65 

60/40 

173.2 

156.2 

0.63 

54/46 

282.6 

232.5 

0.50 

52/48 

324.3 

303.0 

0.34 

50/50 

326.8 

328.4 

0.30 

40/60 

200.6 

202.0 

0.23 

V 
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Table  II.  Calculations  from  the  Low  Temperature  Ceramic  Dielectric  Data. 

all2C  a123C 

(1014  m9  °C/  C4  F)  (1015  m9  °C/  C4  F) 


Zr/Ti  TJn' 

°r  Hu 

Rhombohedral 

^33' 
or  n33 

C 

(105oC) 

94/6 

98.0 

84.1 

2.0031 

90/10 

110.2 

98.4 

2.0823 

70/30 

142.6 

133.1 

2.1648 

60/40 

187.5 

146.8 

2.4243 

54/46 

338.4 

189.4 

3.1714 

52/48 

392.3 

243.3 

4.0965 

50/50 

Tetragonal 

429.7 

254.1 

4.2962 

13.905 

5 

-3.3019  3 

7.0625 

-  1.8033  0 

3.3549 

-  1.5592  =£ 

3.1431 

-  1.5008  | 

2.4256 

-  1.1019  $ 

<3 
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Table  QL  Value*  of  the  Coefficients  of  the  Energy  Function  at  23  °C. 

Mole  Fraction  PbTiOj  in  PZT 


0.0 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

1.0 

Tc(«0 

231.5 

256.5 

300.6 

334.4 

364.3 

392.6 

418.4 

440.2 

459.1 

477.1 

492.1 

cao5  °c) 

2.027 

2.050 

2.083 

2.153 

2.424 

4.247 

2.664 

1.881 

1.642 

1.547 

1.500 

Qu(10*2Qta2) 

4.620 

5.080 

5.574 

6.175 

7.260 

9.660 

8.116 

7.887 

8.142 

8.504 

8.900 

Ql2(10‘2Qta2) 

-  1.391 

-  1.540 

-  1.720 

-  1.997 

-2.708 

-4.600 

-2.950 

-2.480 

-2.446 

-2.507 

-2.600 

Q^dO^Ota2) 

4.664 

4.900 

5.165 

5.522 

6.293 

8.190 

6.710 

6.356 

6.417 

6.569 

6.750 

Oj  (107m/F)al23oC 

-4.582 

-6.376 

-  7.470 

-8.116 

-7.904 

-4.887 

-8.340 

-  12.47 

-  14.84 

-  16.17 

-  17.08 

ond07m5A^F> 

32.35 

41.25 

31.29 

2130 

13.62 

4.764 

3.614 

0.6458 

-3.050 

-  5.845 

-  7.253 

Ojjdfl^mVF) 

- 16.71 

-4.222 

-  0.0345 

1.688 

1391 

1.735 

3.233 

5.109 

6.320 

7.063 

7.500 

^lO8  m^F) 

-34.42 

-0.2897 

9.284 

11.75 

11.26 

6.634 

10.78 

15.52 

18.05 

19.44 

20.32 

alud^mVF) 

5.932 

5.068 

4.288 

3.560 

1713 

1.336 

1.859 

1348 

2.475 

2.518 

1606 

oll2(10?m9yC4F> 

311.2 

34.45 

18.14 

15.27 

12.13 

6.128 

8.503 

10.25 

9.684 

8.099 

6.100 

Ojy  ( 109  m?/C*  F) 

-  104.1 

-8.797 

-7.545 

-  7.052 

-5.690 

-2.894 

-4.063 

-5.003 

-4.901 

-  4.359 

-  3.660 

4(10^  m9/C*F) 

84.41 

13.39 

4.627 

3.176 

2.402 

1.183 

1.596 

1.851 

1.652 

1.256 

0.7818 

PROOUCT  OF  SIXTH  ORDER  DIELECTRIC 


i* 


juuMMJiHwiMMOuoLyTiyiiwMiuataBWV^v^ui^^uvirwi’Tuvvituvu'vw^ 


THERMODYNAMIC  THEORY  OF  THE  LEAD  ZIRCONATE-TITANATE  SOLID 
SOLUTION  SYSTEM,  PART  IV:  TILTING  OF  THE  OXYGEN  OCTAHEDRA 

M J.  Haun,t  E.  Furman, ft  T.R.  Halemane.ttt  and  L.E.  Cross 

Materials  Research  Laboratory 
The  Pennsylvania  State  University 
University  Park,  PA  16802 

Values  of  coefficients  related  to  the  tilting  of  the  oxygen  octahedra  in  the  low-temperature 
rhombohedral  phase  of  the  lead  zirconate-titanate  (PZT)  solid  solution  system  are  needed  in  the 
development  of  a  thermodynamic  theory  of  the  PZT  system.  In  this  paper  these  coefficients  will  be 
determined  from  experimental  spontaneous  strain  and  oxygen  octahedral  tilt  angle  data. 

Values  of  the  tilt  angle  related  coefficients  Pi,  (in  and  <[>  are  assumed  to  be  independent  of 
temperature  and  are  then  related  to  Ph,  Pl  and  Tr  where  Ph  is  the  spontaneous  polarization  in  the 
untilted  phase  and  Pl  the  polarization  in  the  tilted  phase  at  the  transition  temperature  Tr  between 
tilted  and  nontilted  structures. 

Pl  and  0l  are  derived  from  experimental  data  for  the  PZT  90/10  composition,  and  used 
together  with  spontaneous  strain  values  at  this  composition  to  determine  electrostriction  (Q44)  and 
rotostriction  R44  constants.  These  constants  are  then  in  mm  used  to  derive  Pl  and  0l  values  for 
other  PZT  compositions  from  the  spontaneous  strain  vs  temperature  at  each  composition,  and  in 
the  final  paper  in  this  sequence  to  derive  the  Pi,  Pi  1  and  <D  values  for  these  compositions. 


tNow  at  DuPont  Experimental  Station,  Electronics  Dept.,  Bldg.  334,  Wilmington,  DE  19898. 
^Now  at  Allied-Signal  Inc.,  Morristown,  NJ  07960. 
tttNow  at  AT&T  Bell  Laboratories,  Holmdel,  NJ  07733. 
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I.  INTRODUCTION 

This  paper  is  the  fourth  paper  in  a  series  of  five  papers^  describing  the 
development  of  a  thermodynamic  theory  for  the  entire  lead  zirconate-titanate  (PZT)  solid 
solution  system.  In  the  low-temperature  rhombohedral  phase  in  the  PZT  system  the  oxygen 
octahedra  tilt  or  rotate  about  the  polar  axis.  In  the  first  paper  of  this  series  the  oxygen 
octahedral  tilt  angle  was  included  in  the  PZT  energy  function  as  an  order  parameter  [see 
Equation  ( l )  in  Reference  l  j  to  account  for  the  low-temperature  rhombohedral  phase.  In  this 
paper  the  hit  angle  related  coefficients  of  the  energy  function  will  be  determined. 

n.  EVALUATION  OF  THE  TILT  ANGLE  RELATED  COEFFICIENTS 


The  tilt  angle  related  coefficients  will  be  determined  from  the  available  experimental 
low-temperature  rhombohedral  spontaneous  strain  and  tilt  angle  data.  Multiplying  the  energy 
of  the  low -temperature  rhombohedral  state  [Equation  (12)  in  Reference  l  ]  by  the  Curie 
constant  results  in  the  following  relation: 


Fr(LT)  ao  c  =  3  <X,C  p32  +  ;  c  p34  +  %  c  p36  +  3P(c  s32 

+  3PnCe34  +  3<|)CP32932. 

where  $  C  =  (Y(  i  +  2  Yi2  +  Y44)  C 


(1) 

(2) 


£  C  and  \  C  were  defined  by  Equation  (24)  in  Reference  2. 

Equations  for  the  spontaneous  polarization  P3  and  tilt  angle  03  can  be  found  from 
the  first  partial  derivative  relations  [Equations  ( 1 7)  and  ( 18)  in  Reference  1  ]: 

b  +  [b2  -  4  a  c)^ 


P72  and  9-j2  = 


2a 


where  for  P-»2: 


u  =  ;c 

b  =  2/3  ;c  -«DC)2/(2  0llC) 
c  =  ot,  C-p,C#C/(2pllC) 


(3) 


’■>  V*.  * 


and  for  63^: 


a  =  4$C(PuC)2/(4>C)2 
b  =  0  C  +  4/3  Pj  {0/(0  C)  [3  %  C  ptC  A<|>  O  -  ;  Cl 
c  =  04  C  +  PiC/(^  C)  £  C  ptc  /«>  C)  -  2/3  ;  C] 
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(4) 


Values  of  the  ajC.  ^  C.  and  $  C  coefficients  were  determined  in  the  second  paper  in 
this  series.  Thus  the  tilt  angle  related  coefficients  PjC.  Pj  jC,  and  0  C  are  the  only  unknown 
coefficients  in  these  equations.  These  coefficients  were  assumed  to  be  independent  of  ! 
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temperature,  and  were  related  to  three  new  constants  at  the  high  to  low  temperature  j 

rhombohedral  transition  temperature.  Tr  ,  using  the  first  partial  derivative  stability 
conditions  [Equations  (17)  and  ( 18)  in  Reference  1],  and  by  equating  the  energies  of  the  two 
phases  at  the  boundary  [Equations  ( l )  above  and  Equation  (37)  in  Reference  2).  This 
resulted  in  the  following  equations: 


♦  C  =  -  1/0L2  [2  (3  04  tC  +  al2C)  (PH2  -  PL2) 

+  (15  omC+  14  a^2C  +  <*i23C)  (P^4  *  Pl4)1  (5) 

PlC  *  -  *C  PL2  +■  1/9l2  [(Tr  -  T0)/eo  (PH2  -  PL2)  +  2/3  CC  (PH4  -  PL4) 

+-2/3  5C(PH6-PL6)]  (6) 

Pt  iC  =  1/0L4  [(Tr  -  T0)/(2  E0)  (PH2  -  PL2)  +  CC/3  (PH4  -  PL4) 

+  %Q/3  (PH6  -  PL6)  -  PjC  0L2  -  <pC  PL29l21  (7) 


Pj^  and  Pl  are  the  spontaneous  polarization  of  the  high  and  low  temperature  rhombohedral 
states  at  Tr.  0^  is  the  tilt  angle  of  the  low-temperature  rhombohedral  state  at  Tr. 

Pp{  can  be  calculated  from  Equation  (23)  in  Reference  2  using  the  constants 
determined  in  Reference  2.  Thus  by  determining  Tr,  P^  and  0^  the  <J>  C.  PjC,  and  Pi  jC 
constants  can  be  calculated.  When  calculating  these  constants,  the  <|>  C  constant  should  be 
found  first  from  Equation  (5).  This  value  should  then  be  used  in  Equation  (6)  to  calculate  the 
P|C  constant.  Finally  the  values  of  the  <j>  C  and  P^C  constants  should  be  used  in  Equation 
(7)  to  solve  for  the  Pj  constant. 
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The  following  least  squares  fit  of  the  experimental  phase  diagram  was  used  to 
determine  the  compositional  dependence  of  Tr. 

For  x£  0.15:  T^  =  a  +  bx  +  cx^ 

where  a  =  -  105.07.  b  =  2812.6,  c  =  -  7665.9 

For  x£  0.15:  T^  =  a  +  bx  +  cx^  +  dx^+ex4 

where  a  =  5.5465.  b  =  2388.8,  c  =  -  13864..  d  =  32282..  e  =  -  32675.  (8) 

In  fitting  the  data  with  x  >  0. 15  an  additional  data  point  was  added  at  the  PZT  50/50 
composition  at  a  temperature  of  -  273  °C.  This  was  required  to  keep  the  high  to  low 
temperature  rhombohedral  boundary  from  crossing  the  morphotropic  boundary  at  low 
temperatures.  The  fit  of  the  experimental  data  using  these  equations  is  shown  in  Reference  4 
in  Figure  2. 

Glazer  et  al.3  determined  the  spontaneous  tilt  angle  for  the  PZT  90/10  composition  at 
25  and  60  °C.  The  values  of  the  tilt  angle  reported  represent  the  rotation  about  the  cubic 
[111]  direction.  In  the  three-dimensional  energy  function  used  in  this  theory  components 
(O3)  of  the  resultant  tilt  were  included.  Thus  the  component  tilts  were  first  calculated  from 
the  resultant  tilts  that  were  reported  in  the  literature.  The  following  relation  was  derived  to 
relate  the  component  tilts  (83)  about  each  orthogonal  axis  to  the  resultant  tilt  (85)  about  the 
body  diagonal: 

cos  83  =  1/3  (2  +  cos  85)  (9) 

Another  relation  was  also  found  to  give  similar  results  as  Equation  (9)  for  small  angles: 

3l/3  83  =  8s  (10) 

Equation  (10)  was  used  to  calculate  the  component  tilt  (83)  from  Glazer  et  al.’s  data. 

This  data  was  then  used  to  determine  the  P[_  and  8^  constants  for  the  PZT  90/10 
composition  that  gave  the  best  least-squares  fit  using  Equation  (4)  combined  with  Equations 
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(5)  -  (8).  The  resulting  values  of  Pl  and  8j_  are  listed  in  Table  I.  The  experimental  data  and 
the  fit  using  these  constants  is  shown  in  Figure  5  in  Reference  4. 

The  rotostrictive  R44  constant  was  used  to  fit  the  low-temperature  rhombohedral 
spontaneous  strain  X4  data  (from  Reference  6)  for  the  PZT  90/10  composition  using 
Equation  (24)  from  Reference  1.  and  the  values  of  P^  and  determined  above.  Using  the 
value  of  the  electrostrictive  constant  Q44  of  0.049  m^/C^  from  Reference  7,  resulted  in  a 
poor  fit  of  the  experimental  data.  However,  a  good  fit  was  obtained  by  varying  both  the  Q44 
and  R44  constants.  The  best  least-squares  fit  resulted  in  Q44  and  R44  values  of 
0.0433  rnVc^  and  -  1.93  x  10“*  deg'^,  respectively.  The  experimental  strain  data  and  the  fit 
using  these  values  is  shown  in  Figure  1(a). 

The  change  in  the  Q44  constant  from  high  to  low  temperature  rhombohedral  phases 
may  indicate  that  the  tilting  of  the  oxygen  octahedra  influences  the  electrostriction,  resulting 
in  different  electrostrictive  constants  in  these  two  phases.  Another  possible  cause  of  this 
could  be  that  the  tilt  angle  coefficients  are  temperature  dependent  and/or  that  a  higher  order  tilt 
angle  coefficient  is  needed  (0^  term).  However,  to  not  complicate  matters  any  more  than 
they  already  are,  the  Q44  and  R44  constants  of  the  low-temperature  rhombohedral  state  were 
assumed  to  have  the  values  given  above  independent  of  temperature  and  composition.  The 
values  of  the  electrostrictive  constants  of  the  high  temperature  rhombohedral  state  were  still 
kept  the  same  as  given  in  Reference  7. 

Using  the  same  values  of  the  Q44  and  R44  constants  given  above  for  the  PZT  90/10 
composition,  values  of  P^  and  0l  were  found  that  gave  the  best  least-squares  fit  of  the 
experimental  spontaneous  strain  data  for  the  PZT  80/20  composition.  These  values  are  listed 
in  Table  II.  The  experimental  data  and  fit  using  these  values  are  plotted  in  Figure  1(b). 

The  spontaneous  strain  X4  was  measured  at  two  temperatures  in  the  low-temperarure 

rhombohedral  state  for  the  PZT  70/30  composition.  Since  there  were  only  two  data  points,  a 

procedure  was  used  to  combine  these  data  with  a  single  tilt  angle  value  at  9K  that  was 

measured  by  Amin  et  al.®  for  the  PZT  60/40  composition.  ! 
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The  first  step  was  to  make  an  initial  guess  at  a  value  of  P[_  for  the  PZT  70/30 
composition.  This  value  of  Pl  was  used  with  the  Pl  values  listed  in  Table  I  for  the  PZT 
90/10  and  80/20  compositions  to  solve  for  a  quadratic  compositional  dependence  of  Pl 
(three  unknown  constants  and  three  data  points).  This  equation  was  used  to  extrapolate  a 
value  of  Pl  to  the  PZT  60/40  composition,  which  allowed  a  value  of  ©l  t0  b*  found  that  fit 
the  tilt  angle  data  at  9K.  This  value  of  &l  was  used  with  the  0l  values  listed  in  Table  I  for 
the  PZT  90/10  and  80/20  compositions  to  solve  for  a  quadratic  compositional  dependence  of 
0L-  This  equation  was  then  used  to  calculate  a  value  of  6l  back  at  the  PZT  70/30 
composition.  This  value  was  used  with  the  initial  guess  of  Pl  to  determine  the  spontaneous 
strain  X4,  and  to  calculate  the  least-squares  error  with  the  two  data  points.  Another  guess  of 
Pl  for  the  PZT  70/30  composition  was  then  made,  and  the  entire  procedure  repeated  until  the 
values  that  gave  the  best  least-squares  error  were  found. 

A  problem  developed  from  the  procedure  described  above:  Pl  became  less  than  Pj| 
for  compositions  near  PZT  60/40.  The  difference  in  the  values  between  Pl  and  P|^ 
decreased  from  around  PZT  80/20  to  PZT  60/40,  which  seems  to  be  in  agreement  with 
pyroelectric  measurements.^  However,  it  was  felt  that  Pl  should  probably  not  become  less 
than  P(|  (the  polarization  should  increase  when  going  to  the  low -temperature  rhombohedral 
phase).  Thus  an  additional  manipulation  was  used  to  keep  Pl  from  becoming  less  than  Pf^ 
To  allow  the  difference  between  PLand  Pj$  to  decrease  as  the  titanium  content  increased,  but 
to  keep  this  difference  from  going  to  zero  within  the  stability  region  of  the  low-temperature 
rhombohedral  phase,  an  additional  Pl  data  point  was  included  in  determining  the 
compositional  dependence  of  Pl-  This  value  was  arbitrarily  chosen  by  letting  Pl  become 
equal  to  Pj^  at  the  PZT  50/50  composition.  This  effectively  moved  the  composition  where 
Pl  became  equal  to  Pj_j  into  a  metastable  region. 

This  value  of  Pl  at  PZT  50/50  (listed  in  Table  I)  was  then  used  with  the  previous 
values  of  Pl  for  the  PZT  90/10,  80/20,  and  70/30  compositions  to  determine  the 
compositional  dependence  of  Pl  by  fitting  the  data  with  a  cubic  polynomial  equation  (four 
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unknowns  and  four  data  points).  This  equation  was  used  to  calculate  a  P|_  value  back  at  the 
PZT  60/40  composition.  Using  this  value  a  new  9^  value  for  PZT  60/40  was  found  by 
refitting  the  experimental  tilt  angle  data  point  at  9  K.  This  value  of  9l  for  PZT  60/40  was 
then  used  with  the  previous  values  of  9^  for  the  PZT  90/10,  80/20,  and  70/30  compositions 
to  determine  the  compositional  dependence  of  9^  from  a  cubic  polynomial  fit. 

The  final  values  of  Pj_  and  9l  used  in  the  fitting  are  listed  in  Table  I,  and  plotted  in 
Figure  2.  The  cubic  compositional  fits,  also  plotted  in  Figure  2,  of  these  data  are: 

P^  and  9L  =  a  +  bx  +  cx^  +  dx^, 


where  for  Pi 


a  =  0.28079 
b  =  -0.030117 
c  =  -0.46150 
d  =  1.8367 


and  for  9^:  a  =  0.91540 

b  =  12.967 
c  =  -40.255 
d  =  11.417 


The  experimental  spontaneous  strain  data  and  final  fit  for  the  PZT  70/30  composition 
are  plotted  in  Figures  1(c).  The  experimental  tilt  angle  data  point  and  final  fit  for  the  PZT 
60/40  composition  is  plotted  in  Figure  5  in  Reference  4. 

With  the  compositional  dependences  of  Tr,  Pl>  and  9^  determined  from  Equations 
(8),  (11),  and  (12),  the  tilt  angle  related  constants  PjC,  PijC,  and  4>  C  can  be  calculated 
from  Equations  (5)  -  (7). 


HI.  SUMMARY 

The  tilt  angle  related  coefficients,  (3|C,  (3 j  jC,  and  <j>  C,  were  related  to  three  new 
constants,  Tr,  Pl  and  0^,  at  the  transition  between  the  high  and  low  temperature 


rhombohedral  phases.  The  transition  temperature  between  these  phases.  Tr,  was 
determined  from  the  phase  diagram.  Values  of  the  spontaneous  polarization  and  tilt  angle  of 
the  low-temperature  rhombohedral  phase  at  Tr,  and  9^  were  determined  from 
experimental  spontaneous  strain  and  tilt  angle  data.  These  values  were  then  fit  with 
equations  to  determine  the  compositional  dependence  of  Pj_  and  9^  These  equations  will  be 
used  in  the  last  paper  in  this  series  to  calculate  the  values  of  (5^,  j,  and  (p  for  theoretical 
calculations  of  the  phase  stability  and  properties  of  the  low-temperature  rhombohedral  phase. 
As  a  summary  of  this  paper  Table  II  gives  values  of  these  constants  at  four  PZT 
compositions. 
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UST  OF  FIGURES  AND  TABLES 

Figure  I.  The  spontaneous  strain  *4  in  the  low-temperature  rhombohedral  state  plotted 
versus  temperature  for  the  PZT  90/10,  80/20,  and  70/30  compositions.  The 
data  points  are  from  Reference  6.  The  solid  curves  are  the  theoretical  fits  of  the 
data. 

Figure  2.  The  spontaneous  polarization  (PiJ  and  tilt  angle  (&l)  of  the  low-temperature 

rhombohedral  state  at  T^  plotted  versus  composition.  The  data  points  are  from 
Table  I.  The  solid  curves  are  fits  of  the  data  using  Equations  (11)  and  (12). 

Table  I.  Values  of  P^  and 

Table  II.  Values  of  the  Tilt  Angle  Related  Coefficients  at  25  °C. 
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Table  I.  Values  of  Pl  and  0l 


Zr/Ti 

PL  (C/m2) 

0L  (deg.) 

90/10 

0.275 

1.82 

80/20 

0.271 

1.99 

70/30 

0.280 

1.49 

60/40 

50/50 

0.380 

0.392 

Table  II.  Values  of  the  Tilt  Angle  Related  Coefficients  at  25  °C. 


Mole  Fraction  PbTiC^  in  PZT 


0.1 

0.2 

0.3 

0.4 

(105  J/m3  Deg2) 

5.618 

3.395 

4.343 

29.79 

^(lO4  J/m3  Deg4) 

2.506 

1.674 

1.682 

4.174 

(jKlO^m/FDeg2) 

-  9.626 

-  6.428 

-  6.502 

-31.82 
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THERMODYNAMIC  THEORY  OF  THE  LEAD  ZIRCONATE-TITANATE  SOLID 
SOLUTION  SYSTEM,  PART  V:  THEORETICAL  CALCULATIONS 

MJ.  Haun,t  E.  Furman, S J.  Jang,  and  L.E.  Cross 


Materials  Research  Laboratory 
The  Pennsylvania  State  University 
University  Park,  PA  16802 

In  this  final  paper  of  a  sequence  of  five  papers  presenting  details  of  a  thermodynamic 
phenomenology  for  the  whole  PbZr03:PbTi03  solid  solution  family,  the  constants  derived  in  the 
earlier  papers  are  used  to  calculate  the  temperature  and  composition  dependence  of  the  Elastic 
Gibbs  Free  Energy  AGi  at  zero  stress,  the  spontaneous  electric  polarization  and  the  spontaneous 
tilt  of  the  oxygen  octahedra  in  a  single  domain,  the  anisotropic  weak  field  dielectric  permittivity  and 
the  piezoelectric  constants  for  compositions  all  across  the  phase  diagram. 

These  data  can  be  used  to  determine  the  intrinsic  component  of  the  behavior  of  any  pure 
lead  zirconate:lead  titanate  in  any  part  of  the  phase  field  and  at  any  temperature,  and  thus  can  form  a 
first  step  in  the  separation  of  the  experimentally  observed  properties  into  their  intrinsic  (single 
domain)  component  and  their  extrinsic  domain  and  defect  related  responses. 

The  function  proposed  is  certainly  not  the  final  word  for  PZT,  and  will  clearly  be  subjected 
to  continuous  modification  and  improvement  as  new  and  better  experimental  data  becomes 
available.  The  purpose  of  the  exercise  has  been  to  correlate  a  very  wide  panoply  of  experimental 
data  within  a  single  function  and  to  extract  from  it  the  constants  which  would  be  directly  accessible 
if  single  domain  single  crystals  could  be  grown  at  the  interesting  PbZr03:PbTiC>3  compositions. 
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I.  INTRODUCTION 


This  paper  is  the  final  paper  in  a  senes  of  five  papers^  descnbing  the  development 
of  a  thermodynamic  theory  for  the  entire  lead  zirconate-titanate  (PZT)  solid  solution  system. 
In  this  paper  theoretical  calculations  will  be  made  using  the  equations  that  were  denved  in  the 
first  paper^  of  the  series,  and  the  values  of  the  coefficients  determined  in  the  second,  third, 
and  fourth  papers.^  The  antiferroelectric  calculations  will  be  made  using  the  coefficients 
that  were  determined  in  Reference  5. 

The  free  energies  of  each  phase  will  be  calculated  in  the  next  section  and  used  to 
determine  the  phase  diagram.  In  Section  HI  the  spontaneous  polarization  and  tilt  angle  will 
be  calculated.  The  dielectric  and  piezoelectric  properties  will  then  be  calculated  in  Sections 
IV  and  V,  followed  by  a  summary  and  discussion  of  the  future  applications  of  this  theory. 


0.  PHASE  DIAGRAM 

In  the  second  paper  of  the  series  products  of  the  dielectric  stiffness  coefficients  (cqj) 
and  the  Curie  constant  (C)  were  shown  to  be  independent  of  C.  These  new  constants  (cijj  C) 
can  be  used  to  calculate  products  of  the  free  energies  and  the  Curie  constant  (AG  C),  which 
are  also  independent  of  C.  Since  the  Curie  constant  is  the  same  in  all  of  the  solutions  of  the 
energy  function,  the  AG  C  products  can  be  used  to  determine  the  phase  stability  independent 
of  the  effect  of  the  compositional  dependence  of  C. 

In  Figure  1  the  product  of  the  energies  and  Curie  constant  (C)  are  plotted  versus 
composition  at  three  temperatures.  The  solution  with  the  lowest  AG  C  product  corresponds 
to  the  stable  phase,  and  a  phase  transition  occurs  when  the  curves  cross.  At  25  °C 
[Figure  1(c)]  the  tetragonal  solution  has  the  lowest  AG  C  on  the  lead  titanate  side  of  the 
figure,  and  is  thus  the  stable  phase.  Moving  across  the  figure  towards  lead  zirconate  the 
tetragonal  phase  remains  stable  to  the  dashed  line  labeled  5.  At  this  point  the  AG  C  of  the 
high-temperature  rhomboheilral  solution  becomes  less  than  that  of  the  tetragonal  phase,  and  a 
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phase  transition  occurs  between  these  phases  (the  morphotropic  phase  boundary). 

Continuing  across  the  figure  the  high-temperature  rhombohedral  phase  remains 
stable  to  the  dashed  line  labeled  4,  where  a  transition  to  the  low-temperature  rhombohedral 
phase  occurs.  This  phase  is  then  stable  until  the  dashed  line  labeled  2  is  reached,  and  a 
transition  to  the  anuferroelectnc  orthorhombic  phase  occurs.  At  higher  temperatures  an 
additional  region  of  high-temperature  rhombohedral  phase  develops  between  the 
low-temperature  rhombohedral  and  antiferroelectric  orthorhombic  phase  fields.  The 
low-temperature  rhombohedral  phase  field  also  becomes  narrower  as  the  temperature  is 
increased. 

The  AG  C  curves  for  the  low-temperature  rhombohedral  and  antiferroelectric 
orthorhombic  phases  shown  in  Figure  1  were  only  drawn  over  the  regions  where  these 
phases  were  stable.  The  compositional  dependence  of  the  antiferroelectric  orthorhombic 
phase  was  determined  by  equating  the  energies  of  this  phase  with  the  high-temperature 
rhombohedral  phase  at  the  transition  between  these  states  (described  in  Reference  5).  Since 
this  transition  only  extends  over  a  narrow  compositional  region,  the  above  method  will  only 
apply  over  this  region.  Thus  the  energy  of  the  antiferroelectric  orthorhombic  phase  should 
only  be  calculated  within  this  region. 

In  the  case  of  the  low-temperature  rhombohedral  phase  the  AG  C  can  be  calculated 
over  a  very  small  region  into  the  adjacent  phase  fields,  at  which  point  either  the  spontaneous 
polarization  becomes  imaginary,  or  the  low-temperature  rhombohedral  phase  undergoes  a 
metastable  phase  transition  to  the  cubic  state.  At  this  point  additional  work  needs  to  be 
conducted  to  more  completely  understand  what  is  actually  happening.  This  is  particularly 
important  in  gaining  a  better  understanding  of  the  change  in  the  spontaneous  polarization 
from  the  high  to  low  temperature  rhombohedral  phases,  as  will  be  discussed  in  more  detail  in 
the  next  section. 

Calculating  the  AG  C  product  of  the  ferroelectric  orthorhombic  phase  provided  an 
independent  check  on  the  values  of  the  coefficients  of  the  energy  function.  As  shown  in 
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Figure  l  this  phase  was  found  to  be  metastable  across  the  phase  diagram,  indicating  that  the 
ratios  of  the  coefficients  are  at  least  in  the  right  range. 

By  plotting  the  composition  versus  temperature  where  the  AG  C  products  of  the 
solutions  cross,  the  phase  diagram  is  obtained  as  shown  in  Figure  2.  The  data  points  are 
from  the  experimental  phase  diagram.^  and  the  solid  curves  are  from  the  theoretical 
calculations.  The  excellent  agreement  between  the  experimental  data  and  theoretical 
calculations  is  to  be  expected,  because  the  data  was  used  to  determine  values  of  the 
coefficients  in  such  a  way  as  to  cause  the  energies  of  the  adjacent  phases  of  a  transition  to  be 
equal  at  the  transition.  However,  this  agreement  shows  that  the  energy  function  and  values 
of  the  coefficients  that  have  been  determined  will  quantitatively  model  all  of  the  known  phase 
transitions  of  the  PZT  solid  solution  system. 

The  inflection  point  in  the  transition  temperature  Tq  between  the  cubic  and 
high-temperature  rhombohedral  phases  near  the  PZT  90/10  composition  (see  Figure  2) 
appears  to  be  due  to  the  tricritical  behavior.  The  Curie -Weiss  temperature  Tq  should  form  a 
continuous  function  with  a  continuous  slope  across  the  phase  diagram,  because  it  is 
determined  from  the  cubic  phase.  In  the  second-order  region  Tq  is  equal  to  Tq,  but  in  the 
first-order  regions  Tq  becomes  less  than  Tq.  The  Tq  -  Tq  difference  increases  very  quickly 
over  a  narrow  compositional  region,  from  zero  at  the  tricritical  point  to  42  °C  for  lead 
zirconate  (see  Figure  9  in  Reference  2).  This  therefore  provides  an  explanation  for  the 
inflection  point  in  Tq  near  the  PZT  90/10  composition. 

A  question  now  arises  as  to  why  a  similar  change  in  Tq  is  not  seen  at  the  second 
tricritical  point  near  the  PZT  28/72  composition.  This  may  be  because  the  Tq  -  Tq  difference 
only  changes  from  zero  at  the  tricritical  point  to  13.4  °C  for  lead  titanate,  and  this  smaller 
change  occurs  over  a  much  wider  compositional  range.  In  addition  there  is  not  enough  Tq 
data  on  the  phase  diagram  in  this  region  to  really  be  able  to  notice  this  type  of  change  in  Tq 
Additional  high-temperarure  x-ray  diffraction  work  is  needed  to  more  precisely  locate  the 
second  tricritical  point.  In  the  process  of  this  work  it  would  be  interesting  to  try  and  see  if 
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T(;  does  change  in  a  similar  (but  smaller)  way  as  occurs  on  the  lead  zirconate  side. 


m.  SPONTANEOUS  POLARIZATION  AND  TILT  ANGLE 

Figure  3  shows  the  calculated  spontaneous  polarization  plotted  versus  temperature 
for  five  tetragonal  PZT  compositions.  The  first  order  transitions  of  the  PbTi03  and  PZT 
20/80  compositions  can  be  seen  by  the  discontinuous  changes  in  the  polarization  at  Tq  The 
other  three  compositions  have  second  order  transitions  causing  the  polarization  to  change 
continuously  at  T^.  The  dielectric  stiffness  coefficient  aj  j  changes  sign  at  a  uncritical  point 
between  the  PZT  20/80  and  30/70  compositions  causing  the  transition  to  change  from  first  to 
second  order.  This  tricritical  point  was  estimated  to  occur  at  the  PZT  28/72  composition 
from  an  extrapolation  in  Reference  2. 

The  calculated  spontaneous  polarization  data  shown  in  Figure  3  are  in  good 
agreement  with  the  available  experimental  data.  The  value  of  0.75  C/m2  at  25  °C  for  lead 
titanate  is  equal  to  Gavrilyachenko  et  ais7  experimental  single-crystal  measurement.  The 
polarization  values  for  the  PZT  40/60  composition  are  slightly  larger  than  the  single-crystal 
measurements  from  25  to  200  °C  of  Tsuzuki  et  al.^  They  commented  that  the  polarization 
was  not  completely  saturated  in  their  hysteresis  loops,  which  would  result  in  lower  values  of 
the  polarization.  The  close  agreement  between  the  calculated  polarization  and  experimental 
measurements  indicates  that  the  proper  values  of  the  electrostrictive  constants  were. 

The  spontaneous  polarization  of  four  rhombohedral  compositions  with  the  high  to 
low  temperature  rhombohedral  phase  transition  is  plotted  versus  temperature  in  Figure  4. 

The  resultant  polarization  (=  3  ^2  Pj)  along  the  [l  1 1]  direction  was  plotted  in  this  figure. 
The  small  first  order  change  in  the  polarization  of  the  PZT  90/10  composition  at  the 
high-temperature  rhombohedral  to  cubic  transition  (T^)  can  be  seen  in  Figure  4  (a).  The 
other  three  compositions  [Figures  4  (b)  -  (d)J  all  have  a  second  order  transition  at  T^. 

This  is  because  the  fourth  order  rhombohedral  dielectric  stiffness  (a^  i  +  a^)  is  slightly 
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negative  for  the  PZT  90/10  composition,  and  then  changes  sign  at  a  tricriticai  point  (close  to 
PZT  90/10)  and  becomes  positive  for  the  other  three  compositions. 

The  high  to  low  temperature  rhombohedral  phase  transition  is  first  order  for  all 
compositions,  as  can  be  seen  in  Figure  4  by  the  discontinuous  change  in  the  polarization  at 
this  transition.  This  discontinuous  change  in  the  polarization  increases  from  PZT  90/10  to 
80/20,  but  then  decreases  as  the  titanium  content  is  increased  further,  and  becomes  very 
small  for  the  PZT  60/40  composition.  This  may  explain  why  Zhuang  et  al.^  found  from 
pyroelectric  measurements  that  the  polarization  only  slightly  changes  through  the  high  to  low 
temperature  rhombohedral  transition  for  the  PZT  70/30  and  60/40  compositions,  but  a  large 
change  occurs  for  the  PZT  90/10  and  94/6  compositions. 

The  spontaneous  polarization  data  for  the  PZT  90/10  composition  [Figure  4  (a)]  is  in 
fairly  good  agreement  with  the  single-crystal  measurements  by  Clarke  and  Glazer.  ^  Some 
of  their  experimental  data  in  the  high-temperature  rhombohedral  region  was  used  to 
determine  the  electrostrictive  Q44  constant  for  this  composition,  and  thus  good  agreement 
occurs  in  this  region.  However,  their  data  showed  a  decrease  in  the  polarization  on  cooling 
in  the  high-temperature  rhombohedral  state  close  to  the  transition  to  the  low -temperature 
phase,  possibly  due  to  a  leakage  of  charge  during  their  measurements.  In  this  region  the 
theoretically  calculated  polarization  continues  to  increase  as  the  temperature  is  decreased. 

The  spontaneous  tilt  angle  was  also  calculated  versus  temperature  for  several 
compositions  as  shown  in  Figure  5.  The  resultant  tilt  angle  about  the  cubic  [111]  direction 
(=  3^  83)  was  plotted.  The  data  points  are  the  only  available  experimental  values  of  the  tilt 
angle,  and  were  determined  from  neutron  diffraction  by  Glazer  et  al.  ^  and  Amin  et  al.  ^ 
These  data  were  used  in  determining  values  for  two  of  the  coefficients  (see  Reference  4),  and 
thus  the  theory  and  experimental  data  should  agree.  The  theoretical  calculations  can  be  used 
to  determine  values  of  the  tilt  angle  at  other  compositions  and  temperatures,  as 
shown  in  this  figure. 
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The  spontaneous  tilt  angle  was  also  plotted  in  Figure  6  versus  composition  at  several 
different  temperatures.  The  curve  labeled  represents  the  tilt  angle  at  the  high  to  low 
temperature  ihombohedral  transition  (Tr).  The  discontinuity  in  the  slope  of  the  curves  at  the 
PZT  85/15  composition  occurs  due  to  the  relations  that  was  used  to  fit  Tr  [Equation  (8)  in 
Reference  4]. 

Figures  4  and  5  show  that  the  spontaneous  polarization  and  tilt  angle  increase  very 
quickly  after  transforming  to  the  low-temperature  rhombohedral  phase,  especially  the 
compositions  with  greater  titanium  content.  This  may  indicate  a  limitation  of  the  present 
theory.  For  the  PZT  70/30  and  60/40  compositions  the  discontinuous  change  in  the 
polarization  at  the  high  to  low  temperature  rhombohedral  transition  is  very  small  compared  to 
the  sudden  increase  in  polarization  that  occurs  just  below  the  transition.  This  increase  in 
polarization  in  the  low-temperature  rhombohedral  phase  is  even  more  noticeable  when  the 
polarization  is  ploned  versus  composition  at  different  temperatures  as  shown  in  Figure  7. 

If  the  spontaneous  polarization  as  a  function  composition  does  actually  increase  in  the 
low-temperature  rhombohedral  phase  as  much  as  is  shown  in  this  figure,  then  it  seems  that 
the  previous  experimental  measurements  of  the  polarization  on  ceramic  samples  would  have 
found  a  similar  effect. 

The  polarization  should  increase  going  to  the  low -temperature  rhombohedral  phase, 
but  the  increase  calculated  from  this  theory  seems  to  be  too  large.  Two  possible  explanations 
of  this  large  increase  in  polarization  in  the  low-temperature  rhombohedral  phase  will  now  be 
described. 

The  first  possibility  is  that  this  large  increase  may  actually  occur,  and  that  the 
previous  measurements  on  ceramic  samples  missed  this  effect,  because  not  enough 
measurements  were  made.  The  effect  of  the  morphotropic  boundary  on  the  ceramic  data, 
where  a  coexistence  of  tetragonal  and  rhombohedral  phases  occurs,  may  have  also  made  it 
more  difficult  to  measure  the  polarization  in  the  narrow  region  where  the  high-temperature 
rhombohedral  phase  is  stable  Additional  polarization  measurements  on  homogeneous 
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ceramic  samples  would  be  useful  in  trying  to  see  how  the  polarization  changes  in  this  region. 

The  second  possibility,  which  seems  more  likely,  is  that  the  present  theory  does  not 
adequately  model  the  spontaneous  polarization  in  the  low-temperature  rhombohedral  phase 
for  the  compositions  with  greater  titanium  content.  This  could  be  due  to  the  assumptions 
made  in  determining  the  coefficients.  For  example,  the  Q44  and  R44  constants  were 
assumed  to  be  independent  of  composition  and  temperature  in  the  low-temperature 
rhombohedral  phase.  The  values  of  these  constants  were  determined  at  the  PZT  90/10 
composition,  and  then  used  for  the  other  compositions.  The  calculated  polarization  in  the 
low-temperature  rhombohedral  phase  of  the  PZT  90/10  composition  does  not  have  the  large 
increase  that  occurs  in  the  PZT  70/30  and  60/40  compositions.  Therefore,  it  seems  likely 
that  the  Q44  and  R44  constants  are  dependent  on  composition,  and  if  this  was  accounted  for 
the  large  increase  in  polarization  in  the  PZT  70/30  and  60/40  compositions  would  be 
decreased.  However,  at  this  time  there  is  not  enough  experimental  data  available  in  the 
literature  to  determine  these  compositional  dependences  Additional  tilt  angle  data  obtained 
from  neutron  diffraction  would  possibly  solve  this  problem. 

Another  significant  change  could  be  made  to  the  present  theory  by  adding  a 
temperature  dependence  to  the  second-order  tilt  angle  related  coefficient  (p  j ).  and/or  a 
sixth-order  tilt  angle  term.  Again  the  problem  with  this  is  that  there  is  not  enough  data  to 
properly  determine  all  of  these  coefficients.  The  addition  of  a  sixth-order  tilt  angle  term  will 
also  greatly  complicate  the  theory,  because  the  equations  for  the  spontaneous  polarization  and 
lilt  angle  in  the  low-temperature  rhombohedral  phase  will  change  form  quadratic  to  quamc. 

The  above  discussion  has  shown  that  there  may  be  some  limitation  to  how  well  the 
present  theory  can  quantitatively  calculate  the  spontaneous  polarization  of  the 
low-temperature  rhombohedral  phase  for  the  compositions  with  greater  titanium  content. 

Due  to  the  lack  of  tilt  angle  data  and  to  avoid  complicating  the  theory  significantly  more  than 
it  already  is,  the  present  theory  is  probably  the  best  that  it  can  be  ai  this  time,  and  is  probably 
adequate  for  most  applications. 
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IV.  DIELECTRIC  PROPERTIES 

The  calculated  dielectric  susceptibility  coefficients  of  PZT  50/50  were  plotted  versus 
temperature  in  Figure  8.  In  comparison  to  the  properties  of  PbTi03  (see  Figure  8  in 
Reference  13).  the  dielectric  susceptibilities  and  anisotropy  have  become  much  larger.  The 
transition  at  T^  is  also  second  order  for  PZT  50/50,  while  PbTiO^  is  first  order.  The 
second-order  transition  causes  the  dielectric  susceptibility  to  go  to  infinity  at  Tq 

The  compositional  dependence  of  the  dielectric  susceptibility  coefficients  at  25  ~C  is 
plotted  in  Figure  9.  A  peak  occurs  in  these  coefficients  at  the  PZT  50/50  composition  due  to 
the  peak  that  formed  in  the  Curie  constant  (see  Figure  2  in  Reference  3).  The  increase  in 
dielectric  anisotropy  (r|  ^  1/TI33)  that  occurs  from  PbTiC^  to  the  PZT  50/50  composition  is 
due  to  the  change  in  the  ratios  of  the  dielectric  stiffness  coefficients  (0112/04  [.  04  I2^al  1 1- 
and  04  23/04 1 1 )  35  a  function  of  composition. 

On  the  rhombohedral  side  of  the  peak  the  anisotropy  decreases  and  T|  n’  eventually 
even  becomes  less  than  433'.  This  effect  is  at  least  partially  due  to  compositional 
dependence  of  the  ratios  of  the  dielectric  stiffness  coefficients,  but  may  also  be  due  to  an 
internal  pressure  that  may  develop  within  the  perovskite  structure  when  a  larger  zirconium 
ion  is  replaced  by  the  smaller  titanium  ion.  Amin  et  al.  ^  showed  that  for  the  PZT  50/50 
composition  the  application  of  hydrostatic  stress  would  cause  the  anisotropy  to  decrease,  and 
thus  this  might  partially  account  for  the  small  dielectric  anisotropy  on  the  PbZrC>3  side  on  the 
PZT  system.  Using  the  present  theory  of  PZT  the  effect  of  mechanical  boundary  conditions 
on  the  properties  and  phase  stability  should  now  be  investigated  in  more  detail  across  the 
entire  PZT  system. 

The  dielectric  susceptibility  coefficients  of  two  rhombohedral  compositions.  PZT 
60/40  and  90/10.  were  plotted  versus  temperature  in  Figures  10  and  11.  A  similar  change  in 
the  behavior  of  the  dielectric  anisotropy  occurs  from  PZT  60/40  to  90/10  as  occurred  from 
PZT  50/50  to  PbTiC>3  For  PZT  90/10  and  PbTiC^  (see  Figure  8  in  Reference  13)  the 
dielectric  susceptibility  parallel  to  the  polar  axis  (433  or  1133  )  is  larger  than  the  perpendicular 


susceptibility  CH 1 L  or  H 1 1 )  at  hi8h  temperatures,  but  smaller  at  low  temperatures. 

At  high  temperatures  close  to  for  rhombohedral  compositions  close  to  the 
morphotropic  boundary  (Ti  content  >  0.45)  a  problem  develops  when  calculating  the 
dielectric  susceptibility  coefficient  r|  j  ^  using  the  coefficients  from  References  2  and  3.  At  a 
temperature  below  T^  the  r)i  i  coefficient  will  increase  to  infinity,  and  then  become  negative 
between  this  temperature  and  T^.  This  is  simply  due  to  the  values  of  the  coefficients  used, 
and  indicates  that  in  this  region  these  values  are  not  quite  right.  Because  of  this  problem, 
care  should  be  taken  in  calculating  the  rhombohedral  dielectric  susceptibilities  at  high 
temperatures  for  compositions  near  the  morphotropic  boundary.  Considering  that  all  of  the 
coefficients  were  assumed  to  be  independent  of  temperature,  except  aj,  and  that  some  of  the 
coefficients  were  determined  at  (232  -  492  °C)  and  others  at  4.2  K,  it  is  amazing  that 
everything  works  as  well  as  it  does,  and  that  more  problems  have  not  occurred. 

Away  from  this  region  the  calculated  dielectric  properties  are  in  fairly  good 
agreement  with  the  available  experimental  data.  Haun  et  ai.^  showed  that  the  calculated 
dielectric  susceptibility  coefficients  and  anisotropy  of  PbTiOj  were  in  good  agreement  with 
Fesenko  et  al.'s^  experimental  single-crystal  data  at  low  temperatures,  where  the  defect 
contributions  to  the  experimental  measurements  had  "frozen  out".  The  theoretically 
calculated  dielectric  properties  near  the  morphotropic  boundary  have  similar  values  as  the 
experimental  measurements  on  single  crystals  by  Tsuzuki  et  al.^  Good  agreement  was  also 
found  between  the  calculated  dielectric  properties  and  Clarke  and  Whatmore's^  PZT  90/10 
single-crystal  measurements. 


V.  PIEZOELECTRIC  PROPERTIES 

The  calculated  piezoelectric  voltage  coefficients  gy  (assumed  to  be  equal  to  the  by 
coefficients)  for  the  PZT  50/50,  60/40,  and  90/10  compositions  were  plotted  versus 
temperature  in  Figures  ( 1 2 )  - 1 14).  The  gy  coefficients  plotted  in  these  figures  are  based  on 


the  cubic  axes.  Thus  the  X3  axis  is  along  the  polar  axis  for  the  tetragonal  coefficients,  but  not 
for  the  rhombohedral  coefficients. 

Since  the  electrostrictive  constants  were  assumed  to  be  independent  of  temperature, 
the  negative  temperature  dependence  of  the  gjj  coefficients  was  caused  by  the  temperature 
dependence  of  the  spontaneous  polarization.  The  PZT  90/10  composition  has  a  first-order 
change  in  the  spontaneous  polarization,  which  causes  discontinuities  in  the  gjj  coefficients  at 
T(3_  The  PZT  50/50  and  60/40  compositions  undergo  second  order  transitions,  and  thus  the 
gjj  coefficients  change  continuously  from  a  value  of  zero  at  T^. 

The  compositional  dependence  of  the  gjj  coefficients  at  25  °C  is  plotted  in  Figure  15 
Since  the  electrostrictive  constants  were  given  a  compositional  dependence  (see  Figure  2  in 
Reference  17),  the  compositional  dependence  of  the  gjj  coefficients  depends  on  the 
electrostrictive  constants  and  the  spontaneous  polarization  (see  Figure  7).  The  result  is  that 
the  gjj  coefficients  change  with  composition  as  shown  in  Figure  15.  The  compositional 
dependences  of  the  gjj  coefficients  are  much  smaller  than  that  of  the  djj  coefficients,  which 
will  be  presented  next. 

The  piezoelectric  charge  coefficients  djj  of  the  PZT  50/50  composition  were  plotted 
versus  temperature  in  Figure  16.  The  positive  temperature  dependence  of  the  djj  coefficients 
was  caused  by  the  strong  temperature  dependence  of  the  dielectric  susceptibility  coefficients. 
In  Figure  17  the  djj  coefficients  are  plotted  versus  composition  at  25  °C,  showing  the  large 
peak  in  these  coefficients  near  the  morphotropic  boundary.  The  large  increase  of  the  r)  1 1 
coefficient  from  PbTiOj  to  the  PZT  50/50  composition  (see  Figure  9)  caused  a  large  increase 
to  also  occur  in  the  d  [5  coefficient. 

For  the  rhombohedral  compositions  the  djj  coefficients  were  calculated  based  on  the 
original  cubic  axes^jesuiting  in  four  independent  nonzero  coefficients.  The  relations  for  the 
rhombohedral  coefficients  are  more  complicated  than  those  for  the  tetragonal  coefficients, 
because  two  terms  are  involved  in  three  of  the  four  relations  [see  Equations  (44)  and  (45)  in 
Reference  1  ].  These  two  terms  are  often  competing  with  one  term  being  positive,  while  the 
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ocher  is  negative.  In  addition  the  rj  1 2  coefficient  can  be  positive  or  negative  depending  on 
the  composition  and  temperature. 

The  dy  coefficients  are  plotted  versus  temperature  for  two  rhombohedral 
compositions,  PZT  60/40  and  90/10.  in  Figures  18  and  19.  The  djj  and  d  15  coefficients 
form  similar  shaped  curves  for  both  compositions,  but  the  d3  [  and  d^4  coefficients  have 
definite  differences. 

For  the  PZT  60/40  composition  the  dj4  coefficient  is  negative  at  all  temperatures, 
and  becomes  more  negative  with  increasing  temperature.  The  dj4  coefficient  of  the  PZT 
90/10  composition  is  slightly  negative  at  Low  temperatures,  and  then  becomes  positive  at 
about  - 100  °C,  and  continues  to  increase  with  increasing  temperature.  These  differences 
are  due  to  the  t\^2  coefficient,  which  is  negative  for  the  PZT  60/40  composition,  and  positive 
for  the  PZT  90/10  composition  (except  at  low  temperatures  where  it  changes  sign  and 
becomes  negative). 

The  d3i  coefficient  of  the  60/40  composition  is  negative,  and  becomes  more 
negative  with  increasing  temperature.  For  the  PZT  90/10  composition  the  d3  \  coefficient  is 
negative  at  low  temperatures,  becomes  more  negative  with  increasing  temperature,  and  then 
suddenly  increases  and  becomes  positive  just  below  Tq.  This  behavior  is  due  to  the 
coefficient,  and  also  due  to  the  two  terms  in  the  d3  j  relation,  which  are  of  opposite  sign  with 
one  dominating  at  low  temperatures  and  the  other  at  high  temperatures. 

The  rhombohedral  djj  coefficients  described  above  are  based  on  the  original  cubic 
axes.  At  this  point  what  needs  to  be  done  in  the  future  is  to  rotate  the  axes,  so  that  the  new 
X3  axis  is  along  the  rhombohedral  polar  direction.  This  would  reduce  the  number  of 
independent  nonzero  coefficients  to  three  (CI33,  d3[,  and  d^),  similar  to  the  tetragonal  djj 
matrix,  and  thus  simplify  the  analysis  of  these  coefficients. 

The  calculated  d^  coefficients  have  similar  values  as  the  experimental  ceramic  data 
near  the  morphotropic  boundary.  The  only  available  single-crystal  data  was  measured  by 
Gavrilyachenko  and  Fesenko 1  ^  on  lead  titanate.  Their  data  were  compared  with  the 


theoretical  lead  titanate  calculations  in  Reference  13.  Fairly  good  agreement  was  found, 
especially  for  the  dj  j  coefficient. 

VI.  SUMMARY  AND  FUTURE  APPLICATIONS  OF  THE  THEORY  OF  PZT 

The  applications  of  compositions  of  the  PZT  solid  solution  system  as  piezoelectric 
transducers,  pyroelectric  detectors,  electro-optic  devices,  and  explosively  induced  charge 
storage  devices  were  described  in  the  first  paper  of  this  series  to  demonstrate  the 
technological  importance  of  PZT.  These  applications  use  PZT  in  polycrystalline  ceramic 
form,  and  thus  the  properties  of  these  ceramics  are  well  established  in  the  literature. 

However,  the  mechanisms  contributing  to  these  outstanding  ceramic  propenies  of  PZT  are 
not  well  understood,  because  of  the  complexity  of  the  interactions  within  the  ceramic. 

A  first  step  in  the  analysis  of  a  ferroelectric  ceramic  material  is  to  separate  the 
intrinsic  and  extrinsic  contributions  to  the  properties.  The  intrinsic  contributions  result  from 
the  averaging  of  the  single-domain  single-crystal  properties,  while  the  extrinsic  contributions 
arise  from  the  interactions  at  grain  or  phase  boundaries  and  from  the  domain  wall  or  thermal 
defect  motions.  Unfortunately,  due  to  the  difficulty  of  growing  good  quality  single  crystals 
of  PZT,  very  little  single-crystal  data  is  available.  Thus  the  goal  of  this  project  was  to 
develop  a  thermodynamic  phenomenological  theory  to  calculate  the  single-domain  properties 
of  PZT.  This  theory  could  then  be  used  to  separate  the  intrinsic  and  extrinsic  contributions  to 
the  ceramic  properties.  In  addition  there  are  several  other  applications  of  this  theory,  which 
will  be  described  later  in  this  summary. 

In  the  first  paper  of  this  series  a  modified  elastic  Gibbs  free  energy  function  was 
expanded  in  powers  of  the  ferroelectric  and  antiferroelecrric  measures  of  the  polarizauon.  .An 
additional  order  parameter  was  also  included  to  account  for  the  tilting  of  the  oxygen 
octahedra  in  the  low-temperarure  rhombohedral  phase.  The  resulting  energy  function  can  be 
used  to  model  the  phase  transitions  and  single-domain  propenies  of  the  PZT  system. 


Solutions  to  this  energy  function  were  used  to  derive  relations  for  the  energies;  spontaneous 
polarizations,  strains,  and  tilt  angles;  and  dielectric  and  piezoelectric  properties  corresponding 
to  the  different  phases  in  the  PZT  system. 

Ail  of  the  coefficients  of  the  energy  function  were  assumed  to  be  independent  of 
temperature,  except  the  ferroelectric  and  antiferroelectric  dielectric  stiffness  coefficients  04 
and  <Ji  ,  which  was  given  a  linear  temperature  dependence  based  on  the  Curie-Weiss  law. 

The  experimental  phase  diagram  was  extensively  used  in  the  evaluation  of  the  coefficients, 
by  requiring  that  the  energies  of  the  adjacent  phases  of  a  phase  transition  to  be  equal  at  the 
transition.  The  first  partial  derivative  stability  conditions  were  also  used  as  additional 
relations  in  the  evaluation  procedure. 

Without  single-crystal  data,  the  development  of  this  theory  was  complicated  and 
involved  indirect  methods  of  determining  the  coefficients  of  the  energy  function.  Additional 
experimental  data  were  needed  to  determine  the  compositional  dependence  of  the  coefficients. 
A  sol-gel  procedure  was  used  to  prepare  pyre  homogeneous  PZT  powders  across  the  phase 
diagram  to  be  used  in  collecting  additional  PZT  data.  The  lanice  constants  versus 
temperature  were  determined  from  these  powders  from  high-temperature  x-ray  diffraction, 
and  used  to  calculate  the  spontaneous  strain.^’ 13,19  The  sol-gel  derived  powders  were  also 
used  to  fabricate  ceramic  samples  for  dielectric,  piezoelectric,  elastic,  pyroelectric,  and 
electrostrictive  measurements.-^ ’^’20-2 1.22  These  data  were  very  important  in  completing 
the  evaluation  of  the  coefficients  of  the  energy  function. 

The  compositional  dependence  of  the  cubic  electrostrictive  constants  across  the  PZT 
system  were  determined  from  a  combination  of  single  crystal  and  ceramic  data  using  senes 
and  parallel  averaging  relations  (analogous  to  the  Reuss  and  Voigt  methods  of  averaging  the 
elastic  constants).  ^  An  interesting  result  of  this  work  was  that  the  electrostrictive  constants 
form  a  peak  at  the  PZT  50/50  composition.  In  addition,  the  ratios  of  the  electrostrictive 
constants  change  with  composition  in  such  a  way  as  to  explain  the  compositional  dependence 
of  the  electromechanical  anisotropy  found  in  PZT  ceramics.^  The  electromechanical 
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anisotropy  is  very  large  in  compositions  near  the  end  members  PbTi03  and  PbZrC^,  but 
becomes  very  small  for  compositions  in  the  center  of  the  phase  diagram  near  the 
morpho tropic  phase  boundary. 

In  the  second  paper  of  this  series  the  electrostrictive  data  were  used  with  the 
spontaneous  strain  data  to  calculate  the  spontaneous  polarization,  which  was  used  to 
determine  values  of  the  higher-order  dielectric  stiffness  coefficients.  In  the  third  paper  of  this 
series  low-temperature  ceramic  dielectric  data  were  particularly  important  in  determining  the 
compositional  dependence  of  the  Curie  constant,  and  in  separating  the  sixth-order 
polarization  interaction  coefficients  (aj  12  ^d  0^23)-  The  tilt  angle  related  coefficients  were 
determined  in  the  fourth  paper  of  the  series  from  spontaneous  strain  and  tilt  angle  data.  The 
compositional  dependence  of  the  antiferroelectric  orthorhombic  free  energy  was  determined 
by  combining  data  from  a  separate  theory  developed  for  PbZrC>3  with  the  requirement  that 
the  energies  of  the  high-temperature  rhombohedral  and  antiferroelectric  orthorhombic  phases 
must  be  equal  at  the  transition  between  these  phases.5 

After  determining  values  of  the  coefficients  at  several  different  compositions  as 
described  above,  the  compositional  dependence  of  the  coefficients  was  determined  by  fining 
an  equation  to  these  values.  A  set  of  equations  was  then  established  that  can  be  used  to 
calculate  values  of  the  coefficients  at  any  composition. 

Using  these  equations  the  phase  stability,  spontaneous  polarization  and  tilt  angle, 
and  dielectric  and  piezoelectric  properties  were  calculated.  The  theoretically  calculated  phase 
diagram  was  shown  to  quantitatively  model  the  experimental  phase  diagram.  An  independent 
check  of  the  coefficients  was  also  made  by  calculating  the  free  energy  of  the  ferroelectric 
orthorhombic  stare,  which  showed  that  this  phase  was  metastable  across  the  phase  diagram 
as  experimentally  expected.  The  spontaneous  polarization  and  tilt  angle  were  calculated  for 
several  compositions  across  the  phase  diagram.  These  calculations  were  tn  good  agreement 
with  the  available  experimental  single-crystal  data.  The  dielectric  and  piezoelectric  propertu  s 
were  also  calculated  and  found  to  be  m  good  agreement  with  the  experimental  data.  As  a 
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summary  of  these  calculations  the  theoretical  properties  of  PZT  at  25°C  are  listed  in  Table  I. 

Some  Limitations  of  the  present  theory  were  found.  Due  to  the  lack  of  experimental 
data  for  the  low-temperature  rhombohedral  phase,  especially  tilt  angle  data,  some 
assumptions  were  necessary  in  determining  the  compositional  dependence  of  the  tilt  angle 
related  coefficients.  In  addition  these  coefficients  were  assumed  to  be  independent  of 
temperature,  and  only  tilt  angle  terms  up  to  the  fourth  power  were  included.  With  these 
assumptions  the  calculations  of  the  spontaneous  polarization,  tilt  angle,  and  dielectric 
properties  were  still  in  fairly  good  agreement  with  the  experimental  single-crystal  data  at  the 
PZT  90/10  composition. 

The  problems  developed  in  the  calculations  for  compositions  with  greater  titanium 
content.  The  spontaneous  polarization  calculated  in  the  low-temperature  rhombohedral  phase 
for  the  PZT  70/30  and  60/40  compositions  increased  more  than  would  be  expected.  Another 
problem  was  found  when  calculating  the  dielectric  susceptibility  (negative  values  resulted)  in 
the  high-temperature  rhombohedral  phase  at  temperatures  near  Tq  and  at  compositions  close 
to  the  morphotropic  boundary.  However,  away  from  this  small  region  the  dielectric 
properties  could  be  calculated  with  good  agreement  with  the  available  experimental  dqra 

The  present  theory  could  be  improved  as  additional  experimental  data  becomes 
available.  In  particular,  additional  spontaneous  strain  data  on  very  homogeneous  powders 
are  needed  to  more  precisely  locate  the  tricritical  point  on  the  lead  titanate  side  of  the  phase 
diagram.  Strain  data  on  compositions  close  to  lead  zirconate  are  also  needed  to  determine  the 
compositional  dependence  of  the  antiferroelectric  related  coefficients.  As  described  above 
additional  tilt  angle  data  are  needed  to  more  accurately  determine  the  tilt  angle  related 
coefficients. 

In  this  series  of  papers  a  single  energy  function  and  set  of  coefficients  was  presented 
that  can  be  used  to  model  all  of  the  known  phase  transitions  in  the  PZT  solid  solution 
system.  This  theory  can  also  be  used  to  calculate  the  single  domain  properties  of  PZT.  In 
addition  to  the  properties  calculated  in  this  thesis,  the  theory  can  be  used  to  calculate  the 


1  7 


elastic  compliances  at  constant  field,  the  dielectric  susceptibilities  at  constant  strain  (the 
dielectric  susceptibilities  calculated  in  this  paper  were  at  constant  stress),  and  thermal 
properties  such  as  the  entropy  change  and  latent  heat  at  a  transition.  26.27  can 

now  be  used  for  several  different  applications. 

One  application  that  this  theory  has  already  been  used  for  was  to  help  explain  the 
large  electromechanical  anisotropy  that  occurs  in  PbTiC>3  ceramics,  but  not  in  the  single 
crystal  properties.  By  averaging  the  single-crystal  constants  obtained  from  this  theory,  a 
large  electromechanical  anisotropy  was  found  to  occur  due  to  the  particular  values  of  the 
single-crystal  electrostnctive  ratios  of  PbTiC^.23  Other  compositions  in  ceramic  form,  such 
as  PZT  50/50,  have  small  electromechanical  anisotropy,  because  the  values  of  the  single 
crystal  electrostrictive  ratios  have  changed. 

The  intrinsic  and  extrinsic  contributions  to  the  properties  of  polycrystalline  PZT 
ceramics  can  now  be  separated  with  the  results  of  this  theory.  By  using  simple  averaging 
relations  the  intrinsic  contributions  to  the  ceramic  properties  can  be  calculated  from  the 
theoretical  single-domain  properties.  By  then  comparing  these  properties  with  the 
experimental  measurements  on  ceramic  samples,  the  extrinsic  contributions  can  be 
determined.  Comparisons  with  low-temperature  and  high-frequency  measurements  on 


ceramics  will  be  particularly  useful  in  seeing  where  these  extrinsic  contributions  "freeze  out ' 
at  low  temperature  and  "relax  out"  at  high  frequencies.  After  separating  the  extrinsic 
contributions,  comparisons  of  these  results  should  be  made,  with  the  calculations  from  other 
theories,  such  as  that  by  Arlt  and  Peusens,2^  where  the  domain  wall  contributions  to  the 
dielectric  constant  of  BaTi03  ceramics  were  calculated. 

This  theory  can  also  be  used  to  study  the  effects  of  electrical  and  mechanical 
boundary  conditions  on  the  properties  and  phase  stability,  which  should  also  lead  to  a  better 
understanding  of  ferroelectric  polycrystailine  materials  where  the  boundary  conditions  may 
vary  from  grain  to  grain.  Amin  et  al.  ^,2^  did  study  these  effects  using  the  earlier 
phenomenological  theory  that  was  developed  for  the  morphotropic  phase  boundary  region.  '^ 


They  found  that  the  application  of  an  electric  field  would  easily  transform  the  rhombohedral 
state  to  the  tetragonal  state,  but  that  by  applying  a  field  to  the  tetragonal  state  would  not 
transform  it  to  the  rhombohedral  state.  They  then  concluded  "that  this  may  be  pan  of  the 
reason  there  is  a  rapid  escalation  of  coercivity  against  poling  in  the  tetragonal  phase 
compositions."  Amin  et  al. also  studied  the  effect  of  hydrostatic  pressure  on  the  properties 
and  phase  stability,  and  found  that  for  morphotropic  boundary  compositions  the 
rhombohedral  phase  can  easily  be  transformed  to  the  tetragonal  phase  from  hydrostatic 
pressure.  Additional  studies  of  the  effects  of  electrical  and  mechanical  boundary  conditions 
should  now  be  continued  using  the  more  complete  theory  of  PZT  that  has  now  been 
developed. 

The  energy  function  developed  for  PZT  was  based  on  the  elastic  Gibbs  free  energy, 
which  resulted  in  dielectric  stiffness  coefficients  under  constant  stress  conditions.  If  the 
Helmholtz  free  energy  had  been  used,  then  the  dielectric  stiffness  coefficients  would  be 
under  constant  strain  conditions.  The  fourth-order  dielectric  stiffness  coefficients  from  these 
two  energy  functions  can  be  related  through  terms  involving  the  electrostrictive  and  elastic 
constants.26  -3 1  Since  the  fourth-order  dielectric  stiffness  coefficients  changed  signs  and 
caused  the  tricritical  behavior  to  occur,  it  would  now  be  useful  to  compare  these  coefficients 
with  the  coefficients  from  the  Helmoitz  function  to  see  if  the  tncritical  behavior  is  due  to  the 
electrostrictive  coupling,  or  due  to  the  intrinsic  behavior  of  the  material.  Also  by  comparing 
the  coefficients  of  these  energy  functions  additional  understanding  of  the  phase  stability  may 
be  gained. 

With  the  phenomenological  theory  now  developed  for  the  PZT  system  extensions 
into  more  complex  systems,  such  as  lanthanum  modified  PZT  (PLZT),  can  now  be  started. 
This  should  be  very  important,  because  for  most  technological  applications  PZT  is  modified 
with  other  ions.  Extending  the  PZT  theory  into  the  PLZT  system  should  also  be  useful  in 
furthering  the  understanding  relaxor  type  ferroelectric  materials. 
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In  addition  to  the  applications  described  above,  the  phenomenological  theory  of  PZT 
will  provide  a  method  of  collecting  ail  of  the  dielectric,  piezoelectric,  elastic,  and  thermal  data 
into  an  organized  and  easily  tractable  form. 
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LIST  OF  FIGURES  AND  TABLES 

Figure  1.  The  product  of  the  free  energy  AG  and  Curie  constant  plotted  versus 
composition  for  the  different  solutions  of  the  energy  function. 

Figure  2.  Superposition  of  the  theoretical  and  experimental  phase  diagrams.  The  data 
points  are  from  the  experimental  phase  diagram. ^  The  solid  curves  are  the 
theoretical  calcularions. 

Figure  3.  The  theoretical  spontaneous  polarization  P3  plotted  versus  temperature  for  Five 
tetragonal  compositions. 

Figure  4.  The  theoretical  resultant  spontaneous  polarization  P$  (3  ^  P3)  plotted  versus 
temperature  for  four  compositions  with  the  high  to  low  temperature 
rhombohedral  phase  transition. 

Figure  5.  The  resultant  spontaneous  tilt  angle  65  (3l/2 

63)  plotted  versus  temperature. 

The  data  points  are  from  References  1 1  and  12.  The  solid  curves  are  the 
theoretical  calculations. 

Figure  6.  The  theoretical  resultant  spontaneous  tilt  angle  0$  (3^  63)  plotted  versus 
composition. 

Figure  7.  The  theoretical  spontaneous  polarization  plotted  versus  composition.  The 
curves  on  the  rhombohedral  side  of  the  figure  represent  the  resultant 
spontaneous  polarization  P5  (3  ^  P3). 

Figure  8.  The  theoretical  dielectric  susceptibility  coefficients  of  PZT  50/50  plotted  versus 
temperature. 

Figure  9.  The  theoretical  dielectric  susceptibility  coefficients  plotted  versus  composition 
at  25  °C. 

Figure  10.  The  theoretical  dielectric  susceptibility  coefficients  of  PZT  60/40  plotted  versus 
temperature. 

Figure  1 1 .  The  theoretical  dielectric  susceptibility  coefficients  of  PZT  90/10  plotted  versus 


temperature 


24 


Figure  12. 


Figure  13. 


Figure  14. 


Figure  15. 


Figure  16. 


Figure  17. 


Figure  18. 


Figure  19. 


The  theoretical  piezoelectric  voltage  coefficients  gjj  of  PZT  50/50  plotted  versus 
temperature. 

The  theoretical  piezoelectric  voltage  coefficients  gjj  of  PZT  60/40  plotted  versus 
temperarure. 

The  theoretical  piezoelectric  voltage  coefficients  gy  of  PZT  90/10  plotted  versus 
temperature. 

The  theoretical  piezoelectric  voltage  coefficients  gjj  plotted  versus  composition 
at  25  °C. 

The  theoretical  piezoelectric  charge  coefficients  djj  of  PZT  50/50  plotted  versus 
temperature. 

The  theoretical  piezoelectric  charge  coefficients  djj  plotted  versus  composition 
at  25  °C. 

The  theoretical  piezoelectric  charge  coefficients  djj  of  PZT  60/40  plotted  versus 
temperature. 

The  theoretical  piezoelectric  charge  coefficients  djj  of  PZT  90/10  plotted  versus 
temperature. 


Table  I.  The  Theoretical  Properties  of  PZT  at  25  °C. 


Table  l.  The  Theoretical  Properties  of  PZT  at  25  °C. 


Mole  Fraction  PbTiOj  in  PZT 


CLL 

01 

01 

OA 

01 

M 

01 

<LS 

01 

j_0 

Ps  (C/m2) 

0.57 

0.66 

0.65 

0.50 

0.50 

0.57 

0.64 

0.70 

074 

0.75 

0S  (Deg.) 

5.42 

7.36 

6.64 

— 

— 

— 

— 

— 

— 

— 

n33  or  1333’ 

219. 

111. 

280. 

295. 

382. 

197. 

116. 

86.4 

72.9 

66  6 

nn  °rTiii 

195. 

262. 

360. 

529. 

1721. 

498. 

218. 

143. 

121. 

124 

g33  (10-3  Vm/N) 

26.8 

30.9 

34.6 

41.9 

96.6 

92.9 

101. 

114. 

124. 

134. 

g31  (10'3  Vm/N) 

-  8.13 

-9.53 

-  11.2 

-  15.6 

-46.0 

-33.8 

-31.8 

-34.2 

-37.1 

-  39  2 

gl5(lO*3  Vm/N) 

12.9 

14.3 

15.5 

18.1 

41.0 

38.4 

40.8 

44.9 

48.6 

50.9 

d33(10-12C/N) 

47.2 

72.1 

107. 

189. 

327. 

162. 

104. 

87.2 

81.2 

79.2 

d3l(10*l2ON) 

-  13.3 

-  21.6 

-38.5 

-80.5 

-  156. 

-58.9 

-32.6 

-26.2 

-23.9 

-23.1 

d,4(l0-12C/N) 

1.81 

1.27 

-7.26 

-25.1 

— 

— 

— 

— 

— 

— 

d15(10-l2C/N) 

24.2 

34.5 

42.0 

60.0 

624. 

169. 

78.6 

57.0 

52.1 

56.1 

*The  values  listed  from  0. 1  to  0.4  were  calculated  from  the  high- temperature  rbombohedral  equanons.  except 
for  P<j  and  85  from  0. 1  to  0.3  which  were  calculated  from  the  low- temperature  rbombohedral  equanons.  The 
values  listed  Grom  0.5  to  1.0  were  calculated  Grom  the  tetragonal  equations.  P5  *  P3  in  the  tetragonal  state, 
but  P5  *  3 P3  in  the  rbombohedral  states.  85  *  3^  83.  n  j  j  and  ^33  were  used  in  the  tetragonal  state, 
and  q(  and  1^33'  (see  Section  2.5)  were  used  in  the  rbombohedral  state.  The  piezoelectric  coefficients  are  all 
based  on  the  cubic  axes. 
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THERMODYNAMIC  THEORY  OF  PbZr03 


MJ.  Haun,t  T.J.  Harvin,  M.T.  Lanagan.tt 
Z.Q.  Zhuang.ttt  s.J.  Jang,  and  L.E.  Cross 

Materials  Research  Laboratory 
The  Pennsylvania  State  University 
University  Park,  PA  16802 

A  thermodynamic  theory  is  presented  to  model  the  phase  transitions  and  properties  of  lead 
zirconate.  A  two-sublattice  model  is  used  to  account  for  the  antiferroelectric  and  ferroelectric 
phases. 

The  free  energy  AG  is  expressed  as  a  power  series  in  terms  of  the  ferroelectric  polarization 
(Pi  ss  Pai  +  PbO  and  antiferroelectric  polarization  (pi  =  Pai  -  Pbi)  including  all  possible  terms  up  to 
sixth  power,  but  only  first  order  cross  coupling  terms  and  couplings  to  elastic  stress.  Under  the 
assumption  that  only  the  lowest  order  coefficients  of  Pi  and  pi  are  linearly  temperature  dependent 
(Curie  Weiss  behavior)  and  all  other  constants  are  temperature  independent  experimental  data  are 
used  to  define  the  constants  and  permit  calculation  of  ferroelectric  and  antiferroelectric  free  energies 
as  a  function  of  temperature. 

Use  of  the  function  to  define  the  averaged  dielectric  permittivity  at  room  temperature  in  the 
antiferroelectric  phase  gives  a  value  of  6r  =  120  in  good  agreement  with  recent  microwave 
measurements. 

A  simplified  technique  for  modifying  the  function  to  explore  solid  solution  with  lead 
titanate  is  examined  and  shown  to  lead  to  excellent  agreement  with  the  known  phase  diagram. 


^Now  at  DuPont  Experimental  Station,  Electronics  Dept.  Bldg.  334,  Wilmington,  DE  19898. 
t^Now  at  Argonne  National  Lab.,  Materials  and  Components  Technology  Div.,  Argonne,  IL 
60439. 

^Visiting  Scientist  from  the  Department  of  Inorganic  Materials  Science  and  Engineering,  South 
China  Institute  of  Technology,  Guangzhou,  The  People's  Republic  of  China. 


I.  INTRODUCTION 


Lead  zirconate  is  an  end  member  of  the  technologically  important  lead 
zirconate -titanate  (P2T)  solid  solution  system.*  At  room  temperature  PbZtOj  has  an 
antiferroelectric  orthorhombic  (Aq)  perovskite  structure  with  an  antipoiar  arrangement  along 
the  [l  10]  direction.  The  Aq  phase  remains  stable  up  to  =  220  °C,  where  a  transition  occurs 
to  a  ferroelectric  rhombohedral  ( Fr)  phase  with  a  polarization  along  the  [  1 1 1  ]  direction.  - 
The  Fr  phase  is  only  stable  over  a  narrow  temperature  range,  and  transforms  to  a 
paraelectric  cubic  (Pq)  phase  at  =  232  °C.^ 

Lead  zirconate  has  been  recently  studied  as  a  possible  dielectric  material  for  high 
frequency  applications.-*’^  In  this  study  a  dielectric  relax  anon  was  found  to  occur  at 
microwave  frequences.  To  further  the  understanding  of  the  dielectric  properties  of  lead 
zirconate,  a  thermodynamic  theory  has  been  developed  to  calculate  the  intrinsic  dielectric 
response.  The  development  of  a  thermodynamic  theory  of  lead  zirconate  was  also  needed  to 
complete  the  theory  that  was  developed  for  the  PZT  solid  solution  system.  ^  *0 

Whatmore  and  Glazer*  *  used  a  one -dimensional  energy  function  to  model  the  Pb 
ion  displacement  in  lead  zirconate.  Uchino  et  al.*^  included  stress  terms  in  the  energy 
function  to  determine  the  hydrostatic  electrostrictive  coefficient.  However,  in  these  papers 
not  enough  coefficients  were  determined  to  calculate  the  energies  and  dielectric  properties  of 
the  phases. 

In  this  paper  a  more  complete  thermodynamic  theory  of  lead  zirconate  will  be 
developed  using  a  three-dimensional  energy  function.  In  the  next  section  the  energy  function 
will  be  presented,  along  with  the  solutions  and  property  relations  that  can  be  derived  from 
this  energy  function.  Values  of  the  coefficients  will  be  determined  from  experimental  data  in 
Section  ID,  and  used  to  calculate  the  theoretical  properties  in  Section  IV.  The  calculated 
dielectric  properties  will  be  compared  with  experimental  high-frequency  data.  In  Section  V  a 
procedure  will  be  presented  to  calculate  the  energy  of  the  antiferroelectric  orthorhombic  phase 
into  the  PZT  system.  Finally  a  summary  of  this  paper  will  be  given  in  Section  VI. 
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0.  PHENOMENOLOGICAL  THERMODYNAMIC  THEORY 


The  following  energy  function  for  lead  zirconate  was  derived  from  the  rwo-sublattice 
theory  that  was  previously  developed  to  model  antiferroelectric  materials. 1 3- 1 6  The 


coefficients  were  limited  by  the  symmetry  of  the  paraelectric  phase  (m3m).  By  assuming 
isothermal  conditions,  using  reduced  notation,  and  expanding  the  energy  function  in  powers 
of  the  ferroelectric  (P,)  and  antiferroelectric  (pj)  polarizations,  and  including  couplings 
between  these  order  parameters,  and  between  the  stress  (Xrn)  and  the  order  parameters  the 
following  energy  function  resulted: 

AG  =  al[Pl2  +  P22  +  P32]  +  all[Pl4  +  P24  +  P34] 

+  al2  [P^  P22  +  P22  P32  +  P32  Pi2]  +  «i i  i[Pi6  +  P26  +  P36] 

+  all2  tPl4  (p22  +  P32)  +  p24  Cpl2  +  p32)  +  P34  (P[2  +  P22)] 

+  a123pl2p22p32  +  a\  [Pi2  +  P22  +  P32l  +  all  [Pl4  +  P24  +  P34l 
+  °12  lPl2  P22  +  P22  P32  +  P32  Pi21  +  <*l 1 1  tPl6  +  P26  +  P36] 

+  Oi  12  [pi4  (p22  +  p32)  +  p24  (Pl2  +  p32)  +  p34  (Pl2  +  p22)] 

+  al23  Pl2  P22  P32  +  ^11  Pl2  Pl2  +  p22  P22  +  p32  P32l 
♦  M-12  ^l2  (P22  +  P32)  +  p22  ^Pl2  +  P32)  +  P32  (Pl2  +  P22)l 
+  M-44  fP lp2PlP2  +  P2P3P2P3  +  P3PlP3Pl) 

-  1/2  S  i  i [X 1 2+  X22+  X32]  -  S  i 2  [X  i X2+  X2X3+  X3X 1  ] 

-  1/2  S44  [X42+  X62+  X62]  -  Q i i  [XiPi2  +  X2P22  +  X3P32] 

-  Ql2  [Xi  (P22  +  P32)  +  X2  (Pi2  +  P32)  +  X3  (Pi2  +  P22)] 

-  Q44  [X4  P2P3  +  X5  plp3  +  x6  plp2l  -  Z1 1  txlPl2  +  X2P22  +  X3P321 

"  zi2  (xi  (p22  +  P32)  +  X2  (pi2  +  p32)  +  X3  (pi2  +  p22)] 

-  Z44  [X4p2p3  +  X5  pjp3  +  X$  pip2j  (1) 


The  coefficients  of  this  energy  function  are  defined  in  Table  I.  The  energy  function  includes 
all  possible  ferroelectric  and  antiferroelectric  polarization  terms  up  to  the  sixth  order,  and  the 


first -order  coupling  terms. 

The  following  solutions  to  the  energy  function  (Equation  l)  are  of  interest  in  the  lead 
zirconaie  system: 


Paraelectric  Cubic  (P^) 

p1  =  p2  =  p3  =  o. 

?!  =p2  =  P3  =  °, 

(2) 

Ferroelectric  Rhombohedral  (Fr) 

Pl2  =  P22  =  p32  Pi  =P2  =  P3  =  0- 

(3) 

Antiferroelectric  Orthorhombic  1 A  q) 

Pl=P2  =  P3=0. 

Pl  =o.  p22  =  p32*o 

(4) 

Applying  these  solutions  to  Equation  ( 1 )  under  zero  stress  conditions  results  in  the 
following  relations  for  the  energies  of  each  solution: 


Pc 

o 

H 

3 

(5) 

PR 

AG  -  3  ot[P3^  +  3  (<Xj  ai2)  P34  +  (3  aj  1 1+  6  aj  j2  +  al23^  p3^ 

(6) 

A0 

AG  =  2  a!P32  +  (2  !+  a12)  p34  +  2  (aL  i  12)  p36 

(7) 

The  spontaneous  ferroelectric  and  antiferroelectric  polarizations  (P3  and  p3)  in  the 

above  equations  can  be  found  from  the  first  partial  derivative  stability  conditions  (dAG/9P3 

and  dAG/dp3)  as  shown  below: 

PR 

8AG/dP3  =0  =  (3am+6a|[2  +  al23)  p34  +  2  (al  1+  a12>  p32  +  al 

(8) 

A0 

3AG/dp3  =  0  =  3(am+an2)  p34  +  (2  <Ji  i+  <?i2)  P32  +  <*i 

(9) 

The  polarizations  can  be  calculated  by  solving  these  quadratic  equations.  Equations  (6)  -  (9) 
relate  the  energies  of  each  solution  to  the  coefficients  of  the  energy  function.  Thus  by 
determining  these  coefficients,  the  energies  of  each  phase  can  be  calculated. 


x* mm  TrsannrsrBnwmr  mg  'rrrivm  r:  .■  v. 


.  ,.  (as  a^c/ax,!  under  zero  stress  conditions  can  be 

The  spontaneous  elastic  strains  i 

derived  from  Equation  (1)  as  follows: 

(10) 

Pc  xi  =  *2  =  *3  =  x4  =  *5  =  *6  =  ° 

rR  M-2-3-‘Qn^Qi2>^-  *«-*s**6-<lW 

A.  =  ll2> 

,  ■  t;ffn*<ises  v •  •  (=  a2 AG/aPiaPj)  were  derived 

Relations  for  the  relative  dielectnc  stiffnesses  Xy '  J 

from  Equation  (1)  for  the  solutions: 

(13) 

PC  Xlls:X22  =  X33s2E0al'  Xl2  =  X23  =  X31-°  ^ 

,  xu=M2  =  X33-eol«l^-ll-al2)P,na5aul  +  .4au2  +  a123)P3l. 

*  X12 * X23 * X31  “ 4 €o  tal2 ^32  +  3f  all2  +  a123  ^  ^3 

-  _  2i  Y^  =  Y'zx  =  2£olal  +  ^ll  +  JAl2)P3 

Aq  XUs2£0tal  +  2jtl2P3  ’  X22  2X33 
Xl2  =  X3l=0'  X23s£0^44P3 


(15) 


The  multiplication  by  the  permittivity  of  fee  space  £<,  »  *  uscd 

MRV,„  h„m  absolute  ,0  —  — *  — -  *—  axes 

alculate  the  relaove  dtelecntc  stresses  for  each  phase  based  ou  .he  onguul  cub.c  axes, 
to  calculate  tne  reian  <■ ,  <  n>  directions 

In  die  orthorhombic  state  the  polarization  can  be  along  any  o 

..  The  polarization  of  the  thombohedial  state  can  be  along  any  of 

of  the  original  cubtc  axes,  the  polanaat.  the 

a  axes  so  that  for  both  states  the  new  x3  axis  is  aio  * 

1 1>  directions.  By  rotating  these 

u  The  new  dielectric  stiffness  coefficients 
polar  directions,  diagonalized  matrices  will  result.  Tl* 

a  h  nrime)  can  be  related  to  the  old  coefficients  ^defined  by  Equations  l  Id)  and 
(indicated  by  a  prune)  can  oe  rci* 

(15))  with  the  following  relations: 


Fr 


Xll'  =  X22  =  XU  *Xl2- 
X12  =  223'  =  231  =0 


233  =  Xll  +  2  2.12 


(16) 


v'  V  V 


aataaa  m  aaaai  riaeawi  i 


Xu'=Xll.  X22  =  X33  ‘  X23 

X33'  =  X33  +  X23-  X 12  =  X23  =  X31  =  0 
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A0 


(17) 


These  equations  can  be  used  to  calculate  the  dielectric  stiffnesses  of  the  orthorhombic  and 
rhombohedral  phases  parallel  and  perpendicular  to  the  polar  axes. 


UI.  EVALUATION  OF  THE  COEFFICIENTS 

In  this  section  coefficients  of  the  energy  function  will  be  determined  from 
experimental  data.  All  of  the  coefficients  were  assumed  to  be  independent  of  temperature, 
except  the  ferroelectric  (a^)  and  antiferroelectric  (Oj)  dielectric  stiffness  coefficients  which 
were  given  a  linear  temperanire  dependence  based  on  the  Curie-Weiss  law: 

T -  0  T - 0A 

ai  =  -  (18)  CTi  =  -  (19) 

2  £0  C  2  EO  CA 

C  is  the  Curie  constant,  £q  is  the  permittivity  of  free  space.  9  is  the  Curie-Weiss  temperature, 
and  CA  and  0A  are  antiferroelectric  constants  analogous  to  C  and  0. 

By  combining  Equations  (13)  and  (18).  C  and  9  can  be  determined  from  a  linear  fit 
of  experimental  dielectric  stiffness  data  in  the  paraelectric  cubic  state  (the  Curie-Weiss  law). 
The  dielectric  stiffness  is  the  inverse  of  the  dielectric  susceptibility  rj y,  which  will  be 
assumed  to  be  equal  to  the  relative  dielectric  permittivity  £jj  (actually  £tj  =  rj-j  +  1 ).  Values  of 
C  and  8  for  lead  zirconate  were  previously  found  to  be  equal  to  1.5  x  10^  °C  and  190  °C  by 
r rang  the  Cune- Weiss  law  to  single-cTystal  dielectric  data.2,2  These  constants  will  be  used 
10  calculate  the  ot|  coefficient  versus  temperature.  A  different  procedure  was  used  to 
determine  the  antiferroelectric  CA  and  9 A  constants  as  will  be  described  later  in  this  section. 

The  ferroelectric  rhombohedral  spontaneous  polarization  (P3)  can  be  related  to  the 
coefficients  of  the  energy  function  by  solving  the  quadratic  relation  formed  from  the  first 
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partial  derivative  stability  condition  [Equation  (8)]: 

.  -4  +  K2-9oc14)l/2 

p  2  =  -  ,  (20) 

where  4  =  3  (a^  +  a^).  and  4  =  3  04 1 1  +  6  12  +  a123  (21) 

At  T(2>  the  transition  temperature  between  the  ferroelectric  rhombohedral  and 
paraelectric  cubic  phases,  two  relations  must  be  satisfied: 

0  =  3aic  +  4P3c2  +  4P3C4  •  ^d  (22) 

0  =  alc  +  2/3  4  P3C2  +  4  P3C4  ,  (23) 

where  and  P3£  are  ot(  and  P3  at  Tq  Equation  (22)  was  derived  from  the  requirement 
that  the  AG's  of  the  cubic  and  rhombohedral  phases  [Equations  (5)  and  (6)]  must  be  equal  at 
T£.  Equation  (23)  is  the  first  partial  derivative  stability  condition  [Equation  (8)],  which  must 
be  satisfied  so  that  the  stable  state  corresponds  to  the  minima  of  the  energy  function. 

When  the  transition  at  is  first  order  the  spontaneous  polarization  in  the 
ferroelectric  state  will  develop  discontinuously  at  the  transition,  and  thus  P3£  will  be 
nonzero.  In  this  case  Equation  (18)  can  be  substituted  into  Equations  (22)  and  (23)  to  obtain 
relations  for  the  4  and  4  coefficients: 

-  3  (Tc  -  9)  3  (Tc  -  0) 

4  =  - 2  •  5  =  - \  (24) 

«0  c  P3C2  2  £0  C  P3C4 


t 


) 


> 


> 


By  assuming  that  the  4  and  4  coefficients  are  independent  of  temperature,  and  then 
substituting  Equations  (24)  and  (18)  into  Equation  (20).  the  following  relation  results  for  the 
spontaneous  polarization  of  the  ferroelectric  rhombohedral  phase: 

3  (  T  -  9) 


P32  =  v|/P3^2  where  y  =  2  {  1+ [  1  - 


4  (Tc  -  0) 


]l/2} 


(25) 


'M 


s* 

3 


Now  if  P3£,  T^,  and  0  can  be  determined,  the  spontaneous  polarization  of  the  ferroelectric 
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rhombohedral  phase  can  be  calculated  versus  temperature.  Note  that  P3  is  the  X3  component 
of  the  resultant  spontaneous  polarization  (Pg)  along  the  [1 1 1 1  direction,  and  thus 
PS  =  31«P3. 

Using  the  value  of  0  given  above,  with  equal  to  232  °C,^  a  value  of  Pj£  was 
found  from  the  best  least-squares  fit  of  Equation  (25)  to  experimental  spontaneous 
polarization  data  (from  Curve  2  in  Figure  3  of  Reference  17),  as  shown  in  the  insen  of 
Figure  l.  The  resulting  P3C  value  is  listed  in  Table  H  with  the  values  of  the  other  constants 
that  were  used  in  the  calculations.  These  values  were  used  to  calculate  the  ferroelectric 
rhombohedral  spontaneous  polarization  versus  temperature  down  to  -  273  °C  as  shown  in 
Figure  1.  However,  over  most  of  this  temperature  region  the  rhombohedral  phase  is 
metastable  to  the  stable  antifenroelectric  orthorhombic  phase. 

A  similar  relation  to  Equation  (25)  can  be  derived  for  the  spontaneous  strain  X4  by 
substituting  Equation  (25)  into  Equation  (11): 


x4  =  V  X4C  ’  wr*re  X4C  *  044  P3C2 


(26) 


This  equation  was  used  with  the  values  of  the  constants  listed  in  Table  II  to  determine  a  value 
of  X4C  (X4  at  T(-)  that  gave  the  best  least-squares  fit  of  experimental  X4  data,  as  shown  in  the 
insert  of  Figure  2.  The  experimental  X4  data  was  calculated  from  rhombohedral  angle  otR 
data  from  Reference  1 1  using  the  relation:  X4  =  (90  -  <Xr)/90.  A  value  of  the  electrostrictive 
Q44  constant  was  calculated  from  x^  and  P3C  using  Equation  (26).  This  value  is  listed  in 
Table  II.  The  ferroelectric  rhombohedral  spontaneous  strain  X4  was  calculated  versus 
temperature  into  the  andferroelectnc  region  as  shown  in  Figure  2  using  Equation  (26)  and  the 
constants  listed  in  Table  Q. 

The  constants  determined  above  will  be  used  to  calculate  the  AG  of  the  ferroelectric 
rhombohedral  phase  in  the  next  section.  The  procedure  used  to  determine  values  of  the 
antiferroelectric  constants,  which  are  needed  to  calculate  the  AG  of  the  antiferToelectric 
orthorhombic  phase,  will  now  be  presented. 
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The  anciferroelectric  measure  of  the  polarization  P3  for  the  orthorhombic  solution  can 
be  related  to  the  antiferroelectnc  a  coefficients  by  solving  the  quadratic  relation  formed  from 
the  first  partial  derivative  stability  condition  [Equation  (.9)1: 

-  5  +  [52  -6  ai  YlI/2 

P3‘  =  -  (27> 

3  Y 

where  5  =  2ct11+<Jl2.  and  Y=  2(arlll+<Tll2*  (28) 

The  Neel  temperature  (Tjq)  is  the  transition  temperature  from  an  antiferroelectnc 
phase  to  a  paraelectric  phase.  In  lead  zirconate  this  is  a  metastable  phase  transition,  because 
it  occurs  in  the  ferroelectric  rhombohedral  stability  region.  As  in  the  rhombohedral  case  at 
T^,  two  equations  must  be  satisfied  at  Tjq: 

0  =  2alN  +  Sp3N2  +  YP3N4  (29) 

0  =  o^n  +  5p3N^  +  3/2  YP3N4  130) 

(Tin  P3N  ^  <T[and  P3  at  Tjq-  Equation  (29)  was  derived  from  the  requirement  that  the 
AG's  of  the  two  phases  [Equations  (5)  and  (7)]  must  be  equal  at  Tjq.  Equation  (30)  is  the 
first  partial  derivative  stability  condition  [Equation  (9)]. 

Substituting  Equation  ( 19)  into  Equations  (29)  and  (30).  and  solving  for  the  5  and '( 
coefficients  results  in: 


-2(Tn-9a) 
e0  CA  P3N2 


tn  ~9a 

£0  CA  P3N4 


(31) 


where  P3N  is  P3  at  Tjq 


1  0 

By  assuming  that  the  5  and  y  coefficients  are  independent  of  temperature,  and  then 
substituting  Equations  (19)  and  (3 1 )  into  Equation  (27)  yields: 

*  2  r  r  3  ( T  -  9  a  ) 

P32  =  VP3N  .  where  V'aM+L1 - ]  1/2  }  (32) 

3  4(TN-eA) 

Unfortunately,  there  is  not  any  experimental  antiferroelectric  polarization  pj  data  available, 
and  thus  p3pj  will  not  be  determined  in  this  paper. 

However,  antiferroelectnc  spontaneous  strain  data  can  be  calculated  from  cell 
constant  data,  and  used  to  determine  two  of  the  antiferroelectric  constants.  By  substituting 
Equation  (32)  into  Equation  ( 12)  the  following  relation  results  for  the  antiferroelectric 
spontaneous  strain  x4: 

x4  =  yx4N,  where  x4N  =  P3N2  (33) 

x4n  is  x4  at  Tjq.  Values  of  Tjq,  9^,  and  x4pj  are  needed  to  calculate  x4  using  Equations 
(32)  and  (33). 

Spontaneous  strain  x4  data  was  calculated  from  cell  constant  data  from  Reference  1 1 
[interaxial  Yp  data  from  Figure  5  in  Reference  1 1  was  related  to  x4  with  the  relation:  x4  =  (90 
-  Yp)/90j.  By  fitting  this  data  with  all  three  constants  (Tjq,  9^,  and  x4pj)  as  unknowns 
many  combinations  of  the  values  of  these  constants  gave  similar  fits  of  the  data.  For  this 
reason  the  value  of  one  of  these  constants  had  to  be  determined  from  additional  data. 

Without  additional  data,  it  was  necessary  to  make  an  assumption. 

Tpj  must  occur  between  the  antiferroelectric-ferroelectric  transition  temperature 
(T Ap?)  and  the  ferroelectric-paraelectric  transition  temperature  {Tq).  T ^  occurs  at  =  220  °C 
(average  of  the  hearing  and  cooling  DTA  data  in  Reference  2),  and  T^  at  =  232  °C.  Thus 
since  Tp^  must  occur  between  220  and  232  °C,  it  was  assumed  to  be  equal  to  226°C,  the 
average  of  the  two  extremes.  This  was  a  good  assumption,  because  the  constants  were  not 
significantly  affected  even  if  T-q  was  varied  a  few  degrees.  With  Tjq  fixed  at  226  °C  values 


of  9 a  and  were  found  that  gave  the  best  least-squares  fit  of  the  X4  data  as  shown  in 

Figure  3.  These  values  are  listed  in  Table  H. 

Substituting  Equations  ( 19),  (28).  (31),  (32)  and  (33)  into  Equation  (7)  results  in  the 
following  relation  for  the  AG  of  the  antiferroelectric  orthorhombic  phase: 

AG  =  y  *4N^e0  Z44  C^)  (T  -  T^  -2  \\f  (Tjq  -  T^)  (1-1/2  y)J  (34) 

All  of  the  constants  needed  to  calculate  the  AG  of  the  antiferroelectric  orthorhombic  phase 
using  this  equation  have  been  determined,  except  for  the  Z44  and  CA  constants.  At  this 
point  there  is  not  enough  experimental  data  available  to  determine  values  of  both  of  these 
constants.  However,  the  product  of  these  constants  can  be  determined  by  equating  the  AG's 
of  the  Aq  and  phases  at  the  transition  between  these  phases: 

At  Taj;:  AGFr  =  AGAo  (35) 

The  rhombohedral  AG  at  Tap  can  be  calculated  from  Equations  (6),  (19),  (21 ),  (24),  and 
(25)  using  the  constants  listed  in  Table  H.  This  value  can  then  be  used  with  Equations  (34) 
and  (35)  and  the  constants  from  Table  II  to  calculate  a  value  of  the  Z44  Ca  product.  This 
value  is  listed  in  Table  II.  By  assuming  that  Z44  Ca  is  independent  of  temperature,  all  of  the 
constants  necessary  to  calculate  the  AG  of  the  antiferroelectric  orthorhombic  phase  versus 
temperature  have  been  determined.  Calculations  of  the  AG’s  of  the  ferroelectric  and 
antiferroelectric  phases  will  be  presented  in  the  next  section. 

All  of  the  a  coefficients  must  be  determined  to  calculate  the  dielectric  stiffness  Xjj 
coefficients  of  the  ferroelectric  rhombohedral  phase  [see  Equation  ( 14)J.  The  constants 
determined  earlier  in  this  section  can  be  used  in  Equation  (24)  to  calculate  the  £ 

[=  3  (<Xj  j  +  aj2)J  and  ^  [=  3  aj  j  j  +  6  aj  12  +  otj 23]  coefficients,  which  can  be  used  to 
calculate  the  spontaneous  polarization,  strain,  and  AG  of  the  Fr  phase.  At  this  tune  there  are 
not  enough  lead  zirconate  data  present  to  determine  values  of  aj  j,  a^.  ai  u.  04  12*  and 
01 123  However,  these  coefficients  were  determined  in  the  PZT  system,  and  extrapolated  to 
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lead  zirconate.^'^ 

To  calculate  the  dielectric  stiffness  Xij  coefficients  of  the  antiferroelecmc 
orthorhombic  phase,  the  |Xjj  antiferroelectric-ferroelectric  coupling  coefficients  must  be 
determined.  Unfortunately  the  data  necessary  to  determine  these  coefficients  has  not  been 
measured.  One  goal  of  this  project  was  to  compare  calculations  from  this  theory  with 
polycrystalline  high-frequency  dielectric  data.  To  accomplish  this  a  procedure  was  used  to 
combine  the  unknown  constants  in  the  dielectric  relations  into  one  constant  that  could  be 
determined  from  low-temperarure  polycrystalline  dielectric  data.  This  procedure  will  be 
described  below. 

By  combining  Equations  (15)  and  (17)  dielectric  stiffness  relations  based  on  the 
normal  orthorhombic  axes  result: 


Xll  =  2e0  lal  +  2M-12P32]  * 

X22  s  2  £0  [“1  +  4*1 1  +  H2  *  M44)  P32l  ’ 

X33'  =  2  £o  tal  +  <Hl  1  +  Hl2  +  i/2  1144)  P32]  (36) 

These  single-crystal  coefficients  were  related  to  the  polycrystalline  dielectric  stiffness  %  la 
bold-faced  italic  symbol  will  be  used  to  refer  to  polycrystalline  constants)  using  the  following 
parallel  averaging  relation: 

X  =  l/3(xu’+  X22'+  X33')  (37) 

By  substituting  Equation  (36)  into  Equation  (37)  the  following  relation  results: 

X  =2E0(al  +  2/3(11!^  2  iil2)P32]  (38) 

Not  enough  experimental  data  was  available  to  determine  the  constants  necessary  to  calculate 
P3.  However,  the  constants  needed  calculate  the  antifertoelectric  spontaneous  strain  X4  were 
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determined.  By  substituting  Equation  (12)  into  Equation  (38)  P3  can  be  replaced  by  x4: 

2  (M-i  1  +  2  IJ-12) 

1/g  a  %  =  2^(04  +  Xx4)  ,  where  X  =  -  (39) 

3  z44 

Since  and  x4  can  be  determined  from  Equations  (18)  and  (33)  using  the  constants  Listed 
in  Table  II,  the  new  constant  X  can  be  calculated  from  an  experimental  polycrystalline 
dielectric  constant  (£)  measurement  (£  =  l/%  )• 

To  provide  experimental  data  to  evaluate  the  X  constant  the  dielectric  constant  was 
measured  at  lKHz  from  -  268.8  °C  to  room  temperature  on  polycrystalline  lead  zirconate. 
The  procedure  used  to  fabricate  the  lead  zirconate  samples  was  described  in  References  4  and 
5.  The  low-temperature  measurement  apparatus  and  technique  that  was  used  was  described 
in  Reference  18.  At  low-temperatures  (-  268.8  °C)  the  extrinsic  contributions  to  the 
polycrystalline  dielectric  properties,  such  as  domain  wall  motion  or  thermally  activated  defect 
motions,  were  assumed  to  "freeze  out".  The  remaining  dielectric  properties  would  then 
represent  the  intrinsic  contribution,  or  averaging  of  the  single-domain  properties. 

A  dielectric  constant  of  95  was  measured  at  -  268.8  °C,  and  used  to  calculate  the 
value  of  the  X  constant  listed  in  Table  Q.  By  assuming  that  this  constant  is  independent  of 
temperature  the  intrinsic  polycrystalline  dielectric  constant  could  be  calculated  versus 
temperature  using  Equation  (39)  and  the  constants  Listed  in  Table  II.  These  calculations  will 
be  compared  with  experimental  high-frequency  dielectric  measurements  in  the  next  section. 

IV.  THEORETICAL  CALCULATIONS 

In  this  section  theoretical  calculations  will  be  made  using  the  equations  from  Section 
II,  and  the  constants  that  were  determined  in  the  Section  HI.  Figure  4  shows  the  energies 
AG  of  the  antiferroelectnc  orthorhombic  and  ferroelectric  high-temperature  rhombohedral 
phases  plotted  versus  temperature  for  lead  zirconate.  Above  T^  (232  °C)  the  energy  of  the 
ferroelectric  rhombohedral  phase  is  positive,  and  thus  the  cubic  state  with  the  reference 


energy  set  to  zero  is  stable.  At  the  energy  of  the  Fr  phase  decreases  to  zero,  and  then 
becomes  negative  below  causing  the  tetragonal  phase  to  become  stable.  At  the 
energy  of  the  amiferroelectric  orthorhombic  phase  is  zero,  and  a  metastable  transition  occurs 
with  the  cubic  phase.  As  the  temperature  decreases  the  the  larger  temperature  dependence  of 
the  AG  of  the  Aq  phase  compared  to  the  phase  causes  the  energies  of  these  phases  to 
become  equal  at  the  transition  T ^j:.  Below  T ^  the  Aq  phase  becomes  stable,  and  remains 
stable  down  to  -  273  °C. 

The  polycry st ailine  dielectric  constant  (£)  of  the  amiferroelectric  orthorhombic  phase 
was  calculated  from  Equation  (39)  using  the  constants  listed  in  Table  II,  and  plotted  in  Figure 
5  versus  temperature  up  to  T^p.  The  dielectric  constant  of  the  paraelectric  cubic  phase  was 
also  calculated  using  Equation  (13),  and  ploned  in  this  figure  above  Tq  The  ferroelectric 
rhombohedral  dielectric  properties  were  not  investigated  in  this  paper,  but  could  be  calculated 
from  the  values  of  the  coefficients  that  were  recently  determined  in  References  6-10. 

The  experimental  data,  that  was  measured  at  1  KHz  from  -  268.8  °C  to  room 
temperature,  is  ploned  in  the  insert  in  Figure  5.  The  data  point  at  -  268.8  °C  was  used  to 
calculate  the  A.  constant,  as  described  in  Section  ED.  By  assuming  that  the  theory  is 
predicting  the  intrinsic  response  of  the  material,  the  difference  that  develops  between  the 
theoretical  calculations  and  experimental  data  as  the  temperature  is  increased  is  due  to  the 
thermally  activated  (extrinsic)  contributions  to  the  polycrystalline  dielectric  constant.  At 
room  temperature  the  theoretical  dielectric  constant  has  a  value  of  120,  compared  to  the 
experimental  value  of  167.  This  indicates  that  72%  of  the  experimentally  measured 
polycrystalline  dielectric  constant  is  due  to  the  intrinsic  averaging  of  the  single-domain 
constants,  while  28%  is  from  extrinsic  contributions. 

Lanagan  et  al."*^  measured  the  dielectric  constant  of  polycrystalline  lead  zirconate 
from  100  Hz  to  26  GHz.  Their  data  are  plotted  versus  the  log  of  the  frequency  in  Figure  6. 
As  can  be  seen  by  these  data  a  relaxation  occurred  at  microwave  frequences.  Above  this 
relaxation  the  dielectric  constant  agrees  remarkably  well  with  the  calculation  from  the  theory 
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presented  in  this  paper,  indicating  that  the  extrinsic  contributions  to  the  dielectric  properties 
have  "relaxed  out”  above  =  10  GHz. 

The  Han  in  Figures  5  and  6  indicate  that  the  extrinsic  contributions  to  the 
polycrystalline  dielectric  constant  of  lead  zirconate  will  "freeze  out"  at  low  temperatures  and 
"relax  out"  at  high  frequences,  and  that  the  remaining  intrinsic  contributions  can  be 
thermodynamically  modeled.  The  theory  provides  a  method  of  connecting  the 
low-temperature  data  to  the  high-frequency  data.  It  would  now  be  useful  to  measure  the 
high-frequency  dielectric  properties  down  to  low  temperatures. 

The  dielectric  constant  at  10  GHz  was  measured  as  a  function  of  temperature  as 
shown  in  Figure  7.  The  1  KHz  data  from  Figure  5  is  also  plotted  in  this  figure  along  with 
the  theoretical  calculations.  Unfortunately  at  this  time  the  10  GHz  data  was  only  measured 
down  to  -100  °C,  and  thus  it  is  still  not  clear  how  the  data  would  compare  at  lower 
temperatures.  There  is  a  slight  difference  between  the  temperature  dependences  of  the  10 
GHz  data  and  the  theoretical  calculations.  This  may  indicate  that  at  10  GHz  all  of  the 
extrinsic  contributions  had  not  completely  "relaxed  out",  or  that  the  actual  temperature 
dependence  was  not  completely  accounted  for  in  the  theory  (.most  of  the  coefficients  of  the 
energy  function  were  assumed  to  be  independent  of  temperature).  In  addition  there  is  some 
experimental  error  in  the  measurements.  The  calibration  for  these  measurements  was  only 
made  at  room  temperature,  and  thus  an  error  of  ±  10%  has  been  estimated  for  the  10  GHz 
data. 

V.  EXTENSION  OF  THE  THEORY  INTO  THE  PZT  SYSTEM 

In  this  section  a  procedure  will  be  presented  that  can  be  used  to  calculate  the 
compositional  dependence  of  the  AG  of  the  antiferroelectric  orthorhombic  phase  into  the  PZT 
system.  This  procedure  was  needed  to  complete  the  development  of  a  thermodynamic  theory 
to  model  all  of  the  phase  transitions  of  the  PZT  system.^*  ^ 


The  AG  of  the  antiferroelectric  orthorhombic  phase  for  lead  zirconate  was  calculated 
in  Section  IV  using  Equation  (34).  The  constants  needed  for  these  calculations  were 
determined  from  the  available  experimental  data,  as  described  in  Section  HI.  Unfortunately, 
similar  experimental  data  is  not  available  for  PZT  compositions  in  the  antiferroelectric  region, 
and  thus  the  compositional  dependence  of  each  of  the  constants  necessary  to  calculate  the  AG 
of  the  Aq  phase  could  not  be  determined.  However,  the  following  method  was  used  to 
calculate  the  AG  of  the  Aq  phase  into  the  PZT  system  from  the  compositional  dependence  of 
a  single  constant,  which  was  determined  by  equating  the  AG's  of  the  Aq  and  FR  phases  at 
the  transition  between  these  phases. 

The  AG  of  the  antiferroelectric  orthorhombic  phase  was  determined  from  the 
following  relation  by  multiplying  the  AG  of  lead  zirconate  (PZ)  by  a  factor  Z: 

AGAq  =  EAGAq(PZ)  <40) 

The  compositional  dependence  of  Z  was  then  determined  by  combining  this  relation  with 
Equation  (35): 

5  =  AGFR(AtTAF)/AGA0(PZ)  (41) 

To  calculate  Z  using  this  equation,  AG  of  the  FR  phase  (high-temperature 
rhombohedral  phase  in  PZT)  at  T^p  was  calculated  from  the  PZT  theory  developed  in 
References  6-10.  The  antiferroelectric-ferroelectric  transition  temperature  T ^p  was 
determined  by  fitting  the  following  polynomial  equation  to  the  experimental  phase  diagram: 1 

T^p  =  220.84  -  706.15  x  -  26778.  x2  (42) 

where  x  is  the  mole  fraction  PbTiC^  in  PZT. 

AG  of  the  Aq  phase  for  lead  zirconate  was  calculated  from  the  equations  and  data 
presented  in  this  paper,  except  that  the  value  of  the  Curie  constant  was  changed.  The 


compositional  dependence  of  the  Cune  constant  that  was  used  in  the  PZT  theory8  resulted  m 
a  value  of  2.0  x  10^  °C,  instead  of  the  value  of  1.5  x  105  °C  that  was  used  in  the  previous 
calculations  in  this  paper.  To  be  consistent  with  the  PZT  calculations  the  value  of  the  Curie 
constant  was  changed,  which  also  caused  the  value  of  the  Q44  constant  to  change  to 
889.27  m*  °C/C^,  instead  of  the  value  listed  in  Table  II. 

Using  the  calculations  described  above  in  Equation  (41 ),  the  E  constant  was 
calculated  versus  composition,  as  plotted  in  Figure  8.  By  then  assuming  the  Z  constant  to  be 
independent  of  temperature,  the  AG  of  the  antiferroelectric  phase  was  calculated  versus 
composition  into  the  PZT  system,  as  shown  in  Reference  10.  This  method  resulted  in 
excellent  agreement  between  the  experimental  and  theoretical  phase  diagrams.1® 

IV.  SUMMARY 

A  two-sublattice  theory  was  used  to  derive  an  energy  function  to  account  for  the 
ferroelectric  and  antiferroelectric  behavior  of  lead  zirconate.  Solutions  and  property  relations 
were  derived  from  this  energy  function  corresponding  to  the  stable  solid  phases  of  lead 
zirconate.  The  coefficients  necessary  to  calculate  the  ferroelectric  rhombohedral  AG,  and 
spontaneous  polarization  and  strain;  and  antiferroelectric  AG  and  spontaneous  strain  were 
determined  from  the  available  experimental  data  in  the  literature. 

Additional  data  were  needed  to  model  the  intrinsic  polycrystalline  dielectric 
properties  of  lead  zirconate.  To  provide  these  data  the  dielectric  properties  were  measured 
down  to  low  temperatures  (-  268.8  °C),  where  extrinsic  contributions  to  the  properties 
"freeze  out".  These  data  were  then  used  to  determine  the  value  of  a  combination  of 
constants,  which  could  be  used  to  calculate  the  intrinsic  polycrystalline  dielectric  constant 
versus  temperature.  The  calculations  indicate  that  at  room  temperature  72%  of  the 
experimentally  measured  polycrystalline  dielectric  constant  (measured  at  1  KHz)  is  due  to  the 
intrinsic  averaging  of  the  single-domain  constants,  while  28%  is  from  extrinsic 


contributions.  The  calculations  were  found  to  be  in  good  agreement  with  experimental  data 
at  high* frequencies  (>  10  GHz),  indicating  that  the  extrinsic  contributions  had  "relaxed  out” 
A  simple  method  was  presented  to  calculate  the  AG  of  the  antiferroelectnc 
orthorhombic  phase  into  the  PZT  system.  This  method  was  used  to  complete  the 
development  of  a  thermodynamic  theory,  which  quantitatively  accounts  for  all  of  the  known 
phase  transitions  in  the  PZT  system.  ^*10 

Additional  data  is  still  needed  to  determine  values  of  the  coefficients  of  the  energy 
function  that  were  not  determined  in  this  paper.  The  effects  of  the  application  of  electric 
field  ^  or  mechanical  stress'-®  may  provide  some  of  this  dam 
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